THE BELL SYSTEM 


echnical Journal 


DEVOTED TO THE SCIENTIFIC AND ENGINEERING 


ASPECTS OF ELECTRICAL COMMUNICATION 





VOLUME XXXxXI JANUARY 1952 NUMBER 1 





The Ferromagnetic Faraday Effect at Microwave Frequencies and 
its Applications—The Microwave Gyrator C. L. HOGAN 1 


Dialing Habits of Telephone Customers a 
CHARLES CLOS AND ROGER I. WILKINSON 32 


Selective Fading of Microwaves 
A. B. CRAWFORD AND W. C. JAKES, JR. 68 


Propagation Studies at Microwave Frequencies by Means of Very 
Short Pulses 0. E. DE LANGE 91 


Properties of Ionic Bombarded Silicon RUSSELL S. OHL 104 


Mechanical Properties of Polymers at Ultrasonic Frequencies 
WARREN P. MASON AND H. J. MCSKIMIN 122 


Relay Armature Rebound Analysis ERIC EDEN SUMNER 172 


Abstracts of Bell System Technical Papers Not Published in This 
Journal 201 


Contributors to This Issue ; 213 





COPYRIGHT 1952 AMERICAN TELEPHONE AND TELEGRAPH COMPANY 


4 
fo ge 
~Se 


on 


’ 


THE BELL SYSTEM TECHNICAL JOURNAL 
PUBLISHED SIX TIMES A YEAR BY THE 
AMERIGAN TELEPHONE AND TELEGRAPH COMPANY 


195 BROADWAY, NEW YORK 7, N.Y. 


CLEO F. CRAIG, President 
CARROLL O. BICKELHAUPT, Secretary 


DONALD R. BELCHER, Treasurer 





EDITORIAL BOARD 


F.R. KAPPEL 0. E. BUCKLEY 
H.S. OSBORNE M. J. KELLY 
J.J. PILLIOD A. B. CLARK 
R. BOWN D. A. QUARLES 


F.J. FEELY 


PHILIP C.JONES, Editor 


M. E. STRIEBY, Managing Editor 





SUBSCRIPTIONS 


Subscriptions are accepted at $3.00 per year. Single copies are 75 cents each. 


The foreign postage is 65 cents per year or 11 cents per copy. 


PRINTED IN U.S.A. 


THE BELL SYSTEM 


TECGCHNICGCAL JOURNAL 





VOLUME XxxI JANUARY 1952. NUMBER 1 





The Ferromagnetic Faraday Effect 
at Microwave Frequencies 
and its Applications 


The Microwave Gyrator 
BY C. L. HOGAN 


A new microwave circuit element dependent on the Faraday rotation of a 
polarized wave has been developed. The element violates the reciprocity 
theorem and, because it shares this property with a gyroscope and because 
at 1s dependent on gyromagnetic resonance absorption, it has been termed 
a microwave gyrator. It 1s a low-loss broadband device with many applica- 
ttons. Among these are one-way transmission systems, microwave circula- 
tors, microwave switches, electrically controlled variable attenuators and 
modulators. 

The microwave gyrator has been realized by making use of the Faraday 
rotatton in preces of ferrite placed in the waveguide. Polder has previously 
shown, in his analysis of the gyromagnetic resonance phenomenon, that 
ferromagnetic substances should show appreciable Faraday rotations at 
microwave frequencies. In the present study, Polder’s analysis has been 
extended to include a wave being propagated through a ferromagnetic sub- 
stance with dielectric and magnetic loss, and data are presented which give 
experimental verification of the theory. In addition an experimental tech- 
nique is described which may be of some interest in studying the properties 
of ferrites at microwave frequencies. 
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Photograph of the experimental setup shown diagrammatically in Fig. 5. 


INTRODUCTION 


In a recent series of articles, Tellegen’ has discussed the possible 
applications of a new circuit element which he calls a gyrator. He defines 
the ideal gyrator, in principle, as a passive four-pole element which is 
described by: (see Fig. 1) 


y= — Ste 2= Str (1) 


Since the coefficients above are of opposite sign, the gyrator violates the 
theorem of reciprocity. Any network composed of the usual electrical 
circuit elements—resistors, inductors, capacitors, and transformers— 
will satisfy the theorem of reciprocity. In simple terms, this theorem 
states that if one inserts a voltage at one point in the network and 
measures the current at some other point, their ratio (calied the transfer 
impedance) will be the same if the positions of voltage and current are 
interchanged. In the gyrator, however, this transfer impedance for one 


L, Le 
— < 





Fig. 1—General four-pole. 
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direction of propagation is the negative of that for the other direction of 
propagation. Essentially this means that a 180° phase difference exists 
between the two directions of propagation. For this reason it has been 
suggested that the element could be more aptly called a directional 
inverter.” 

Network synthesis today is based upon the existence of four basic 
circuit elements: the capacitor, the resistor, the inductor, and the 
ideal transformer. It is apparent that the introduction of a fifth circuit 
element, the gyrator, would lead to considerably improved solutions for 
many network problems. In fact, Tellegen’ has shown that the synthesis 
of resistanceless four-pole networks would be much simplified by its 
introduction. In addition, McMillan’ has shown that it would be possible 
to construct a one-way transmission system if a gyrator were available, 
and Miles* has shown that it would be possible by use of a gyrator to 
construct a network which is equivalent to a Class A vacuum tube 
amplifier circuit. While the realizable power gain of these gyrator circuits 
is necessarily always not greater than unity, many other networks in- 
cluding gyrators are possible which have properties analogous to vacuum 
tube circuits and some of these may be of practical importance. Since 
this new element offers such interesting possibilities in network synthesis, 
a study has recently been made in these Laboratories of possible methods 
for realizing the gyrator. 

A gyrator was employed by Bloch® in his measurement of the mag- 
netic moment of the proton. Bloch made use of the phenomenon that 
if two crossed coils with a mutual core are adjusted so that there is zero 
mutual inductance between them and if a steady magnetic field is 
.applied perpendicular to the axes of both coils, then an ac voltage ap- 
plied to one of the coils will induce a voltage in the second due to the 
gyromagnetic resonance phenomenon. This induction is ordinarily 
extremely small unless the magnetic field is adjusted so that the exciting 
frequency coincides with a gyromagnetic resonance frequency of the 
material which forms the mutual core of the two coils. In Bloch’s ex- 
periment, the magnetic field was held constant and the exciting fre- 
quency was adjusted until it coincided with the gyromagnetic resonance 
frequency of the proton. If they were wound over a paramagnetic or 
ferromagnetic material, the two crossed coils would form a gyrator 
when the magnetic field was adjusted so that the frequency of the 
exciting field coincided with the gyromagnetic resonance of the unpaired 
electrons. The fact that this structure constituted a gyrator was first 
recognized by Tellegen’ and has been discussed by Beljers and Snoek’ 
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in a paper which gives a very satisfying physical model with which to 
interpret gyromagnetic phenomena, occurring within ferrites. 

Physical analysis indicates that the properties of ferromagnetic 
materials can be explained by assuming that the electron behaves as if 
it were a negatively charged sphere which is spinning about its own 
axis with a fixed angular momentum. This rotation of charge imparts 
to the electron a magnetic moment which is a function of the electric 
charge on the electron, the angular velocity of the electron, and its 
size. Thus the electron behaves as if it were a spinning magnetic top, 
whose magnetic moment lies along the axis of rotation, and its behavior 
can be understood by considering a spinning gyroscope suspended in 
gimbal rings at a point not coinciding with its center of gravity. If a 
gyroscope, thus supported in a gravitational field, is lifted away from 
its position of minimum potential energy and then released, it will not 
return to the position of minimum energy but will precess about the 
vertical axis. This is illustrated in Fig. 2 where the spinning gyroscope 
makes an angle 6 with the vertical zg, axis. Its equilibrium motion, in 
the absence of damping, is a precessional motion about the vertical 
axis with a velocity wp. 

If the gyroscope be regarded as initially hanging vertically downward 


Zz 






Xx Xx 


Fig. 2 (left)—Precessional motion of a gyroscopic pendulum in a gravitational 
field. Fig. 3 (right)—Precessional motion of a gyroscopic pendulum in a gravita- 
tional ficld which oscillates between the directions A and B. 
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as indicated in Fig. 3 and then a gravitational force is suddenly made to 
act along the y axis so that the net gravitational force acts along A, 
it is obvious that the gyroscope will begin to precess about the gravita- 
tional field direction as indicated by the small dotted circle. However, if 
after completing a half cycle, the horizontal component of the gravita- 
tional field is reversed so that now the net gravitational field acts along 
the vector B, the gyroscope will begin to precess about B as indicated 
by the intermediate size dotted circle. If the horizontal component of 
the gravitational field is again reversed after the gyroscope completes 
another half cycle in its precession, the gyroscope will again begin to 
precess about the direction A and the actual path of the precessional 
motion will be along the path a-b-c-d. If this process is continued in- 
definitely, the gyroscope will precess in larger and larger circles around 
the vertical until the damping becomes large enough to contain the 
gyroscope in some equilibrium circle (assuming that the damping is 
large enough to accomplish this). 

The above model affords a classical picture which can be used quite 
readily to describe the motion of the electrons in a ferrite. If the ferrite 
is initially saturated along the z axis by a steady magnetic field, the 
electrons will come to rest with their magnetic moments lying along 
the ¢ axis, as the gyroscope in Fig. 3. If now an alternating magnetic 
field is applied along the y axis, the electrons will begin to precess in 
larger and larger circles about the g axis until they finally reach some 
equilibrium position under the influence of the magnetic fields and the 
damping. Thus it is apparent in the gyromagnetic resonance experi- 
ments described above why an alternating field applied perpendicular 
to a steady magnetic field in a ferrite will give rise to a varying flux 
perpendicular to both the steady field and the alternating field. It is 
also apparent why the alternating flux along the x axis is 90° out of 
phase with the alternating flux along the y axis. Since precession of the 
top will always be in the same direction regardless of whether the alter- 
nating field is applied along the x or y axes, consideration of Fig. 3 
makes it apparent how the two crossed coils with ferrite at their center 
can constitute a gyrator which violates the reciprocity relation in a 
manner described by Equations (1). To the present time, however, no 
practical circuit element making use of this phenomenon has been 
constructed because the coefficient of coupling between the coils is 
always small, even in the vicinity of the resonant frequency, and also 
because the losses in the materials available are so high in the vicinity 
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of resonance that the insertion loss of such a device would be prohibitively 
large. 

MeMillan in his original article,’ showed that a gyrator could be 
realized by means of mechanically coupled piezo-electric and electro- 
magnetic transducers. Later, McMillan’ pointed out that a gyrator 
could be realized by means of the Hall effect in a square plate of bismuth, 
as was also predicted by Casimir.’ Another similar possibility would be 
an electrical-electrical coupling through a gyroscopic link. A gyrator has 
been built by W. P. Mason of these Laboratories which makes use 
of the Hall effect in a crystal of germanium.® This gyrator showed an 
insertion loss somewhat higher than the theoretical loss of 12.3 db. 
R. O. Grisdale of these Laboratories suggested that these losses could 
be greatly reduced if the same Hall effect principle were applied to a 
vacuum tube which contained four electrodes which could both emit 
and collect electrons. This device is no longer passive, but such a struc- 
ture has been built and showed an insertion loss of about 7 db, only 
slightly higher than the theoretical loss which would be expected from 
this geometry. 

In view of the substantial losses found to exist in the earlier forms of 
gyrator discussed above, a study of other “anti-reciprocal’’ phenomena 
which might lead to the realization of a relatively low loss gyrator was 
undertaken. . 

It has long been known that the Faraday rotation of the plane of 
polarization in optics is anti-reciprocal. In order to observe the Faraday 
rotation, polarized electromagnetic waves must be transmitted through 
a transparent isotropic medium parallel to the direction of the lines of 
force of a magnetic field. The effect is usually produced by placing the 
material along the axis of a solenoid. The rotation is “positive” if it is 
in the direction of the positive electric current which produces the field 
and “negative” if in the opposite direction. All optically transparent 
substances show the Faraday rotation. 

Its anti-reciprocal property distinguishes the Faraday effect from 
optical rotations caused by birefringent crystals, or by the Cotton- 
Mouton effect, which are reciprocal. That is, if a plane polarized light- 
wave is incident upon a birefringent crystal in such a manner that the 
plane of polarization is rotated through an angle # in passing through 
the crystal, then this rotation will be cancelled if the wave is reflected 
back through the crystal to its source. In the Faraday rotation, however, 
the angle of rotation is doubled if the wave is reflected back along its 
path. Hence, if the length of path through the “active” material is 
adjusted so as to give a 90° original rotation, the beam on being reflected 
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will have its plane of polarization rotated a total of 180° in passing in 
both directions through the material. Thus, the Faraday rotation in 
optics affords an anti-reciprocal relation quite analogous to the anti- 
reciprocal property of the gyrator postulated by Tellegen. 

Lord Rayleigh’ described a one-way transmission system in optics 
which makes use of the Faraday rotation. Lord Rayleigh’s “‘one-way”’ 
system consisted of two polarizing Nicol prisms (oriented so that their 
planes of acceptance made an angle of 45° with each other), with the 
material causing the Faraday rotation placed between them. Thus, 
light which was passed by the first crystal and whose plane of polariza- 
tion was rotated 45° would be passed by the second crystal also. But, 
in the reverse direction, the rotation would be in such a sense that light 
which was admitted to the system by the second crystal would not be 
passed by the first. 

Although Rayleigh’s one-way transmission system can be actually 
realized, it is experimentally difficult since most substances show ex- 
tremely small Faraday rotations. In fact, large rotations for transparent 
substances in the optical region are of the order of one degree per cm 
path length for an applied magnetic field of 1000 oersteds. To realize 
a rotation of 45° would require maintaining a field of 1000 oersteds over 
a distance of approximately one-half meter. The Faraday effect in 
ferromagnetic substances, however, is unique in that it shows rotations 
many orders of magnitude greater than the rotations exhibited by any 
other substances. For instance, Kénig’’ reports rotations of 382,000°/em 
by passing light through thin layers of magnetized iron. These data, of 
necessity, however, were taken on extremely thin sections and the total 
rotation obtained for any specimen did not exceed 10°. In order to 
obtain appreciable rotations in a device of practical size, it 1s necessary 
to obtain a material which shows a rotation at least intermediate be- 
tween those reported for iron and other ordinary materials. In addition, 
in order to make effective use of these rotations, the material must be 
transparent to the radiation which is being used. 


THEORY OF THE FERROMAGNETIC FARADAY EFFECT 


Polder” has shown in his analysis of the ferromagnetic resonance 
phenomenon, that a plane electromagnetic wave at microwave fre- 
quencies should show appreciable Faraday rotation when propagated 
through a ferromagnetic material which is magnetized in a direction 
parallel to the direction of propagation of the wave. Polder has neglected 
both magnetic and dielectric losses in his analysis and although for the 
ferrites which are of greatest interest, this approximation is quite 
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valid, nevertheless the more complete theory is developed below. The 
exact theory of this phenomenon should, of course, be approached 
through quantum mechanics, but since the classical theory, in this 
particular case, gives a result as satisfactory as the quantum theory 
and since it lends itself more aptly to a fundamental physical inter- 
pretation of the phenomenon, it is the classical theory which is developed 
here. Quantum mechanically, Faraday rotations in the optical region are 
accounted for by the Zeeman splitting of the spectral lines. 

The classical model which proves quite adequate for the description 
of ferromagnetic resonance is that illustrated in Figs. 2 and 3 which 
regards the electrons of the material which contribute to the magnetism 
as being spinning magnetic tops. The angular momentum of each 
electron is: 





[J] = 3(h/2m) (2) 
J = Angular momentum of electron (gm em’/sec) 
h = Planck’s constant (6.62 X 10°” erg sec) 
The magnetic moment which arises due to this rotation is: 
eh 
|u| a 4arme (3) 


where: 
us = Magnetic moment of electron (Bohr magneton) 
e = Charge on electron (4.80 X 10° E.S.U.) 
m = Mass of electron (9.10 X 10°” gm) 
c = Velocity of light (8 X 10°° cm/sec) 


The so-called gyromagnetic ratio of the electron is the ratio of these 
quantities and is given by: 
é lve 
Bae OR See as ME 4 
D2 one a 4) 
If a steady magnetic field is applied to the sample such that the elec- 
tron sees an effective field H, then a torque will be applied to the electron 
which tries to turn the electron so that its magnetic moment lies along 
the field direction. However, as indicated in Fig. 2, the electron will 
precess around the field direction until damping forces dissipate the 
energy of precession. The equation of motion of the electron is: 


2, FU ~1 dus 


us X H =~, Y dt (5) 
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The equation of motion of the magnetization per unit volume can 
thus be written: 


= yM XH (6) 
where: 
M = Magnetization of medium 
H = Macroscopic internal magnetic field 


The above equation, however, does not include damping. The damping 
force, regardless of its origin, must be so introduced into the above 
equation that it tends to cause the electron’s axis of rotation to line up 
with the field direction. It has been shown by Yager, Galt, Merritt 
and Wood” that the shape of the resonance absorption line can be 
accounted for if the damping term is introduced in the following way: 

dM ae gees ya = Tee ee 
aq YM XH par] We x (it x HDI (7) 

The vector M X (M X H) is simply a vector which is in the proper 
direction to act as a damping force (torque) and the coefficient is chosen 
so as to give the correct units along with the parameter, a, which must 
be determined experimentally and which gives the magnitude of the 
damping torque. 

Equation 7 then is the equation of motion of the magnetization of 
an arbitrarily shaped body under the action of an arbitrary internal 
field, H. In the appendix, it is shown that if a steady magnetic field, 
H, , is applied along the zg axis and then a small alternating field is ap- 
plied in an arbitrary direction to a sample which is infinite in size, the 
equation relating the resulting alternating flux density, b, and the ap- 
plied alternating field, h is: 


be = ph, — jKhy 


by =jKhs + phy (8) 
b, = h, 
where 
w= pl — jy” (9) 
K = K' — jk" (10) 


Equations which give » and K in terms of the applied magnetic field 
and fundamental atomic constants are given in the appendix. 
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Equations (8) are easily interpreted in terms of the spinning gyroscope 
model of Fig. 3. If magnetic losses had been ignored (i.e. a = 0) then 
both » and K would have been real. Under this condition, it is seen that 
if an alternating field, h, , is applied along the y axis, then an alternating 
flux, b, , is created along the y axis which is in phase with h, , and an 
alternating flux, b, , is created which is 90° out of phase with h, . Reci- 
procity between the x and y directions would demand that both terms 
containing 7K should have the same sign. Thus, Equations (8) give a 
quantitative expression for the results which were previously qualita- 
tively deduced by means of the electronic model illustrated in Figs. 2 
and 3. . 

If a waveguide is filled with a ferromagnetic material such as a ferrite 
and if then a steady magnetic field is applied along the axis of the 
waveguide, it is necessary in order to describe this wave to find a solu- 
tion to Maxwell’s equations which is consistent with Equations (8) 
and in which b, h, # and D are all proportional to exp [jwt — Ts]. This 
problem is not solved exactly. However, in the appendix a solution is 
obtained for an infinite plane wave. It is found that the ferromagnetic 
medium can support only a positive or a negative* circularly polarized 
wave or a combination of both. It is also shown in the appendix that the 
propagation constants for these two circularly polarized waves are dif- 
ferent and are given by the following expressions: 


ry =" Vu + Bel (11) 
and 
r=" Va — Kel (12) 
where ~ 
T's = Propagation constant 
w = Angular frequency of wave 


c = Velocity of light in unbounded space (3 X 10°° cm/sec) 
e¢ = Complex dielectric constant of medium 
In Equations (11) and (12) it is apparent that the effective perme- 
ability of the medium to a positive circularly polarized wave, for in- 
- * The usual notation is used here, where the positive component is the com- 


ponent which is rotating in the direction of the positive electric current which 
creates the steady longitudinal field. 
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stance, is given by the expression (u + K), and not by the usual perme- 
ability, bz/hz = yw. It is also apparent that the quantity u» + K can 
vary over wide limits in the vicinity of the ferromagnetic resonance. 
For this reason, care must be taken in interpreting permeability data 
for ferromagnetic materials which now occur in the literature and which 
were obtained by means of impedance measurements at microwave 
frequencies, since the above equations indicate that this method does 
not measure the same quantity that is measured at low frequencies by 
means of a toroidal sample overwound with two coils. The low fre- 
quency measurement of permeability obviously measures the quantity 
which is designated as » in Equation (8). 

If Equations (11) and (12) are solved for the attenuation constants, 
a,, and the phase constants, 63, the following results are obtained: 


a fe + Ke! 
ee eee =e 


; H{4/ (1 + tan 6,[4 tan 62 + tan 6n(1 + tan? 6z)] + tan? i} 
} (13) 





- — 1 — tan 6, tan ,| 


and 
poo ,/C ERE 
= Cc 2 


4 H4/a + tan 6n[4 tan 6g + tan 6,(1 + tan? 6,)[ + tan ” i} 





, (14) 
- + 1-+ tan 6, tan | 
where: 
uv uv 
we tk 
tan dn = ———>=> 
an aaa % 


(The + sign must be used for a positive circularly polarized wave; 
the negative sign for the negative circularly polarized wave.) 
oe” 
oF 


tan 6a = = dielectric loss tangent 


lop] 


e = ef — je” = complex dielectric constant 
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It is almost impossible to get a feeling for what these equations mean 
with respect to a wave travelling through the medium, especially since 
» and K are given by equations which are almost as difficult to per- 
ceive. An appreciation of these equations can be obtained however, by 
reference to Fig. 4 which gives qualitatively the behavior predicted by 
these expressions. Essentially, a and 8 are functions of two variables. 
These are w, the frequency of the wave, and H,, the applied magnetic. 
field. In Fig. 4, the index of refraction and attenuation of the positive 
circularly polarized component are given relative to these values for 


INDEX OF REFRACTION ABSORPTION OF 
OF POSITIVE COMPONENT.) _-POSITIVE COMPONENT 


(ARBITRARY UNITS) 


fF 
Ha = 3. ~(N-Nz) Mz 


FOR SAMPLE 
WITH CYLINDRICAL SYMMETRY 


ABSORPTION & INDEX OF REFRACTION 





APPLIED MAGNETIC FIELD (ARBITRARY UNITS) 


Fig. 4—Index of refraction and absorption of a positive circularly polarized 
wave relative to the same quantities for a negative circularly polarized wave 
being propagated through a magnetized medium. 


the negative component. Hence, both the index of refraction and 
attenuation of the negative component are represented by the abscissa 
of the graph. In Fig. 4 these quantities are plotted as a function of the 
applied magnetic field for a wave of a fixed frequency. Many of the 
properties of the medium are clearly displayed in this graph. In partic- 
ular, as the field necessary for ferromagnetic resonance is approached, 
the attenuation of the positive component becomes larger and larger. 
Eventually this component will be substantially completely absorbed 
and only the negative circularly polarized component will be propagated. 
Hence it should be possible to establish a circularly polarized wave in a 
waveguide simply by passing the dominant mode through a ferromag- 
netic material which is subjected to a longitudinal magnetic field of the 
proper amplitude. However, there will be an absorption of one-half 
of the power being propagated. If Fig. 4 had been plotted as a function 
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of the frequency of the wave for a fixed magnetic field, a similar set of 
curves would have resulted. This set would indicate the frequency 
dependence of the Faraday rotation. If the frequency of the wave is far 
removed from the resonance frequency, the difference between the 
indices of refraction of the positive and negative component is not 
frequency dependent. However, near resonance, this difference is a 
very rapidly varying function of the frequency. It is to be remembered 
that these equations were derived for an infinite plane wave. However, 
it would be expected that these equations would describe quite accurately 
the propagation of the dominant mode in a waveguide. The approxima- 
tion would, of course, be better when the cut-off wavelength was much 
greater than the unbounded wavelength. This condition is met when 
the waveguide is filled with ferrite and for these cases quantitative 
agreement is obtained. 

The above analysis shows that if a dominant mode wave (plane 
polarized) is incident upon a ferromagnetic material which is magnetized 
along the length of the waveguide, the wave will split into positive and 
negative circularly polarized waves whose phase constants are given by 
Equation (14).) Since the two circular components travel with different 
velocities in the medium, they will upon emerging from it unite to form 
a plane polarized wave whose plane of polarization has been rotated with 
respect to the incident polarization. The angle of rotation of the polariza- 
tion is given by: 


§ = = [B_ — B,] (15) 


bor 


where: 
£ = path length through ferromagnetic material (cm) 


Tn order to evaluate Equation (15), it must be combined with Equation 
(14). However, a few approximations are valid in Equation (14) which 
make it much simpler. In particular many ferrites exist for which the 
magnetic losses are extremely small as long as the internal field within 
the body is kept small so that the frequency of the wave does not ap- 
proach the ferromagnetic resonance frequency. This field can be kept 
small if the magnetic field is not raised above the point necessary to 
saturate the ferrite. Kittel has shown that for a finite body the effective 
internal magnetic field that determines the resonant frequency is given 
by: 


H; = [Ha oe (N; a N.)M,\[Ha st (Ny a N.)M,] 
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where H, is in oersteds. The ferromagnetic resonance frequency is 
given by: 


fres = 2.8H. \/1 + a2 megacycles (16) 


It is easily shown that the following formula is approximately valid 
for a ferromagnetic body with a circular or square cross-section 


_ 4r + Nw — 1) 


~ Ae + Nu — 1) = ou 


H, 


at saturation. Where: 
ww = true de permeability at saturation 
_ N = demagnetizing factor in x and y directions. 


_ If an average value of 1000 is assumed for the de permeability, then 
H,.can be readily computed for various shapes. 
For a®thin disc: 


N = 0 Nz = 4a 
and 
Hi, 
He = 7500 
If a thin disc saturates at 1500 gauss, then: 
H, = 1.5 oersteds 


froes & 4.2 megacycles 


For a long thin pencil: 


and 
H, = 1000H. 
For this case the body could be saturated with a field of about1.5 
oersteds, so: 
H, = 1500 oersteds 
and 


fros = 4200 megacycles 


If, for this case, the resonance frequency is so close to the operating 
frequency that losses due to ferromagnetic resonance become pro- 
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hibitive, it is wise to then raise the applied field to some high value, 
so that the resonance frequency will fall well above the operating fre- 
quency. Thus, for many cases of interest it is possible by various means 
to place the ferromagnetic resonance absorption frequency sufficiently 
far from the operating frequency so that magnetic losses due to this 
phenomenon are negligible.* The data accumulated to date indicate 
that the major component of the magnetic losses at microwave fre- 
quencies is due to this phenomenon. Only in a few cases have data been 
taken which have indicated that other factors, such as domain wall 
relaxation, contribute to the magnetic loss at microwave frequencies. 
If then, the magnetic field is controlled so that the ferromagnetic 
resonance absorption is negligible, Equations (13) and (14) can be 


simplified to: 
a4 eee 4/ 4/1 + tants — 1 ~ (18) 
c 2 i 
pe = 2 4/E SDI VF tame +d (19) 


which can be written as: 


at 4/E Wwak - (20) 








and: 


and | 
-ty/elt! ror @1) 


where y’ and K’ are given in the appendix. 

If Equation (21) is now inserted into Equation (15) a formula for 
rotation is obtained which is valid within the limits of the above ap- 
proximations. If in addition, the frequency of the wave is sufficiently 
greater than the resonance frequency, so that: 


Wres Kw (22) 
then eee (15) takes the particularly simple form: 


y bel Fel 4/1 4 Mele 4/1 = 
15 1+ sae (23) 


Most ferrites saturate at 2,000 gauss or less. Hence, for a frequency of 





* This is not always possible, for some ferrites, in the polycrystalline state, 
exhibit extremely broad ferromagnetic resonance absorption lines and it is diffi. 
cult to operate at any frequency without appreciable absorption. 
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9,000 megacycles, 


4nMy — 2000 X 17.6 X 10° 
w ~ 9000 X 27 X 108 


Hence, the following approximation will be valid to within 5 per cent. 


/ 1 aM 1 Ge 
a) 2 w 


With this approximation, Equation (23) reduces to: 


OF cl le| +e 
(=a / 5 [47M 7] (25) 


Equation (25) is quite remarkable. Not only does it predict large 
rotations, but it also predicts that, within the above approximations 
the rotation will not depend upon the frequency of the incident radia- 
tion. For the assumed values, 





= 0.622 (24) 





e = 15 
e” = 0 
4nM, = 1000, 


Equation (25) predicts rotations of, 


; = 65°/cm. 


DESCRIPTION OF EQUIPMENT AND MEASURING TECHNIQUES 


The Faraday rotation has been measured in a large number of ferrites 
in order to verify the above theory and in an effort to improve the 
characteristics of the microwave gyrator. A diagram of the experimental 
equipment is given in Fig. 5, and a diagram of the test chamber in which 
the rotations were measured is given in Fig. 6. In the test chamber, 
two rectangular waveguides are separated by a circular waveguide, 
the proper nonreflective transitions being made at each end of the 
circular section, which is about twelve inches long. One rectangular guide 
is supported so that it can be rotated about the longitudinal axis of the 
system. The dominant 7 mode is excited in one rectangular guide, 
and by means of the smooth transition this goes over into the dominant 
TE, mode in the circular guide. The rectangular guide on the opposite 
end will accept only that component of the polarization which coincides 
with the TE mode in that guide, the other component being reflected 
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at the transition. Absorbing vanes, inserted in the circular section, 
absorb this reflected component. The circular guide is placed in a solenoid 
to establish an axial magnetic field along its length. 

The ferrite cylinders to be measured were placed at the mid-section 
of the circular guide. When a cylinder was used which did not fill the 
cross-section of the guide, it was supported along the axis of the guide 
by means of a hollow polystyrene cylinder which did fill the guide. 

In addition to measuring the Faraday rotation, measurements of 
insertion loss were made by determining the power transmitted under 
identical conditions with the ferrite cylinder removed, and the ellipticity 
of the transmitted wave was determined by measuring the power trans- 
mitted when the rectangular guide on the detector side was rotated to’ 
both positions of maximum and minimum transmission. Power trans- 
mission measurements could be repeated within 0.2 db. Measurements 
of the angle of rotation of the plane of polarization could be repeated 
within 4° except in the region close to the gyromagnetic resonance where 
rotations were large and ellipticity so great that it was difficult to decide 
the positions of maximum and minimum transmission. These errors 
increased up to the point where the transmitted wave was circularly 
polarized where it was impossible to measure the angle of rotation. 


EXPERIMENTAL RESULTS 


Equation (25) indicates that the rotation per unit path length through 
the ferromagnetic material is proportional to the magnetization of the 
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Fig. 5—Experimental equipment set-up used to measure Faraday rotations. 
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sample and is not dependent directly on the applied magnetic field. 
Fig. 7 shows the dependence of rotation upon magnetization for a 
sample of manganese zinc ferrite, and indicates that after the sample is 
saturated, the rotation is sensibly independent of the applied magnetic 
field. In addition, the complex dielectric constant and the saturation 
magnetization of this sample were measured. From these the rotation 
per centimeter path can be computed from the above theory using 
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Fig. 6—Detail of test chamber in which rotations were measured. 


FARADAY ROTATION IN DEGREES 
PER CENTIMETER 





MAGNETIC FIELD IN OERSTEDS 
Fig. 7—Angle of rotation versus applied magnetic field for a thin disc of man- 
ganese zinc ferrite. : 
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Equation (25). For a particular sample of manganese zine ferrite, the 
following measurements were made: 


(Sample No. 1) 
e = 17 
" = 24 . 
4nrM.x = 1500 gauss 
Using this data, equation (25) predicts: 


m 
I 


— = 121.2°/cm. 


It is seen in Fig. 6, that the actual measured rotation at saturation is 
approximately 123°/em. Hence an extremely good agreement with 
theory has been obtained for this particular sample. 

Equation (25) also indicates that the rotation per unit path length 
should be sensibly independent of frequency within the above ap- 
proximations. The data are shown in Fig. 8. However, it will be noticed 
that the frequency difference between these two sets of data is relatively 
small (3 per cent), and the cumulative experimental error in measuring 
angles is such that it is difficult to state that the rotation is closer 
than 1° between the two sets of data. This represents a possible differ- 
ence of 5 per cent in the rotation for a change of 3 per cent in the fre- 
quency. Thus, even though these preliminary data support Equation 
(25), it cannot be accepted as conclusive evidence until more measure- 
ments can be made over a wider band width. 
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Fig. 8—Dependence of Faraday rotation upon frequency. 


20 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1952 


The loss characteristics of different ferrites as a function of the ap- 
plied magnetic field differed distinctly from each other. Some ferrites, 
such as manganese zinc ferrite showed extremely high loss which was 
associated with the imaginary part of the dielectric constant. This loss 
was not affected by the application of a magnetic field but remained 
substantially constant as the field was applied. However, as the field 
approached that necessary for ferromagnetic resonance, the total power 
absorbed by the ferrite increased, since the positive circularly polarized 
component was almost completely absorbed by the sample. In fact by 
measuring the ellipticity of the transmitted wave, it is possible to com- 
pute the difference between the absorption of the positive and negative 
circularly polarized components. This has been done for Sample No. 1 
and the result is indicated in Fig. 9. If the curve were continued to 
higher fields, it would represent the shape of the ferromagnetic resonance 
absorption line. 

Some ferrites, such as Ferramic G, showed an almost zero dielectric 
loss but on the other hand caused an extremely large absorption at 9000 
megacycles due to magnetic losses. The major contributions to magnetic 
loss at this frequency should be either losses associated with a domain 
wall relaxation or ferromagnetic resonance absorption due to anisotropy 
fields. Unequivocal data can be obtained by the above techniques to 
identify which loss is predominant. If the loss were due to domain wall 
relaxation (or resonance) it would absorb both the negative and positive 
circularly polarized components equally. Thus as the magnetic field 
was applied and as the ferrite became saturated, the losses in both 
components should decrease as the domain walls disappeared. However, 
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Fig. 9—Ferromagnetic resonance absorption curve determined by measuring 
the ellipticity of a wave transmitted through a cylinder of ferrite in a waveguide. 
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Fig. 10—Rotation of the plane of polarization and absorption of the positive 
and negative circularly polarized components when an electromagnetic wave 
(X-band) is propagated through a cylinder of Ferramic G. 
if instead the loss is associated with ferromagnetic resonance absorption, 
the absorption of the positive component should begin to increase as 
soon as more domains are lined up in a direction where they can absorb 
the positive component. Thus even before the sample is saturated the 
absorption of the positive component should be much larger than the 
absorption of the negative component. Of course, in a polycrystalline 
sample with a large anisotropy both components can be absorbed by 
ferromagnetic resonance absorption when the sample is not completely 
saturated, since the random orientation of the domains which occurred in 
zero field has not been completely eliminated until complete saturation 
occurs. Fig. 10 illustrates the rotation per cm path length versus applied 
magnetic field for a sample of Ferramic G. Superimposed on the same 
figure are curves showing the absorption of the negative and positive 
circularly polarized components. It will be noticed that as soon as the 
sample is saturated, the sample becomes transparent to the negative 
component but almost completely absorbs the positive component. 
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Hence the transmitted wave at this point is almost completely circularly . 

polarized, even though the applied magnetic field would indicate that 

the resonance absorption frequency was far removed from 9000 mc. 
Table I gives the data taken on several ferrites at 9000 me. 


APPLICATIONS OF THE FERROMAGNETIC FARADAY EFFECT—THE MICRO- 
WAVE GYRATOR 


As pointed out in the introduction, the Faraday rotation affords an 
anti-reciprocal phenomenon from which a microwave gyrator can be 
constructed. Such a gyrator is illustrated in Fig. 11 along with diagrams 
which help explain its action. Beneath the gyrator are construction lines 
which indicate the plane of polarization of a wave as it travels through 
the gyrator in either direction. On each diagram is a dotted sine wave 
which is for reference purpose only and indicates the constant plane of 
polarization of an unrotated wave. It is noticed that for propagation 
from left to right in Fig. 11, the screw rotation introduced by the twisted 
rectangular guide adds to the 90° rotation given to the wave by the 
ferrite element making a total rotation of 180°. For a wave travelling 
in the reverse direction, these two rotations cancel each other, producing 
a net zero rotation through the complete element. The unique property 
of the Faraday rotation becomes immediately apparent from this dia- 
gram. In the case of the rotation induced by the twisted rectangular 
guide, the wave rotates in one direction in going from left to right through 
the twisted section, and rotates in the opposite direction when it trans- 
verses the section from right to left. For the case of the rotation induced 
by the ferrite element, the direction of rotation is indicated by the 
arrow in the upper figure for either direction of propagation. The im- 
portant characteristic of the element is the time phase relation between 
two points such as A and B in the upper diagram. It is seen with the help 
of the diagrams illustrating the rotating waves that the field variations 
are in phase at points A and B for propagation from left to right, and 
they are 180° out of phase for propagation from right to left. In other 
words the transmission line is an integral number of wavelengths long 
‘between A and B for propagation from left to right and is an odd integral 
number of half wavelengths long for propagation from right to left. 

From the above description of the properties of the gyrator, many 
of its applications in microwave technology become immediately ap- 
parent. Before discussing these applications in more detail, however, 
it is advantageous to introduce standardized terminology and circuit 
symbols which apply to the gyrator and to other circuit elements 
derivable from it. 
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TABLE I 





Material Dimension (cm) 


Insertion Loss (db) 
Line (db) 


Sample Number 
Applied Magnetic 
Field (oersteds) 
Rotation/em 
Ellipticity* (db) 
SWR on Input 








—_ 


BTL 0.447 X 2.28 0 
Mn;Zn,_;Fe.0, (length x dia.) 245 








735) 58.2 


1225] 107 10. 
1470} 120 |10 10. 
1715) 125 {11 
1960] 123 {11.2 
2206) 121 11.3 
2450) 123 {11.4 
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0 





2695} — 12. 
2940) — 
3185} — 
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2 BTL 0}. 0 
Ni;Zn, _3Fe.0, 1.36 x 2.28 245} 25 

490} 44 
735| 56 
980) 61 
1225} 68 
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2450] 118 


3 Ferramic A 2.54 X 0.635 0 0 
245) 34. 
490} 43. 
735} 48. 
980} 51. 

1225] 54 

1715) 57. 

1960} 60. 

2450} 63. 

2695) 64. 


4 Ferramic G 1.77 X 2.28 0} 60 
245) 38 
490| 77 
735| 124 
980) 157 
1225) 170 
1470} 180 

3430| ¢.p. 
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* Data given is the difference in db between the major and minor components 
of the elliptically polarized transmitted wave. 
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Fig. 11—The microwave gyrator with diagrams which help to explain its 
operation. 
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The “active” elemient of the device, the ferrite cylinder, has been 
termed a ‘Faraday Plate.” 

As was pointed out earlier, the fundamental property of the gyrator 
is the 180° phase difference introduced between the two directions of 
propagation through it. Thus the gyrator may be thought of as a four 
terminal circuit element having no phase shift for one direction of trans- 
mission, and having a 180° phase shift for the opposite direction of 
transmission. A convenient circuit symbol for the gyrator, which indi- 
cates this property, is shown in Fig. 12. 

If the rectangular waveguides on each side of the Faraday Plate are 
rotated about their common axis so as to make an angle of 45° with 


— 
an T =e 


Fig. 12—Circuit symbol for gyrator. 
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Fig. 13—Schematic diagram of polarization circulator. 


each other, then a one-way transmission system can be created which 
is similar to Lord Rayleigh’s one-way transmission system of optics, but 
with the important difference that this one-way transmission system 
does not depend upon frequency but is broad band. This one-way trans- 
mission system can be used, for example, to isolate the generator or 
detector from the waveguide in microwave systems. In this application 
it has the great advantage over the attenuators which are presently 
used for this purpose in that it can be made practically lossless for the 
direction of propagation which is desired but the reflected wave will be 
completely absorbed and hence more complete isolation can be effected. 

A more complex and more useful circuit element, than this simple 
one-way transmission property would at first indicate, is obtained by 
adding a second connection on each side of the 45° Faraday Plate. It is 
suggested that this device be called a polarization circulator. Thus, the 
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polarization circulator actually has four output branches corresponding 
to the two different polarizations at each end. The polarizations of the 
four output branches are indicated in Fig. 13. It is noticed that power 
sent into the polarization circulator with polarization a is turned into 
polarization b,' also b is turned into ¢, ¢ is turned into d, and d is turned 
into minus a. This property is indicated very clearly by the circuit sym- 
bol suggested in Fig. 13, the phase inversion between arms d and a being 
indicated by the minus sign between the d and a arms. 

Another one-way transmission system can be created by combining 
the gyrator with two-normal hybrids. This combination is indicated in 
Fig. 14. Since this device has all of the fundamental properties of the 
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Fig. 14—Schematic diagram of circulator. 


polarization circulator with the exception of the phase inversion be- 
tween arms d and a it is suggested that it be called a “circulator” and 
the circuit symbol suggested which indicates its properties is also given 
in Fig. 14. 

This list of applications is obviously not complete since it includes 
only the fundamental elements from which innumerable specific applica- 
tions can be made. 

In addition to the applications discussed above, which depend upon 
the anti-reciprocal property of the element for their operation there 
are several simple applications which are based only upon the fact that 
the amount of rotation can be controlled externally by adjusting the 
magnetic field. Among these uses are electrically controlled attenuators, 
modulators, and microwave switches. 
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APPENDIX 
The equation of motion of the magnetization of a ferromagnetic 
material is: 
— = 7y(M X A) -— — [M MXH 1 
ap VM XD — AF Lt x GE x HD) (1) 
where 
H = internal magnetic field (oersteds) 


4nM = magnetization of medium (gauss) 


a = parameter which measures the magnitude of the damping force 
on the precessing dipole moment of the sample 


y = gyromagnetic ratio of the electron (y = ge/2mc where g is the 
Landé g factor for the electron). 


If a ferromagnetic material is subjected to a steady magnetic field, 
H,, along the z axis and if then an alternating field is applied in an 
arbitrary direction, Equation (1) must be solved in order to find the 
behavior of the magnetization of the material. To solve this problem, 
the following notation is introduced: 


4rM, = magnetization of medium in absence of alternating field 
H, = externally applied steady magnetic field (oersteds) 

hz, hy, hz = components of applied alternating magnetic field 
Mz, My, mz = alternating components of magnetization 

hi, hi , Hi = components of internal magnetic field 

hi = he — Namz 

hi = hy — Nyn, 

Hi = Ha +h. — NAM, + m.) 

N bs N,, Nz = demagnetizing factors of body. 
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Hence the magnetic field, H, occurring in Equation (1) is defined by: 
Pee fa 


and 


M — tet + mi + (M, + mk 

In solving Equation (1), an exponential, exp [jwt], time dependence 
is assumed for the alternating magnetic field and magnetization, and 
if the following assumption is made: 


hz, hy, he K Ha 


it is easily shown that the alternating components of the magnetization 
of the medium are given by (neglecting terms of the second order in 
small quantities) : 


_ WMGHL + @) + jvoMolh: — jyMeoh, 





MN, = YH (1 oe a’) = ar ue jQwyoH | 5 
c fH? (1 + a’) — w* + jl2wyalli] 
m, = 0 
where: 
j=v-l1 
Since 
b = hi + 4am, (3) 


it is possible by means of Equations (2) and (8) to find the relation 
between the alternating flux density b and the internal alternating field 
h’. If the ferromagnetic body is considered as being infinite, the internal 
fields and applied fields are equal. Hence, for this case: 


by Sungai 
by = 7Khz + ph, (4) 
bz = hz 


where: 
b= p! oz ju” 
= K’ — jK” 
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and: 
Poj4 Aa + @) — w'[4eM yA + o)] + Sa w'y'o’H, 
" [yHa(L + a?) — wf? + 40720? H2 
a 4nM yoly’Ha(l + 0”) — «| 
K! = r 2772/4.) .2). 22 ~«4)~CO «iz 
ly Ao ss a’) = w| + 40"7"o°H, 
K" _ 87M .wy’aHa 
~ [yHG. + a”) — «+ 40y*a" He 


» _ _4nMyooly Hal + @’) + a] 
(ER + a) — of + dota 
In order to find the behavior of a wave being propagated in this 
medium, it is necessary to find a solution to Maxwell’s equations which 
are consistent with the above set of equations and in which, ), h, E, 
and D are of the following form: 
b bo exp [ jeot = T(n-7)] 
ho exp [jot — I'(n-7)] 
= EF, exp [jot — P-7)] 
D = Dy exp [jot — T'(n-r)] 
where Hp, and ho are complex vector functions of the coordinates and 


which satisfy the boundary conditions imposed by the waveguide. 
Further: Ps 


os 


nm = unit vector in the direction of propagation 


ll 


(5) 


Rl solo 


r 


propagation constant 


Maxwell’s equations are: 


_ 18 
VRE ae, 
1 aD 
Ve 2 ae 


Inserting the values given in Equations (5), these become: 


VX Hy — TG x B) = = | (7) 


VX fo — TG x fy) =P (8) 
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which can be combined to: 
why = VX (VX ho —In Xho) — Pn X (VX ho —ITXX Io) (9) 
Writing Equation (9) in component form gives: 


2 2 2 2 2 : 
we, hy , dhe The The | yy, oho 4, Se 
c Oyox dzdx dy 02" 
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ce) «COzaxCzdy =? By 
os Ohe , Ohy dh , oe) 
o (E45 fy) +r (n. ah ” ay (12) 
—T tie +n 2) 4 al + nyhy + nzhz) 
* Oe ” 02 ae an " 
— The 


where the subscript 0 has been dropped from all components for con- 
venience. 

If the wave in question is an infinite plane wave being propagated 
along the z axis, then: 


Nz = Ny = 0, Nz = 1 


and the components of ho are constants, and h, is zero. For this particular 
case, Equations (10), (11) and (12) become: 


an |i be | (13) 


Say by =, = Thy (14) 
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Equations (13) and (14) are general differential equations derivable 
from Maxwell’s equations and do not yet contain the properties of any 
particular medium. In order to find the behavior of a wave travelling 
through an infinite ferromagnetic medium which is magnetized along the 
direction of propagation, it is necessary to combine these equations with 
Equations (4) which describe the relation between b and hin the medium. 

This gives: 


(uhz — jKhy,) = —JYh, (15) 





we 

Ce 
2: 

(uhy + jKhe) oS =—Th, (16) 


The only possible solution to this set of equations is a circularly 
polarized wave where: 


fee ah, 


The positive sign above represents a so-called positive circularly 
polarized wave and the negative sign a negative circularly polarized 
wave. The propagation constants for these waves is given by 


5 
Ps == Vel + K) (17) 
Pea | 
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Dialing Habits of Telephone Customers 


BY CHARLES CLOS AND ROGER I. WILKINSON 
(Manuscript received October 3, 1951) 


This paper considers the behavior of customers waiting to dial calls, when 
dial tone is delayed. Tests were made in a panel dial central office, from 
which were determined: relationship between load carried by a group of 
line finders and the resultant dial tone delay; measures, by classes of ser- 
vice, of the magnitude of the generalized trunking formula’s ‘4” factor 
describing the degree to which customers wait when dial tone is delayed; 
comparisons of observed and theoretical distributions of the number of 
simultaneous calls on line finder groups; and statistical accounts of the 
actions of customers when dial tone is delayed. 


Following World War II the conversion of great quantities of manual 
telephone equipment to dial, and the addition of large numbers of new 
telephones, mostly dial, in the Bell System has directed increasing atten- 
tion to those service problems peculiar to automatic operation. These 
problems concern chiefly the provision of adequate amounts of equip- 
ment to give satisfactory service at all times. One of the important 
factors affecting the amount of equipment needed is the action of the 
customers themselves when their calls are momentarily blocked due to 
these equipment shortages. The actions of subscribers whose calls are 
blocked due to a shortage of trunk equipment have been reported pre- 
viously.! This paper considers the behavior of subscribers, waiting to 
dial calls, when dial tone is delayed. 

During 1949, Bell Telephone Laboratories conducted a series of tests 
at the New York Telephone Company’s Sterling-3 panel dial central 
office in Brooklyn, N. Y., with the object of increasing the knowledge 
available regarding subscribers’ actions and their effects when dial tone 
is delayed. 

The following principal results were obtained from the Sterling-3 tests: 

1. The relationship between the load carried by a group of line finders 
and the resultant dial tone delay. 

2. Measures, by classes of service, of the magnitude of the general- 


1 Charles Clos, ‘‘An Aspect of the Dialing Behavior of Subscribers and Its 
Effect on the Trunk Plant,’’ Bell System Tech. J., 27, July 1948. 
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ized trunking formula’s ‘‘j”’ factor describing the degree to which cus- 
tomers wait when dial tone is not immediate.? 

3. Comparisons of observed and theoretical distributions of -the num- 
bers of simultaneous calls on the line finder groups. 

4, Statistical accounts of the actions of subscribers when dial tone 
is delayed. 


GENERAL PLAN OF THE TESTS 


The general plan of the tests was developed around two specially 
constructed devices: 

1. A 100-pen recorder capable of recording observations continuously 
for two hours and sensitive enough to record individual dial pulses. 

2. A speed of dial tone measuring set comprising a means for manu- 
ally originating test calls, and an electric timer which automatically 
stopped when dial tone was received. 

Tests were conducted on a weekly schedule, one or two line finder 
groups being studied each week. Two 100-pen recorder tapes were run 
daily. The first was run from. 10 A.M, to noon for message rate individ- 
ual, flat rate individual and message rate two-party classes of service to 
include the morning busy periods for these customers. For coin customers 
the first tape was run for two hours during the noon coin busy period. 
The second tape was run in the afternoon from 2 to 4 P.M. for all 
classes of service except coin. For the coin class the second tape was run 
for two hours during the early evening coin busy period. Three message 
rate two-party tapes were run in the early evening busy period, and two 
coin tapes were run in the afternoons when World Series baseball games 
were being played in Brooklyn. A summary of the number of line finder 
groups observed, the number of tapes taken, the number of line finders 
made available and the maximum per cent dial tone delays over three 
seconds are given in Table I. 

Except for the morning runs on the message rate individual line 
groups where an effort was made to maintain a fixed number of twenty 
line finders throughout, the number of line finders made available was 
selected by close observation of the flow of traffic. All studies were by: 
half hours during which the number of line finders was held constant 
as far as possible. At the end of a half-hour period the number of line 

2 This formula for both finite and infinite sources was developed by R. I. Wil- 
kinson in 1980, and appeared in the 1936 Bell Telephone Laboratories Out-of-Hour 
Course “The Theory of Probability as Applied to Telephone Trunking Prob- 


lems.’’ This formula for infinite sources was also developed by Conny Palm and 
appeared in ‘‘Etude des delais d’attente’”’ in Erickson Technics—No. 2—1987. 
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finders left in service was adjusted in order to obtain a reasonable num- 
ber of dial tone delays in the next half hour without producing a severe 
reaction from the subscribers served by the line groups under study. 

The data recorded on the tapes showed continuously the busy or idle 
conditions of certain circuits associated with the line groups under 
study. In some cases the receipt of dial pulses and the operation of 
registers were also recorded. These circuits included line finders, a few 
subscriber lines, trip circuit sub-groups, the all trunks busy register, 
the peg count register and the speed of dial tone measuring device. In 











TaBLE I 
Numb Maximum 
Line F fader Number ne Pordets ae aoe n 
Groups of Tapes ade Delays. Over 
Observed Available 3 Seconds 
Message rate individual............. 11 
Morning tapes.................... 25 19 to 20 53.3 
Afternoon tapes.................. 29 10 to 15 55.6 
Message rate two-party............. 3 
Morning tapes.................... 8: 18 to 19 40.0 
Afternoon tapes............ ..... ; 12 10 to 12 88.6 
Evening tapes............... 0... 3 19 to 21 71.2 
Flat rate individual................ 2 
Morning tapes.................... 9 14 to 33 48.9 
Afternoon tapes. ................. 10 10 to 17 35.6 
COU iit a Mok and tis. wai gopnd tales ards 3 
MOrninig tADE co fied aos Ga tdane sak 18 30 to 39 71.1 
Afternoon tapes. ................. 2 30 to 39 35.6 
PMVening tapes..6 eke ei aa ee 10 25 to 39 68.9 














addition the busy and idle conditions of a sample of senders was ob- 
served in order to note the general load level on the senders. 


RELATIONSHIP BETWEEN LOAD CARRIED AND PER CENT DIAL TONE DELAY 


One of the principal objectives of these tests was to establish as far 
as possible the relationship between the average load carried by a line- 
finder group and the corresponding dial tone service when there is a 
shortage of line finders but not of senders.’ The average load carried 
was obtained by making a switch count every thirty seconds of the 
number of line finders busy as indicated by the 100-pen recorder tape. 
The dial tone tests were made with the speed of dial tone measuring 
device. Forty-five dial tone tests were made each half hour for each 


3 The restriction of avoiding dial-tone delays due to a sender shortage was to 
eliminate a factor external to the line-finder group and to make the results of the 
tests applicable to both common-control and non-common-control type dial 
systems. 
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line group studied. Additional dial tone tests were made on all other 
line groups in the office as a check that the delays experienced on the 
line groups under study were not due to a sender shortage. The sender 
data on the tapes were also used for this purpose. Figs. 1(a) to 1(d), 
and 2(a) to 2(d), inclusive, show for various amounts of load carried, 
the per cent of dial tone tests encountering delays greater than three 
and greater than ten seconds for half-hour study periods for the most 
frequent number of line finders in the tests for each class of service. 

Plotted on each of these figures is a theoretical fitting dial tone tester 
delay curve computed for the indicated dial tone delays and for the 
following 7 factor values in the generalized trunking formula, determined 
in a manner to be explained later: 











Class of Service j factor 
MRI 6.6 
MR 2-party 5.8 
FRI 6.5 
Coin 2.1 


To indicate the effect of varying 7, several curves have been added 
to Fig. 1(a). Selections of curves for 7 = 0, 1 and ~ (which correspond 
to the three commonly used infinite source congestion formulae, Erlang 
C, Poisson and Erlang B when adapted to the tester’s delay problem) 
are shown. It is clear that with the wide differences in delays which they 
give for specified loads carried, it is highly desirable to select that 7 
formula for engineering use which most nearly describes the customer 
actions in any situation being dealt with. In the field of curves shown 
on Fig. 1(a), the one labelled 7 = 6.6 was derived in a logical manner 
from the data, and shows an agreeably satisfactory fit. For example, 
during a heavy load period when, say, 20 per cent of a dial tone tester’s 
calls are meeting delays greater than 3 seconds, an actual average load 
of 16.6 erlangs (as shown by the 7 = 6.6 curve) would likely be carried. 
(Load in erlangs equals average number of simultaneous calls.) The 
Erlang C (7 = 0) and Poisson (j = 1) theories would indicate the 
presence of loads of 15.6 and 16.0 erlangs, respectively, figures clearly 
too small for the circumstances shown by the data of Fig. 1(a). On the 
other hand, use of the Erlang B (j = «) theory would predict a con- 
siderably larger load carried, about 17.9 erlangs, than one would prob- 
ably be justified in assuming here for engineering purposes. 

By grouping the dial tone delay data by bands of load carried, rela- 
tionships of per cent of test calls encountering varying dial tone delays 
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up to twelve seconds were obtained. These data are shown on Figs. 3 
to 6. Fig. 3 is for the message rate individual class of service with twenty 
line finders available. The data on this figure correspond to those on 
Figs. 1(a) and 1(b). Fig. 4 is for the message rate two-party class of 
service with ten line finders and corresponds to Figs. l(c) and 1(d). 
Fig. 5 is for the flat rate class of service with ten line finders and cor- 
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responds to Figs. 2(a) and 2(b). Fig. 6 is for the coin class of service 
with 34 line finders and corresponds to Figs. 2(c) and 2(d). 

Plotted on Figs. 3 to 6 are theoretical fitting dial tone tester delay 
curves, curves A, determined by means of the following formulae: 

1. The generalized trunking formula for determining the proportion 
of calls that encounter congestion, i.e., find all line finders busy. 
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OBSERVED DIAL TONE TESTER DELAYS 


mena THEORETICAL DIAL TONE TESTER FITTING CURVE - 
TEST CALL WAITS ITS TURN IN LINE 


THEORETICAL DIAL TONE TESTER FITTING CURVE — 
TEST CALL IS ALWAYS FIRST IN LINE 


Pr —--—— THEORETICAL DIAL TONE SUBSCRIBER DELAY CURVE 


| MESSAGE RATE INDIVIDUAL 
(20 LINE FINDERS AVAILABLE) 


ite Set 





3 4 5 6 7 8 9 10 14 
TIME, t,IN SECONDS FROM START OF DIAL TONE TEST 
Fig. 3—Results of dial tone tests. 

2. A delay formula! for determining the proportion of those test 
calls encountering congestion which will have a delay in obtaining dial 
tone of at least time ¢. 

Some of the theoretical aspects of these formulae are considered in 
the two sections that follow. 


GENERALIZED TRUNKING FORMULA 


The generalized trunking formula combines in one expression the 
various assumptions underlying the Erlang B, Poisson, and Erlang C 


4A development by John Riordan paralleling that of Conny Palm’s. 
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0 1 2 5 7 8& 9 10 1 {2 13 14 
TIME, t; IN SECONDS PROM START OF DIAL TONE TEST 


Fig. 4—Results of dial tone tests. 


trunking formulae and a large field of intermediate assumptions regard- 
ing the disposition of calls which fail to obtain immediate service. These 
assumptions are given in Table II on page 42. 

One method for developing the generalized trunking formula® is: to 
consider the probabilities of the existence of certain states and to de- 
termine the number of transitions during some convenient interval of 


5 In this article only the unlimited sources trunking formula is considered. 
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time from one state to another. By equating certain of these transitions, 
a series of simultaneous equations evolves, which when solved yields one 
overall expression. Of interest is the development of the transition equa- 
tions. Thus for a case of ¢ line finders arranged in a simple group, let 
f(x) represent the probability that x (where x < c) line finders are 
occupied and f(« + 1) represent the probability that x + 1 line finders 
are occupied. Let n be the average number of calls offered to the line 
finders during a long interval of time, 7’. Let 7 be the unit of time and h 
be the average holding time per call measured in terms of 7. Over a 


OBSERVED DIAL TONE TESTER DELAYS 


——— THEORETICAL DIAL TONE TESTER FITTING CURVE - 
TEST CALL WAITS ITS TURN IN LINE 


THEORETICAL DIAL TONE TESTER FITTING CURVE - 
TEST CALL IS ALWAYS FIRST IN LINE 
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TIME, t. IN SECONDS FROM START OF DIAL TONE TEST 


Fig. 5—Results of dial tone tests. 
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period 7’ the number of transits from state x to state x + 1 must equal 
(or differ by no more than one) the number of transits in the reverse 
direction. That is: 


nfo) = =F * fe +1) (1) 


It may be noted that nh = a, where a is the average offered traffic 
load in erlangs. 
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Fig. 6—Results of dial tone tests. 
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TABLE II 





Assumption Concerning the Disposal of Calls that do not Obtain 


Formula a Line Finder Immediately 








Generalized........ Waiting calls are cleared out at a rate j times the rate at 
which calls are terminated when served by the line 
finders. 

Erlang B.......... Calls are cleared out of the system immediately, that is 
no calls wait (j = ©). 

Poisson............ Waiting calls are cleared out at a rate equal to that with 
which calls are terminated when served by the line 
finders (j = 1). 

Erlang Cjccc.si es: Calls wait until served (j = 0). 








When x 2 ¢ a new situation is encountered. ‘‘c” calls are engaged in 
conversation and x — ¢ calls are waiting for service. If the waiting calls 
are forced to wait for an unduly long period of time so that in effect 
they are being denied service, it can be expected that they will wait for 
some average period, say H, and then abandon their attempts. On this 
basis the corresponding equation is: 


nfla) = fet 1) +247 —9 ge 41) @) 


It has been assumed in the above equations that the distribution of 
the holding times is exponential, an assumption which is found in most 
local systems to be reasonably justified. The distribution of the waiting 
times is also taken to be exponential. By introducing a factor 7, where 
j = h/H, equation (2) can be written in the simpler form: 


af) = [e+ je +1—oOlfe +1) (3) 


Solving this system of simultaneous equations, we obtain: 
when x < ¢, 


when «2 ¢, 
a 
Ite) 7 cl(c + 9) (c + 2]) wd eis hee sie [c +- (x = ol f(0) (5) 


where 


z=c 2£ r7=0 x 


a ‘* 
10 = (25+ 3 aepnerm hte) ” 
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The probability of a call encountering congestion, which is equivalent 
to the probability of a call having a delay greater than zero units of 
time is: 


ZO 


P(>0) = 2 fe) (7) 


DELAY FORMULA FOR THE DIAL TONE TESTER 


The probability that a dial tone test call encounters congestion is 
given by expression (7). Once a test call has encountered congestion it 
will experience a delay depending upon a number of variables. The 
assumptions underlying the dial tone tester formula are: . 

1. A dial tone test call when encountering a delay waits until served. 

2. A dial tone test call does not add to the load offered and carried 
by the line finders. 

3. Upon encountering a delay, a dial tone test call is served in the 
order of its arrival with respect to all other waiting calls. For example, 
if the test call finds three other calls waiting, it waits fourth in line. 

Under the third assumption as calls drop out, due to conversations 
terminating on the occupied line finders or due to waiting calls aban- 
doning their requests for service, the test call advances from an initial 
position of say fourth in line to third in line, then to second, then to 
first in line, and finally is served. The overall delay distribution of the 
test calls depends therefore upon the number of calls they find waiting 
ahead of them. The delay distribution for each such number must be 
weighted by the probability of its occurrence in order to obtain the 
overall distribution. The delay distribution for a test call which finds 
zero calls waiting is: 


po(>t) = exp (—ct/jH) (8) 


The probability is f(c) that a call made at random will find all line 
finders busy with no calls waiting. Hence the weighted delay distribution 
P o( >t) ) iS: 


Po(>t) = f(e)po(>t) = f(c) exp (—ct/jH) (9) 


The delay distribution for a test call which finds one call waiting 
ahead of it is: 


pi(>t) = [1 + ¢/j — (6/9) exp (—t/H)] exp (—ct/jH) (10) 
The probability is f(¢c + 1) that a call made at random will find all 
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line finders busy and one call waiting. Hence the weighted delay distri- 
bution, Pi(>1), is: 
Pi(>t) = f(e+ 1[1 + c/j — (c¢/9) exp (—t/H)] exp (—et/jH) (11) 
In the general case p,n(>#) is given by the following formula: 
Pr( >t) = Fraii(t) exp (t(/H) (12) 


where F’,,4:(t) is given by Conny Palm.® The over-all delay distribution 
is then: | 


P(>t) = Po(>t) + Pi(>t) + P(>t) + --- (13) 


By making appropriate substitutions and summing the result, ex- 
pression (13) becomes: 


4 aexp (—t/H) , a exp (—2t/H) 
Pood =F 401 + rig. Caer ap | 


exp | —ot/jH + (a/j)[L — exp (-/H)] (14) 


Expression (14) is equivalent to that of Riordan involving two in- 
complete gamma functions as follows: 


P(>t) = P(>0) yle/i, (a/3) exp (~t/H)| (15) 
y(c/j, 4/9) 
where the incomplete gamma function, 


y(N, x) = [ ate dx (16) 
0 


The theoretical dial tone tester delay curves shown on Figs. l(a) to — 
1(d), 2(a) to 2(d), and 3 to 6 were computed from expression (14), 
using the following values of 7 and H for the classes of service studied, 
these values being determined in a manner explained later: 








Class of Service j factor | H 

MRI 6.6 24 seconds 
MR 2-party 5.8 42 seconds 
FRI : 6.5 27 seconds 
Coin 2.1 74 seconds 








On Figs. 1(a), 1(c), 2(a), and 2(c), which show the per cent of dial tone 
tests encountering delays greater than three seconds for various amounts 


§ Hquation 53, loc. cit. 
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of load carried, it may be noted that most of the theoretical dial tone 
tester delay curves are in close agreement with the observed data, with 
a tendency perhaps to be slightly high. On Figs. 1(b), 1(d), 2(b), and 
2(d), which are for dial tone delays greater than ten seconds, it may be 
noted that the theoretical curves have a slightly stronger tendency to 
lie on the high side of the observed data. On Figs. 3 to 6 the theoretical 
dial tone tester delay, curves A, again lie in the proximity of the curves 
of the observed data, with a tendency to lie higher than these latter 
curves, especially at the ends where the dial tone delays are greatest. 
Among the factors which account for this discrepancy are: 

1. A feature is present in panel line finder circuits for momentarily 
releasing trip circuits with waiting calls to prevent the orphaning of 
calls under certain trouble conditions. The release occurs after a call 
has been waiting from 5 to 12 seconds and reoccurs every 7 seconds 
thereafter. When such a release occurs the call yields whatever waiting 
preference it may have had to a subsequently placed call which is not 
yet affected by such a release. The dial tone test calls did not wait 
beyond 12 seconds. Hence for these test calls there was only one possi- 
bility of such a release and for many of them the release occurred near 
the end of their waiting period. Hence they were more likely to gain 
preference over cther calls than to lose their preference. 

2. Subscribers while waiting for dial tone frequently become impatient 
and proceed to flash (move their switchhook up and down). While 
flashing, a subscriber may lose preference to a subsequently placed test 
call (the latter of course does not flash). 

3. Many subscribers fail to observe dial tone and proceed to dial. 
During such dialing, a subscriber may lose preference to a subsequently 
placed test call. 

4. Line finders serve a large proportion of call attempts of short 
holding time whose presence may militate against the occurrence of 
the longer delays. In connection with the measurement of the j factor, 
the following proportions of call attempts and average holding times 
were noted on which no dialing occurred or where no more than two 
digits were dialed. 

















Class of Service SE ene Coe with No Average Holding Time 
MRI 35.2%* 4.4 seconds* 
MR 2-Pty. 33.3%* 5.8 seconds* 
FRI 25.1% 5.8 seconds 
Coin 9.8% 8.3 seconds 








* Partly estimated 
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The individual contributions of these four factors to discrepancies 
between theory and observation are not easy to assess. The first three 
explain a tendency for test calls to get ahead of calls already waiting 
for dial tone. On Figs. 3 to 6, inclusive, additional theoretical dial 
tone delay curves, curves B, for the case where a dial tone tester 
always gets first in line are shown. Even these curves tend to lie above 
the curves of the observed data on Figs. 4 and 5 where ten line finders 
were available; they more nearly agree with the observed data on Fig. 3 
where twenty line finders were available, and they lie below the observed 
data on Fig. 6 where 34 line finders were available. This is an indication 
of the fact that with higher traffic loads (which occurred on the larger 
line finder groups) a test call will encounter more competition from 
other calls and therefore will have a lesser chance of gaining precedence 
over all of the other calls. The fourth factor indicates that the call 
attempts served on line finders consist of two distinct holding time 
universes and not just one, as was assumed in the development of the 
dial tone tester formula. The effect of the presence of both a short and 
long holding time universe of calls would be to introduce a change of 
slope in the delay curves which may be seen in Figs. 3 to 6 to be at about 
t = 4 seconds. There is reason to believe that the same cause may have 
been responsible for the tendency of the observed delay curves to fall 
away from the theoretical at the lower levels of load carried. 

Due to the reasons given above and to the fact that the dial tone 
delay observations were made by the test call method, the above re- 
sults may not directly describe service from the customer’s point of 
view. Conny Palm has developed the following formula which gives a 
slightly different measure of customers’ dial tone service. It indicates 
the proportion of calls which have neither received dial tone nor have 
dropped out at time ¢. 





P(>t) = P(>o) Ue Gi) exp (A) ong (_y) (17) 
y(c/j, 4/9) 7 

Curves for this formula are shown plotted on Fig. 1(a) and at C on 
Figs. 3 to 6. They are quite close in many cases to the observed dial 
tone tester results. It would appear that a sufficiently good estimate 
of the customer’s dial tone service, whatever its precise definition, can 
be obtained by the dial tone tester method. . 

Recently revised tables for the capacity of step-by-step line finders 
have been published for Bell System use based on Palm’s formula using 
a factor of 7 = 5. This was selected as being slightly conservative for 
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most applications after reviewing the above Sterling-3 results and other 
line finder data collected in step-by-step offices. 


MEASUREMENT OF THE J FACTOR BY CLASSES OF SERVICE 


As indicated previously, the data recorded on the tapes showed the 
states of being busy or idle and of changes in these states for line finders 
and the associated trip circuits. A fully equipped line finder group of 
400 lines has ten trip circuits each of which serves two sub-groups of 
twenty subscriber lines in the following manner. When a line originates 
a call its line relay is operated. This causes a ground to appear on a 
lead which is common to all twenty line relays in the sub-group and 
starts a line finder hunting for the calling subscriber’s line. As soon as 
this hunt is completed the cutoff relay associated with the calling line 
operates and disconnects the line relay, removing the ground (unless, 
of course, another line in the sub-group has originated a call in the 
meantime). During periods of overload when line finders are not im- 
mediately available, the ground due to a single subscriber will persist 
until: 

1. A line finder is obtained, or 

2. The subscriber abandons the attempt, or 

3. The subscriber receives an incoming call which operates the cut- 
off relay. 

The twenty leads from the trip circuit sub-groups were brought out 
to the pen recorder and a record taken of the grounds that occurred on 
each lead. Except for the possibility that more than one subscriber is 
waiting for service at the same time on a given trip circuit sub-group, 
the record of the occurrences of the grounds gives a substantially ac- 
curate’ record of the demands for service and of the number of calls 
waiting for service. Hence an analysis of the events occurring on the 
trip circuit sub-groups and on the line finders as recorded on the tapes 
gives a means for determining H. The quantity H was introduced in 
equation (2) in the term 


(x +1-—e) 


i fe + 1) (18) 


For convenience in the ensuing discussion this term will be replaced by 
7 To obtain absolute accuracy would require the use of a pen recorder with one 


pen for each of the 400 subscribers served on a line finder group plus one for each 
line finder. 
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the equivalent expression: 
Ny = aie Yj) (19) 


where N, = The average number of waiting calls that drop out per unit 
of time during the state (z + y). 
y = The number of waiting calls. 
z = The number of line finders occupied with calls. . 
1/H = A measure of the rate at which calls tend to abandon 
waiting. 
f(e + y) = The proportion of time that the state (2 + y) exists. 


On the tapes we can measure f(z + y) and count N,. Hence H can 
be determined. The result is a statistical quantity subject to many 
chance factors. In the actual analysis of a tape, the composite average 
value of H was determined for all possible observed states where calls 
were waiting. By an analogous process, the composite average value 
of h for all possible observed states where calls were being served by 
line finders was determined. Also as a side computation, a composite 
average value of h’ for calls that were served by line finders but for 
which no peg counts were scored was determined. This value of h’ is 
included in h on the basis that data for engineering line finders consist 
of estimated calls based on peg counts and of holding times which in- 
clude an allowance for these short holding time calls. The average 
values of H, h and of the 7 factor for the four classes of service studied 
are given in Table III. 

The results for H by individual half hours and by various percentages 
of dial tone delays greater than three seconds are shown on Figs. 7(a) 
to 7(d) respectively for the four classes of service. On some of these 
figures an upward bulge may be noted in the center. This is not con- 
sidered to be characteristic of the habits of the subscribers but is the ~ 
overall effect resulting from a number of arbitrary rules followed in 
making the analysis in order to simplify the work and to offset par- 














TasueE III 
Average Values in Seconds 
‘g=h/H 
H h 
Message rate individual..................... 24 159 6.6 
Message rate two-party. ..............000005 42 243 5.8 
Flat rate individual......................... 27 176 6.5 
COMM cise. abe teseite athgia idee Guede ks 74 153 2.1 
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Fig. 7—Average customer waiting time (ff). 


tially the effect of occasionally having two or more calls waiting on one 
trip circuit sub-group. The rules and the reasons for them will be de- 
scribed with the aid of Fig. 8, which shows a hypothetical section of 
one of the tapes. The rules were as follows: 


1. Initial Overlap 


Referring to Fig. 8, at t, a subscriber has initiated a request for 
service. At ¢: a line finder rises to serve the subscriber. At ¢3 the sub- 
scriber receives service. This case is typical of a subscriber receiving 
prompt dial tone service. 

The span from ¢; to fg was difficult to measure accurately because, 
for the usual case, it was about the same as the maximum error due to 
misalignment of the recorder pens. It was not measured unless the com- 
bined span from t; to t3 exceeded one second. 

The span from t, to ¢; involves an overlap, it represents a period 
when a line finder is busy hunting for the terminal of the subscriber who 
originated the request for service. It also represents a period when a 
subscriber is waiting for service. In: the analysis this span was treated 
as a case where a line finder was busy with a call and not as a call wait- 
ing for service. 

If the span from t, to tz; and all similar cases had been treated as 
calls waiting for service and if in addition all spans from ¢; to t2 which 
were not measured had also been treated as calls waiting for service, 
the average values for H would have increased slightly for each class 
of service. : 
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2. Three Second Rule for the Bridging of Calls 


Referring to Fig. 8, again, at t; a request for service is originated on 
trip circuit 5 and at t; this request is withdrawn. At t» apparently a 
new request for service is initiated which is then withdrawn at tu. From 
manual service observations it is found that subscribers often flash when 
dial tone is slow. A few pens were used to observe individual subscribers, 
and Fig. 9 shows a case where a subscriber made several flashes when 
his tone was slow. When a subscriber flashes it appears as though he 
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Fig. 8—Section of a hypothetical tape showing activities on trip circuit sub- 
groups and on line finders. 
were making several bids for service. Actually he is making only one 
real bid. On the trip circuit sub-group pens it was generally impossible 
to distinguish between flashes and requests for service by two or more 
subscribers. To resolve this problem many observations of subscriber 
lines recorded on the tapes were examined from which it was concluded 
that no great error would result if a break in the demand for service 
on a trip circuit sub-group of less than three seconds were considered 
as a flash and was to be bridged, and a break greater than three seconds 


was to be considered as the termination of one call attempt and the 
start of another 
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8. Treatment of Cases Where Two or More Calls Were Found to be 
Waiting on One Trip Circuit Sub-Group 


The occurrence of several calls waiting on one trip circuit was occasion- 
ally noted in the analysis. Referring to Fig. 8, a case is shown on trip 
circuit sub-group 12. At ty a line finder is seized. Trip circuit sub-group 
18 shows that a subscriber is dialing before tone. The appearance of 
dial pulses on this trip circuit indicates that only one subscriber is 
demanding service otherwise the dialing would not show. Trip circuit 
sub-group 12 however appears to have two or more requests for service. 
One of these requests for service began at ty. The start of the second 
request occurred somewhere between t, and t, perhaps half-way be- 
tween. At t, one of the requests was served by a line finder. To sim- 
plify the handling of such cases, the assumption was made that the 
first attempt started at t,; and ended at ty) and the second request started 
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Fig. 9—Example of a customer flashing for dial tone (Tape made October 18). 


at t). The effect of this is to understate by an indeterminate amount 
the average value of H for each class of service. This understatement 
should be noticeable for the higher degrees of overload because the 
occurrence of several calls on one trip circuit sub-group is then most 
likely to occur. 

The effect of several calls simultaneously waiting in a trip circuit 
sub-group and of one or more calls dropping out is to overstate the 
magnitude of H. For instance if two simultaneous call attempts of five 
seconds each overlap for one second and both attempts are abandoned, 
the apparent average waiting period is nine seconds, whereas it should 
be five seconds. It is believed that the above three rules tend to create 
understatements which roughly balance this type of overstatement. 


DISTRIBUTION CURVES OF SIMULTANEOUS CALLS 


The detailed analysis of the tapes provided distributions of simul- 
taneous calls. For each class of service studied these distributions can 
be compared with theoretical distributions derived from the generalized 
trunking formula using the j factors developed in the analysis. Several 
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such comparisons are shown on Figs. 10 and 11. The agreement is quite 
good in most cases. 


SUBSCRIBER DIALING HABITS AS OBSERVED WITH A MONITORING CIRCUIT 
ON A SENDER WITH INDUCED DIAL TONE DELAYS® 


' As a separate study a series of tests was made by means of a moni- 
toring circuit on one of the senders serving in common the subscribers 
in the Sterling-8 and Main 2 central offices, for the purpose of obtaining 
further information on subscriber dialing characteristics under overload 
conditions. A large amount of data was collected on the time intervals 
from the seizure of the sender to the first action taken by subscribers 
when encountering dial tone delays, the latter being introduced under 
the control of the observer. 

The monitoring circuit was wired to a particular sender in a group 
of 100 serving all classes of subscribers. When the circuit was in use, the 
only irregularity introduced was that the dial tone could be delayed 
even though the sender was actually available to the subscriber. The 
delay did not affect the sender in its functions if the subscriber elected 
to dial before tone. 

The sender monitoring circuit provided the following four features: 

1. A receiver was bridged across the tip and ring leads in the sender 
so that an observer could hear certain actions taken by a subscriber 
connected to the sender. The sender was of course disconnected before 
conversation. 

2. The observer was able to preselect one of several intervals by 
which dial tone was delayed on successive calls served by the sender. 
This was accomplished with a capacitance-resistance-vacuum tube cir- 
cuit. 

3. By means of a timer which started when the sender was seized, 
_the observer was enabled to note elapsed time intervals to the occurrence 
of the various actions of the subscribers. The reading of the time of 
the first action of a subscriber had to be made when the second hand 
was in motion, which introduced certain errors later to be discussed. 

4, By means of colored lamps the observer was able to classify all 
calls observed as being message rate, flat rate or coin. 

During the sender dial tone delay tests, observations were made only 
during the afternoons when the flow of traffic was light and the prob- 
ability of a subscriber obtaining a delay before reaching the sender was 
a minimum. 


8 Based on an unpublished report by W. A. Reenstra. 
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DIALING HABITS OF TELEPHONE CUSTOMERS DO. 


The observer was provided: with a means for introducing either no 
delay or one of four values of delay 2, 5, 10 or 15 seconds into the sender 
dial tone circuit. The observer took 50 observations using a particular 
value of dial tone delay and then shifted to another so that no particu- 
lar value of delay would become evident to the customers during an 
afternoon’s test. Each group of 50 observations comprised a mixture of 
message rate, coin and flat rate calls in the approximate proportions of 
13 to 6 to 1, representing the respective volumes of traffic from these 
classes of service during the afternoon periods. It was not possible to 
distinguish PBX lines or two-party lines from the bulk of the message 
rate data nor PBX lines in the flat rate data, although to a limited 
extent the observer could identify PBX dialing by the generally faster 
pulsing. The coin data represent both public and semi-public customers. 

Fig. 12 is a diagram for explaining the results shown on Figs. 13, 14 
and 15 for the message rate, flat rate and coin classes of service, respec- 
tively as obtained with the sender monitoring circuit. Fig. 12 was ob- 
tained by the application of fitting curves to those message rate data of 
Fig. 13 for which a dial tone delay of five seconds was introduced by the 
observer. In the interval from ¢ = 0 to £ = 5 seconds, three curves A, B 
and C represent the per cent of subscribers still waiting at time ¢ for dial 
tone. Curve A and its extension beyond ¢ = 5 seconds represents the 
~ action of subscribers who would dial their calls before tone if dial tone 
were delayed indefinitely. Curve B and its extension beyond t = 5 
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Fig. 12—Explanatory chart for sender monitoring observations; dial tone at 
T = 5 seconds. 
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seconds represents actions of subscribers who would disconnect if dial 
tone were delayed indefinitely. Curve C and its extension is the sum of 
the other two. Curve D in the region beyond ¢ = 5 seconds represents the 
actions of subscribers who disconnect after tone, curve E represents the 
actions of subscribers who will dial their calls after tone and curve F 
represents the sum of the lower curves. Of interest is the fact that for an 
interval of about two seconds following dial tone (at ¢ = 5 seconds), 
the observed total survivor curve F lies above the extended portion of - 
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the hypothetical survivor curve C for infinitely delayed dial tone. This 
indicates that most of the subscribers who would have abandoned their 
attempts during this interval abruptly changed their minds and then 
consumed a noticeable interval of time after hearing tone before start- 
ing to dial. Thus, as might be anticipated, the subscribers exhibit a 
reaction time. . 

Tig. 13 shows the results in terms of survivor curves that were ob- 
served for the message rate class of service. Five sets of curves are 
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presented, namely for 0, 2, 5, 10 and 15 seconds of dial tone delays from 
the instant of sender seizure. It should be noted that the general con- 
tour of the various curves up to the receipt of dial tone and when ex- 
tended beyond gives an estimate of the survivor curve for dial tone 
delayed indefinitely. Fig. 14 shows the results for the flat rate class of 
service (the data here are relatively meager), and Fig. 15 shows the 
results for coin customers. 

Fig. 16 indicates for the three classes of service the progressive changes, 
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Fig. 16—Results of sender monitoring observations. 


with increasing dial tone delay, of the per cent sender seizures resulting 
in dialing before tone, in dialing after tone and in disconnections for the 
five dial tone delay intervals studied. 

Some general comments concerning Figs. 13 to 16 may be made. 

1. There is a striking contrast between the message rate and coin 
classes of service. This may be due to the immediate financial stake 
that a coin customer has in his call. He is reluctant to disconnect before 
dial tone. 

2. The results for the message rate and flat rate classes of service 
appear to be similar at the shorter dial tone delays; at the longer delays 
a higher proportion of flat rate subscribers have already dialed before 
dial tone. This apparent discrepancy may be due to the relatively small 
number of flat rate observations. Flat rate service in Sterling-3 and Main- 
2 was principally for professional people, such as doctors and nurses, 
who were thought to be more demanding than ordinary subscribers for 
prompt dial tone service. These results therefore should not be con- 
sidered characteristic of flat rate customers generally. 

3. Detailed analysis of the data (results not presented in this article) 
indicated that the distributions of time to first pulse for subscribers not 
observing tone, and for those waiting for tone, are quite similar for the 
three classes of service. 

4. Only the unsmoothed raw data have been shown on Figs. 13 to 
15 since certain inadequacies were detected in the observations. These 
were due to observer reactions and reading errors discovered as a result 
of comparing preliminary practice sender monitor test results with a 
simultaneous record obtained by the 100 pen tape recorder. The overall 
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effect is that the data obtained by the observer are generally displaced 
outward along the time axis by about 0.8 second. 

5. The message rate data were for individual, PBX and two-party 

subscribers and the flat rate data were for individual and PBX sub- 
scribers. Furthermore, certain of the flat rate subscribers had auxiliary 
message rate service. It seems likely that different characteristics would 
be obtained for the individual, PBX and two-party subscribers since 
there appear to be reasons for expecting significant differences in their 
dialing habits. The PBX operator is in a position to ‘‘shop”’ for tele- 
phone service. If she fails to get dial tone on one outgoing line, she can 
try any other free line. This can also be done by subscribers with multi- 
line service. This ‘“‘shopping”’ for service tends to produce a large volume 
of disconnections when dial tone is slow. The individual and the two- 
party subscribers cannot do this and hence they can be expected to 
show fewer disconnections. 

6. The results are in terms of intervals of time from the instant the 
sender is seized. It would, of course, be preferable to have these results 
in terms of time from receiver off hook. Since on the average the sender 
is seized in a time interval of about the same magnitude as that of the 
reaction time of the observer, Figs. 18 to 15 can be read approximately 
correctly when the abscissas are redesignated “time in seconds from 
receiver off hook.’’ The foregoing results have been presented to furnish 
an increased understanding of subscriber dialing habits. In the next 
section additional results based on individual line records taken on the 
tapes are presented. 


SUBSCRIBER DIALING HABITS OBSERVED BY INDIVIDUAL LINE RECORDS 


As indicated in the previous section, the results obtained by means 
of the sender monitor tests were subject to certain shortcomings, hence 
data taken on individual lines with the 100-pen recorder have been 
analyzed to augment the information concerning the dialing habits of 
subscribers. 

As noted heretofore, several of the pens on the 100-pen recorder were 
available for taking observations on subscribers lines. Two pens were 
used per subscriber line, one recorded the operation of the subscriber’s 
line relay while the other pen marked whenever the subscriber’s line 
was busy. On an originating call both pens started marking when the 
subscriber initiated a call. When a line finder was obtained, the line 
relay pen ceased marking and it was presumed that the subscriber 
obtained dial tone at that instant. Dialing, hang-up and flashing by a 
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subscriber after receipt of dial tone are noted by breaks in the markings 
of the line busy pen. Dialing, hang-up or flashing before dial tone are 
noted by simultaneous breaks in the markings of both pens. Various 
intervals can be measured and the call attempts classified accordingly. 

Observations were obtained in the foregoing manner during the course 
of the Sterling-3 line finder tests on the following numbers of subscriber 
lines: 
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The observed data were classified for each of the above six types of 
subscribers in terms of the following categories: 

1. Time to subscriber action before receipt of dial tone. 

a. Time from receiver off hook® to first digit dialed by subscriber. 
b. Time from receiver off hook to disconnecting action by the 
subscriber. 

2. Time from receiver off hook to receipt of dial tone. 

3. Time to subscriber action after receipt of dial tone. 

a. Time from receipt of dial tone to first digit dialed by subscriber. 
b. Time from receipt of dial tone to disconnecting action by the 
subscriber. 

Because the data developed in this section are compared with both 
the j factor analysis and the sender monitor test results and because 
the treatment of subscriber dialing habits before dial tone for each of 
these items is different, categories la and 1b are analyzed in two ways. 
In the 7 factor analysis, all actions of a subscriber prior to dial tone, 
except a disconnect, were considered to be one continuing demand for 
service. Thus for the first analysis, cases of dialing, flashing and short 
disconnections before tone lasting less than three seconds were ignored. 
In the sender monitor tests the only items considered were the time to 
the first digit dialed by a subscriber and the time, if no dialing occurred- | 
to the release of the sender by the subscriber. Thus for the second analy- 
sis, cases of dialing before tone and flashing or disconnections before tone 

? When a customer initiates a call, the line relay operates. For individual line 
and two-party subscribers this occurs when the subscriber takes his receiver off 
the hook. For coin customers this occurs when the customer has taken his re- 
ceiver off the hook and made a proper deposit. For a PBX line this occurs when 


the PBX attendant -has established a connection to an outside line. All of this 
is collectively termed “receiver off hook.’’ 
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Fig. 17—Subscriber dialing habits based on individual line records when dia 
tone is delayed indefinitely. 


that caused the sender to release were each counted as separate attempts. 
The results arrived at are shown as survivor curves on Figs. 17(a) to 
19(b) for each of the six types of subscribers studied. 

The survivor curves were developed by considering events during 
0.1 second intervals. The number of cases of dialing before tone and 
the number of cases of disconnection before tone occurring during a 
0.1 second interval were divided by the number of cases waiting for dial 
tone at the start of the interval. This ratio was considered to be a 
retirement rate. The complement of this rate gave a survival rate. By 
a progressive multiplication of survival rates from time, t = 0, the 
resulting survivor curves were obtained. Those cases receiving dial tone 
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during a particular 0.1 second interval are omitted from the number of 
cases waiting for dial tone at the start of the next interval. 

In the development of the generalized trunking formula, the assump- 
tion was made that the waiting times of calls infinitely delayed have an 
exponential distribution. By assuming that the plots for the survivor 
curves in the development comparable to the j factor analysis are ex- 
ponential distributions, it is possible by reading the value of t correspond- 
ing to 36.8 per cent of the subscribers still waiting for dial tone to obtain 
estimates of the values of H for the six types of subscribers. These esti- 
mates, most of which were obtained by extrapolation, are compared in 
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Fig. 18—Subscriber dialing habits based on individual line records when dial 
tone is delayed indefinitely. 
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Fig. 19—Subscriber dialing habits based on individual line records when dial 
tone is delayed indefinitely. 

Table IV with the values developed earlier in connection with the 
j factor analysis. 

These results are not considered to be jnccnaiatont since the tails of 
the survivor curves were constructed from very meager data. Conclu- 
sions based on Figs. 17(a) to 19(b) should therefore be regarded as 
having qualitative value only. The results principally indicate that the 
waiting-for-dial-tone characteristics of subscribers clearly vary with the 
different classes of service. 

Comparisons between the lower survivor curves on Figs. 17(a) to 
19(b) and the curves on Figs. 13 to 15 of the sender monitor tests are 
indicated by the percentages given in Table V of subscribers waiting for 
dial tone 5, 10 and 15 seconds from the time they requested service. 
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TABLE IV 





Estimated Values of H in Seconds 
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Percentages of Subscribers Waiting for Dial Tone 
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Fig. 20—Distributions of the start-to-dial times of subscribers who dial after 
the receipt of dial tone. 
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Fig. 21—Distributions of the start-to-dial times of subscribers who dial after 
the receipt of dial tone. 


These results agree reasonably well when it is recalled that parts of 
the individual line data were scanty and that the sender monitor tests 
included the effects of observer reactions. 

Once dial tone is received, it appears that all types of subscribers 
tend to follow a uniform dialing pattern. Figs. 20 to 23 show for a class 
interval of 0.5 second the distributions of the per cent of subscribers who 
dial at time ¢ for the six types of subscribers studied. Figs. 20, 21 and 22 
show the distributions when dial tone is received from 0.0 to 0.5, 0.6 to 
0.9 and 1.0 to 1.9 seconds after dial. tone, respectively. These curves 
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Fig. 22—Distributions of the start-to-dial times of subscribers who dial after 
the receipt of dial tone. 
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Fig. 23—-Distributions of the start-to-dial times of subscribers who dial after 
the receipt of dial tone. 


indicate that a strong similarity exists among the six types of subscribers 
with regard to their dialing patterns once dial tone is received. 

Fig. 23 shows composite dialing distributions of the six types of sub- 
scribers for four dial tone delay intervals. For dial tone delays less than 
3 seconds the dialing patterns are seen to be similar. Beyond 3 seconds 
delay the start-to-dial curve shifts outward, although the data are too 
scattered to indicate closely where the movement begins. 


CONCLUSION 


The foregoing report of the results of the tests conducted at the 
Sterling-3 central office indicates that the dialing habits of subscribers 
waiting for dial tone can be observed and analyzed to develop so-called 
j factors for use in a more general trunking formula than has been 
employed until recently in the Bell System. The report also presents 
descriptive data regarding the patterns subscribers follow when waiting 
for dial tone. 
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Selective Fading of Microwaves 


BY A. B. CRAWFORD AND W. C. JAKES, JR. 


(Manuscript received October 25, 1951) 


The results of an extended survey of microwave propagation over two line- 
of-sight paths in New Jersey are described. Angle-of-arrival measurements 
at 1.25-cm wavelength and selective fading observations in a 450-mce fre- 
quency band centered at 8950-mce show that the severe fading can be explained 
in terms of multtple-path transmission. A computer of the analogue type 
was built to simulate the more complicated selective fading patterns. 


INTRODUCTION 


During the past few years, studies of microwave propagation have 
been made by the Radio Research group at the Holmdel Laboratory 
over two paths located in eastern New Jersey. Both of these are line- 
of-sight paths which might be considered to be typical links in a cross- 
country microwave radio relay circuit. 

In conducting these studies, the usual continuous recordings of signal 
levels were made but the greater interest was centered in special experi- 
ments designed to reveal more of the processes which can cause fading. 
The most relevant information has been obtained by exploring the inci- 
dent wavefronts with a narrow-beam scanning antenna (angle-of-arrival 
studies) and, more recently, by observing the transmission characteris- 
tics of the paths by means of a frequency-sweep technique and also by 
the use of very short pulses. 

Some results of angle-of-arrival observations have been reported pre- 
viously! and a companion paper describes the transmission tests con- 
ducted with very short pulses.? The present paper describes some of the 
observed mechanisms associated with fading, presents typical data ob- 
tained with the narrow-beam scanning antenna and gives examples of 
the frequency-sweep observations, illustrating the frequency selective 

1W. M. Sharpless, ‘Measurement of the Angle of Arrival of Microwaves,”’ 
Proc. I.R.E., 34, Nov. 1946, pp. 837-845. A. B. Crawford and W. M. Sharpless, 
“Further Observations of the Angle of Arrival of Microwaves,” Proc. I.R.E., 
34, Nov. 1946, pp. 845-848. H. T. Friis, “Microwave Repeater Research,”’ Bell 
System Tech. J., 27, Part I, ‘Propagation Studies” by A. B. Crawford, Apr. 
1948, pp. 183-246. 


20. E. DeLange, ‘‘Propagation Studies at Microwave Frequencies by Means 
of Very Short Pulses,’’ Bell System Tech. J., 31, Jan. 1952, pp. 91-193. 
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nature of the fading. Some data derived from the continuous recordings 
of signal levels are presented in an appendix. 

The angle-of-arrival observations were made at a frequency of 24,000 
megacycles. The frequency-sweep experiment and the recordings of signal 
levels were made in the 3700 to 4200 megacycle frequency band as were 
the short pulse observations described in the companion paper. 


GENERAL DISCUSSION OF PROPAGATION PHENOMENA 


The map of Fig. 1 shows the location of the experimental transmission 
paths. The path between Crawford Hill and Southard Hill is 17 miles 
long and clears the intervening terrain by 65 feet, approximately one 
Fresnel zone at a frequency of 4000 megacycles. The other path, between 
Crawford Hill and a 100-foot tower on the Murray Hill Laboratory 
property, is 22.8 miles long and has clearance of 280 feet. Fig. 2 shows 
the profiles of these two paths. 

The general characteristics of over-land microwave transmission are 
well known and need be reviewed only briefly. During the daytime 
hours, when the lower atmosphere is thoroughly mixed by rising con- 
vection currents and winds, the signals are normally stable and are near 
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Tig. 1—Map showing location of the experimental transmission path. 


70 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1952 

the free space levels. Also during the winter months, when the humidity 
content of the atmosphere is low, signal variations are usually very small. 
However, on clear summer nights with little or no wind, non-uniform 
distributions of temperature and humidity can create steep dielectric 
constant gradients in the lower atmosphere, thus causing anomalous 
propagation and fading. 

When fading occurred on our experimental transmission paths, an 
alarm circuit connected into the continuously recording equipment was 
arranged to operate when the signal level dropped below a predetermined 
value. This enabled observers to be present during severe fading periods 
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Fig. 2—Profiles of the transmission paths. 


and to seek, by means of the special experiments, to determine the causes 
of the fading. 

Although it has not been possible to provide satisfactory explanations 
for all of the observed fading phenomena, much of the fading (occasions 
when the signals are depressed to levels 15 to 20 decibels or more below 
the normal daytime value) can be explained qualitatively in terms of 
simple ray pictures. Fig. 3 is intended to illustrate some of the observed 
fading phenomena. The case of multiple path transmission, the most 
common cause of fading on either transmission path, is shown in Fig. 3(a). 
Two, three and sometimes more signal components are found to arrive 
- at varlous angles in the vertical plane, usually above the line of sight. 
Wave interference among these components produces fading, the severity 
of which depends upon the relative amplitudes and delays of the com- 
ponents. At these times, different frequencies fade differently and the 
signals received on two vertically spaced antennas also fade differently. 
The use of either frequency or space diversity would be effective in this 
type of fading. . 

A relatively rare type of fading, observed only on the Murray Hill 
path, is believed to be caused by the mechanism illustrated in Fig. 3(b). 
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Here a reflecting layer is situated between the heights of the transmitter 
and receiver. The signal then suffers attenuation due to reflection of 
part of the energy from the direct path. Widely separated frequencies 
are affected in like fashion and the outputs of antennas spaced for 
diversity reception tend to be in agreement although the fine structure 
fading is usually different. 

On neither of the experimental transmission paths is there a rgeular 
ground-reflected component of any consequence. Due to the roughness 
of the ground and the presence of vegetation, the effective reflection coef- 
ficient is of the order of 0.2 for either path. Ground reflections thus play 
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Fig. 3—Possible ray paths involved in severe fading. (a) Multiple path trans- 
mission. (b) Attenuation by reflection from an elevated layer. (c) Abnormal 
water reflection on the Murray Hill-Crawford Hill path. (d) Substandard con- 
ditions on the Southard Hill-Crawford Hill path. 
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no significant part in the fading picture with the exception of the situa- 
tion illustrated in Fig. 3(c). Occasionally on the Murray Hill path, 
conditions of atmospheric refraction are such that a strong signal com- 
ponent is received by virtue of reflection from the water surface of 
Raritan Bay. Under normal conditions, the geometry of the path does 
not permit such a reflection. 

Normally the dielectric constant of the atmosphere decreases with 
height above ground so that the ray path usually has a curvature in 
the same direction as the earth curvature. However, it is possible for the 
dielectric constant of the atmosphere to increase with height above 
ground (sub-standard conditions) so that the ray path has a curvature 
opposite that of the earth. This results in the condition illustrated in 
Fig. 3(d) where the limiting or tangent ray does not reach the receiver 
and only a weak signal is received by virtue of diffraction. Widely sepa- 
rated frequencies and vertically spaced antennas are affected alike as 
regards the average signal level but not the fine structure fading. This 
effect has been observed only on the Southard Hill-Crawford Hill path 
which has small clearance to begin with. It has been observed on several | 
‘ nights in late summer or early autumn after a radiation type ground © 
fog has formed in the late evening and usually persists until the fog is 
dispelled by winds or by the morning sun. . 

There are, of course, times when the transmission conditions are con- 
siderably more complicated than those described above. Some of these 
apparently are due to a combination of the situations illustrated in Fig. 3 
while others may be the result of an atmospheric focussing or trapping 
phenomenon. In addition to the various phenomena just described, which, 
fortunately, occur rather infrequently, there are numberous occasions 
when the signal varies plus and minus a few decibels relative to the free 
space level. It has not been possible actually to demonstrate the mecha- 
nism responsible but it seems most likely that these smaller variations 
are due to non-linear dielectric constant gradients which give the atmos- 
phere the properties of a convergent or divergent lens. 

An important result of the observations made to date is the con- 
viction that the severe fades, signal excursions to levels 30 decibels or 
more below the free space field, were all caused by wave interference. It 
appears that, as the average signal level is depressed by any mechanism, 
it becomes more and more vulnerable to the effects of extra signal 
components of small amplitude that often may be present but go un- 
noticed when the signal is near normal levels. Thus, while the average 
signal level during the conditions illustrated in Figs. 3(b) and 3(d) may 
be no more than 15 to 20 decibels below the normal daytime level, there 
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is usually superimposed a fine structure fading in which short duration 
fades to levels as much as 45 decibels below free space have been ob- 
served. For this reason it is desirable to avoid paths having small clear- 
ance over intervening terrain and also paths which have a permanent 
ground reflection of sufficient magnitude to depress the signal to critical 
levels when, due to variable atmospheric refraction, the direct and re- 
flected components are in phase opposition. 

The following sections describe the angle-of-arrival and frequency- 
sweep experiments on which much of the preceding discussion was 
based. 


ANGLE-OF-ARRIVAL OBSERVATIONS 


A photograph of the Crawford Hill receiving site is shown in Fig. 4. 
The building housing the receiving equipment and the associated an- 
tennas are mounted on a framework which can be rotated on a concrete 
track, permitting investigation of the transmission characteristics of 
either path. The parabolic antennas on the tower are used for continuous © 
recording of 4195 megacycle signals. The long object at the left of the 





Fig. 4—The Crawford Hill receiving site. 
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picture is the metal-lens antenna used for making the angle-of-arrival 
observations. Its half-power beamwidth is 0.12 degree at the operating 
frequency of 24,000 me. The focal length of the lens is 48 feet and its 
feed is located in the little cupola on top of the building. The feed is 
held fixed, while the lens is moved vertically by a motor-driven mecha- 
nism; thus the antenna beam also moves vertically. The antenna scans a 
total angle of two degrees in ten seconds. It is fed by a 24,000-mce radar 
set which is gated to receive only the pulses reflected from a corner 
reflector located at the distant terminal of the transmission path. The 
spot on the radar cathode ray tube moves vertically in synchronism 
with the scanning antenna, and the horizontal deflection is proportional 
to the amplitude of the pulse received from the corner reflector. The 
display thus shows amplitude of the various incoming signal components 
as a function of their angles of arrival. 

The antenna installation on Southard Hill is shown in Fig. 5. At the 
left is the transmitting paraboloid for the 4195-mec continuous wave 
transmitter, the radar corner reflector is in the center, and on the right 
is the horn-reflector antenna used in the frequency-sweep experiments 
described below. Similar equipment is located at the Murray Hill termi- 
nus. The corner reflector is 5.5 feet on a side, and at 24,000-mc has 
sufficient gain to override reflections from other nearby objects, and 
thus becomes easily identifiable on the radar screen. 

The radar oscilloscope for typical propagation conditions is shown in 
Fig. 6. These pictures were obtained by leaving the camera shutter open 
- during the ten-second interval required for the antenna beam to scan 
through the angular range of 2°. All of these representative photographs 
were taken on the Murray Hill-Crawford Hill path although similar 
results were obtained on the Southard Hill-Crawford Hill path with the 
exception of Fig. 6(f). The normal daytime transmission is shown in Fig. 
6(a) to consist of a single path arriving at an angle of —0.2° with respect 
to a fixed reference angle. The horizontal lines represent intervals of 
0.1°, so that changes of 0.05° can be estimated. The other pictures in 
Fig. 6 were all taken during fading conditions. 

Figs. 6(b) and 6(c) are good examples of the multiple-path condition 
shown in Fig. 3(a) in which the individual components are almost equal 
in amplitude and well separated in angle. In Fig. 6(b) there are two 
components arriving at angles of 0.1° and 0.6° above the normal line-of- 
sight while in 6(c) there are three components with angles of 0.05°, 0.35° 
and 0.7° above the normal angle. The position and amplitude of the 
signal components may change radically in a matter of minutes, and 
often there is no component that can be identified as the “normal’’ one. 
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Fig. 6—Representative photographs of angle-of-arrival observations on the 
Murray Hill-Crawford Hill path. (a) Normal day. (b) Two elevated paths. Sept. 
8, 1950; 9:23 p.m. (c) Three elevated paths. Aug. 27, 1950; 1:11 a.m. (d) Multiple 
paths. August 26, 1950; 11:00 p.m. (e) Wide angle “‘fill-in’”’. Aug. 26, 1950; 11:04 
p.m. (f) Abnormal water reflection. Sept. 8, 1950; 11:28 p.m. 





During these multiple-path conditions, the recordings of the 4195-mc 
transmission generally show the broad maxima and sharp minima charac- 
teristic of wave interference. 

Figure 6(d) shows a case in which the various paths are not completely 
separated while Fig. 6(e) (taken four minutes later) shows that energy 
is being received almost without variation over a vertical angle of 0.4°. 
This may represent a number of ray paths which would be separable by 
a narrower-beam antenna, or it may indicate a focussing or trapping 
phenomenon. Often when the type of transmission illustrated by 6(e) 
is present, the recorded 4195-mc signal may be as much as 12 to 15 | 
decibels above the free space levels. 
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Fig. 6(f) illustrates the case of abnormal reflection from the water of 
Raritan Bay on the Murray Hill path as indicated in Fig. 3(c). Here the 
“normal” signal component is arriving at 0.1° above the line-of-sight 
while another component, almost equal in amplitude, is arriving at the 
very bottom of the scan, about 0.8° below thé line-of-sight. It is quite 
probable that there have been times when this component was present 
but was outside the range of the scanning antenna. . 

The mechanisms discussed in connection with Fig. 3(b) and 3(d) can- 
not be demonstrated by photographs such as those just presented al- 
though the angle-of-arrival radar was instrumental in furnishing the 
clues to the phenomena. Due to the two-way attenuation of the radar- 
corner reflector technique, the signal at these times rapidly falls below 
the noise level of the receiver. For the same reason, it is not possible to 
detect the extra signal components of small amplitude which were postu- 
lated to account for the very deep fades sometimes observed under these 
transmission conditions. 


FREQUENCY-SWEEP OBSERVATIONS 


Since most of the fading is due to interference between waves which 
travel over different paths of, presumably, different lengths it was 
realized that the fading was likely to be frequency selective. Just how 
selective would depend on the relative lengths of the individual trans- 
mission paths. The usual methods for determining path length differ- 
ences are to use short pulses, or to sweep the frequency. Since it was 
likely that the path-length differences would be measured in feet rather 
than yards, very short pulses or a wide frequency-sweep were required. 
An oscillator? was available whose frequency could be swept over the 
licensed band of 500 mc between 3700 me and 4200 mc. The frequency- 
sweep experiment was set up on the Murray Hill-Crawford Hill path 
for the summer of 1949. The following summer, the milli-microsecond 
pulse transmission tests described in the companion paper were con- 
ducted over the same path. As might be expected, simultaneous observa- 
tions showed good agreement between the two methods. 

The frequency of the transmitter, located at Murray Hill, is swept 
over a 450-mc band centered at 3950 me at a 60-cycle rate. At the 
receiver, a similar oscillator is used for the beating oscillator except that 
its frequency is swept linearly through the same frequency band in one 

3 This oscillator was developed by M. Ei. Hines and is described in his paper 
published in the Bell System Technical Journal, Vol. 29, Oct. 1950. It uses a 416A 
close-spaced triode in a wave-guide cavity. The frequency is changed by means 


of a plunger which is capacity-coupled to the plate of the tube and which is actu- 
ated by a modified loud speaker unit. 
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second. Since the intermediate frequency amplifier of the receiver is only 
350 ke wide, (centered at 600 kc) narrow pulses are generated each time 
the frequency of the transmitter crosses the frequency to which the 
receiver is tuned. These intermediate frequency pulses are displayed 
vertically on a cathode ray tube. The horizontal trace is synchronized 
with the one-second sweep rate of the beating oscillator. 

The normal daytime frequency-sweep pattern is shown in Fig. 7(a). 
The vertical scale is linear in amplitude and the horizontal scale is almost 
linear in frequency, with frequency decreasing from left to right. Visible 
at the extreme left is the signal used for continuous recording. Since — 
there is only one transmission path involved, the amplitude of the re- 
ceived signal is nearly constant over the 450-mc band. If another signal 
were present which had travelled over a path of different length, the two 
signals would add when the frequency is such that the path length 
difference is an even multiple of half-wavelengths and subtract when 
the path length difference is an odd multiple of half-wavelengths. Simple 
calculation shows that if the path length difference is one foot, the 
frequencies at which the signals add and subtract are separated about 
500 me. Thus the limit of resolution for the frequency-sweep experiment 
is a little more than one foot. 

Photographs taken on a night when the angle-of-arrival radar indicated 
two almost equal components separated about 0.4 degrees in angle are 
shown in Fig. 7(b). The time interval between the two pictures is 30 
seconds, during which the minimum had shifted about 150 mc. The 
pictures can be interpreted as simple two-path transmission with an 
indicated path difference of about two feet and an amplitude ratio of 
0.7 to 1. On this night the minimum shifted back and forth across the 
frequency band—sometimes slowly and sometimes rapidly. At times the 
position of the minimum might remain fixed but its depth would change. 

Photographs taken on a night when there were abnormal reflections 
from the water of Raritan Bay are shown in Fig. 7(c). There are evidently 
two main components with path difference of about six feet, with a small 
third component causing the slight decrease in amplitude of the peaks 
from left to right. These pictures were taken 9 minutes apart, but this 
type of pattern was observed over a period of about three hours on this 
night. 

Usually the frequency sweep patterns are considerably more compli- 
cated than the ones shown so far. Fig. 7(d) shows two photographs which 
indicate that at least three signal components and perhaps more were 
present. The time interval between the two pictures was about 30 
seconds, 
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Fig. 7—Representative frequency-sweep patterns observed on the Murray 
Hill-Crawford Hill path. (Summer 1949.) (a) Normal day. (b) Two components 
with a path difference of two feet. (¢) Two main components with a path differ- 
ence of about six feet plus a small third component. (d) Multiple component 


pattern. 
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SYNTHESIS OF FREQUENCY-SWEEP PATTERNS 


To aid in the interpretation of the complicated frequency sweep pat- 
terns, a computor of the analogue type was built. This apparatus com- 
bines four signal components, three of which are variable in delay and 
amplitude, and presents the result on a cathode ray tube in the same 
form as the actual frequency sweep patterns. Thus a particular pattern 
can be synthesized on the computor and the number of components, 
together with their path differences and relative amplitudes, read directly 
from the computor dials. This is accomplished by generating four 600-ke 
signals, three of which are phase modulated at 60 cycles per second. The 
total phase deviation and relative signal amplitude are variable. The 
four signals are then summed and displayed in vertical deflection on a 
cathode ray tube having a 60-cycle horizontal sweep. 

The synthesis of the patterns of Figures 7(b) and 7(c) are shown in 
Fig. 8. The upper synthesized pattern is simply a combination. of two 
components with relative amplitudes of 0.7 and 1 and a path difference 
of two feet. The lower pattern consists of the reference component with 
unity amplitude, a second component with an amplitude of 0.5 and a 
path difference of 5.7 feet, and a third component with an amplitude of 
0.2 and path difference of 0.8 feet. The similarity between the actual. 
and synthesized patterns is obvious. 








(b) 


Fig. 8—Synthesis of the frequency-sweep patterns of Figs. 7(b) and 7(c). 











Fig. 9—Synthesis of complicated frequency-sweep patterns using four com- 
ponents. See Table I for values of the relative amplitudes and path differences. 
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Examples of attempts to synthesize some more complicated frequency 
sweep patterns, taken on the night of August 2, 1949, are shown in 
Fig. 9. Four components were required in each case, with the path 
differences and delays being summarized in Table I. Although the pic- 
tures all appear very different, in general major changes were required 
only-in component No. 2 to go from one pattern from the next, as Table I 
shows. All the remaining components had to be very carefully trimmed 
in both amplitude and delay to get good synthesis (especially in the case 
of Fig. 9(d), but these changes were relatively small. 


CONCLUDING REMARKS 


The special experiments just described have led to the conclusion, 
expressed earlier, that the severe fading observed on the two test paths 
is the result of multiple-path transmission in which several components 
may be involved. These components may arrive at the receiver at various 
angles up to three quarters of a degree above the normal daytime angle- 
of-arrival and, in the case of abnormal water reflection on the Murray 
Hill path, as much as 0.8 degree below the normal angle. The path 
differences among these components may vary from a fraction of a foot 
to about ten feet. The long-delay components are usually small in ampli- 
tude. . 

In all cases where observations were made during periods of excep- 
tionally high signal levels, say 10 to 15 decibels above free space level, 
the frequency-sweep patterns were substantially flat, suggesting a focus- 
sing or trapping phenomenon. The frequency-sweep patterns were also 
flat on those nights when the signal excursions were only a few decibels 
above and below the normal daytime level. However, the severe fades 


TaBuie | 





Fig. 9 (a) Fig. 9(b) Fig. 9(c) Fig. 9(d) Fig. 9(e) 











Component No. 1| Amplitude 1 1 1 1 1 
(Reference) Path diff. (ft.) 0 0 0 0 0 
Component No, 2 | Amplitude 0.9 i 0.7 1.1 0.4 
Path diff. (ft.) 1.1 0.5 1.7 0.5 1.1 
Component No, 3| Amplitude 0.2 0.1 0.2 0.2 0.2 
Path diff. (ft.) 5.2 5.7 5.6 5.7 5.7 
Component No. 4| Amplitude ; 0.05 0.1 0.15 0.1 0.1 
Path diff. (ft.) 9.2 9.2 11.0 9.3 8.7 

















Time 12:08; AM 12:09 AM } 12:18 AM | 12:24 AM | 12:25 AM 
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Fig. 10—A sequence of angle-of-arrival and frequency-sweep patterns taken 
at ten-second intervals, 


THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1952 


84 







































































oe 166) 12-92 isnonv- . 
TH auOsMWY2 = 1H AVEHOW 









































: ose « ce ¥380190 
oa TH GuOsMVe 





oc- 
S¢e- 


O€- 
Se- 


O¢2- 


ge- 
Se- 


$1381930 


$13981930 


siagioad 


DECIBELS 
t 
3 


DECIBELS 
' 
ro) ro) 


' 
i) 
fe) 


trot 
WW fh 
WOW 


[Fig 11.—Typical recordings of 4195-mc transmission during severe f 
mal daytime (free space) level. (a) Substa 
(b) Multiple path transmission on the Mu 
sion on a night when abnormal water reflec 
(d) Simultaneous recording of the outputs of two simil 























- 

















—— 
3: 


MU 


RRAY HILL - CRAWFORD 
“UPPER ANTENNA. 


HILL & 





- AUGUST 14-15, 1951 


ae 























ndard conditions on the So 



































ading conditions. The zero decibel line represents the nor- 
uthard Hill-Crawford Hill path. Night of Oct. 21-22, 1950. 
rray Hill-Crawford Hill path. Night of Aug. 26-27, 1950. (c) Multiple path transmis- 
tions were present on the Murray Hill-Crawford Hill path. Night of Sept. 8-9, 1950. 
ar antennas; vertical spacing of 30 feet. Night of Aug. 14-15, 1951. 


SHMAVMOUOIN AO DNIGVA AAILOGTAS 


G8 


é 


86 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1952 





wei 
ALU 


% 








i a la a 


SOUTHARD HILL ~ CRAWFORD HILL PATH 








1950 


10 1948 
ce) 
0 
0 


oo oOo 
~ iY 


$1391930 $13819530 $13981930 





-1 
~30 
“1 
-30 
-40 


1948 








10 
n 0 
- 
uw 
M -10 
(s) 
w 
a 

-20 


DECIBELS 





DECIBELS 


Fig. 12—Chart showing the daily fading range observed on both transmission paths for the years 1948, 1949, and 1950. The 
vertical lines terminate on points indicating the highest and lowest signals observed during the twenty-four hour period of noon 
to noon. 





SHUAVMOYOIN FO ONIGVA DALLODTAS 


~] 


88 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1952 


of short duration sometimes observed when the average signal level was 
depressed by the mechanisms of Figs. 3(b) or 38(d) were found to be fre- 
quency selective. 

Some of the studies described in this paper were made with vertically 
polarized waves and some with horizontally polarized waves; at times, 
45° polarization was used. In so far as it was possible to determine, the 
propagation characteristics of both paths were independent of the polari- 
zation used. 

No meteorological soundings were made in connection with this work. 
Considering the rapid changes usually observed with the angle-of-arrival 
and frequency sweep apparatus, it is doubtful that meteorological meas- 
urements made in the usual manner would show much correlation with 
the radio observations except, perhaps, in a general way. The sequence 
of pictures in Fig. 10 is included to show how the angle-of-arrival and 
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Fig. 13—Time distribution curves of the signal levels observed on the Murray 
Hill-Crawford Hill path. Data of 1947, 1948, 1949 and 1950. 
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frequency-sweep patterns change with time. These pictures were taken 
at 10-second intervals. On this occasion there was good correlation be- 
tween the angle-of-arrival and frequency-sweep data. Such was not 
always the case, however, and considering the wide difference in operat- 
ing frequencies, 24,000 mc and 4000 me, instantaneous correlation should 
not necessarily be expected. 

Although all the studies described in this paper were made on the two 
local paths, the results are compatible with propagation measurements 
made by another group in the Laboratories during a survey for the 
transcontinental radio relay system. 
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APPENDIX 


This appendix is included to illustrate some of the characteristics of 
the propagation as shown by the recordings of 4195-me signal levels. 
Fig. 11 is a reproduction of some typical records obtained during severe 
fading periods. Fig. 11(a) is an example of transmission during the sub- 
standard conditions illustrated by the ray diagram of Fig. 3(d). Fig. 
11(b) is typical of multiple-path type fading in which the signal com- 
ponents arrive from elevated angles as shown in Fig. 3(a), while Fig. 11(c) 
was recorded on a night when, for a time, there were abnormal reflections 
from the water of Raritan Bay on the Murray Hill path; see Figs. 3(c), 
6(f) and 7(c). The records of Fig. 11(d) show how the outputs of two 
similar antennas, spaced vertically about 30 feet, differ in regard to the 
deep fades of short duration. 

The chart of Fig. 12 shows how the fading varies with the time of year. 
On this chart, the vertical lines represent the extremes in signal level 
observed during the twenty four hour period from noon to noon. The 
large signal variations are concentrated mainly in the summer months. 

The time distribution of the signal levels recorded on the Murray 
Hill-Crawford Hill path are shown in Fig. 13. Each of the curves is for 
a four-month period: the period of least fading, December, January, 
February and March; the period of most fading, June, July, August 
and September; and the in-between period consisting of April, May, 
October and November. Data obtained in the years 1947, 1948, 1949 and 
1950 are included. Fig. 14 shows similar data for the Southard Hill- 
Crawford Hill path. The hump in the time distribution curve for the 
months of April, May, October and November is due to substandard 
conditions, illustrated by the ray diagram of Fig. 3(d) and the typical 
record of Fig. 11(a), which affected transmission on this path during 
several nights in October, particularly in the years 1947 and 1950. When 
it occurred, this type of transmission usually persisted for a period of 
several hours. 


Propagation Studies at Microwave 
Frequencies by Means of Very 
Short Pulses 


BY O. E. DE LANGE 
(Manuscript received March 27, 1951) 


Microwave pulses with a duration of about 0.003 microseconds were trans- 
mitted over a 22-mile path from Murray Hill, N. J., to Holmdel, N. J., 
in order to determine the effects of the transmission medium upon such 
pulses. During “fading” periods multi-path transmission effects with path 
differences as great as 7 feet were observed, as well as some other effects. 
A microwave frequency of 4000 megacycles was employed. 


INTRODUCTION 


This experiment was set up with two main purposes in view: First, 
as a means of studying microwave propagation, especially with regard to 
multi-path transmission effects and second, to determine the effect of 
a transmission path upon the shapes of very short pulses, particularly 
to learn what restrictions might be imposed upon minimum pulse length 
or spacing between pulses by distortions produced in the transmission 
medium. 

In regard to multi-path transmission the pulse method seems to be 
the most straightforward way of studying such effects. For example, if 
there is transmission by more than one path, and if the pulses are suffi- 
ciently short in comparison to the path length differences involved, then 
there will be received a separate pulse for each path. Under these condi- 
tions the number of paths involved, path length differences and other 
information become directly evident. If pulse duration is too great with 
respect to the path differences involved, the pulses received via the various 
paths will overlap in time and the resultant multi-path effect will be 
pulse distortion rather than reception of individual pulses. This situa- 
tion is much more difficult to analyze. 


TRANSMISSION PATH 


The transmission path is the same as that used by A. B. Crawford for 
microwave propagation studies by means of the frequency sweep method, 
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i.e., the path from Murray Hill, N. J., to Crawford’s Hill (near Holmdel), 
N. J The path length is approximately 22 miles, and is partly over 
water and partly over rough land terrain. The frequency sweep studies 
had indicated that the path differences involved in multi-path fading 
were of the order of one or two to about seven feet. In terms of delay 
times this means differences of about 1 to 7 millimicroseconds. In order 
to resolve the paths when the path differences were only one or two 
feet, we should have liked to have pulses of about 1 millimicrosecond 
duration. Because of the difficulties involved in generating, amplifying 
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Fig. 1—Transmitting equipment. 





and detecting such short pulses, we accepted pulses which, when dis- 
played on our final indicating equipment, had a length of 3 millimicro- 
seconds at half amplitude. (About 6 millimicroseconds at the base.) In 
free space this pulse would be just about 6 feet long at the base. 


TRANSMITTING EQUIPMENT 


The transmitter was mounted on top of a 100-foot tower at Murray 
Hill. As can be seen from Fig. 1, it consisted of a c-w reflex oscillator 
operating at 4000 megacycles, a baseband pulse generator, a modulator, 
or gate, for modulating these pulses on the microwave carrier, a single 
stage traveling-wave amplifier and finally a horn antenna. Approximately 
one watt of power was obtained from the transmitter at the peaks of 
the pulses. The antenna area was 25 square feet and its gain 32 db above 
that of a dipole. A.pulse repetition frequency of 10 mc was employed. 


1A.B. Crawford and W. C. Jakes, Jr., “Selective Fading of Microwaves,”’’ 
Bell System Tech. J., 31, Jan. 1952, pp. 68-90. 
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RECEIVING EQUIPMENT 


The receiving antenna, a large horn, was-‘ mounted between two poles 
guyed for support. It had an aperture of about 90 square feet and a gain 
of approximately 38 db over a dipole. The receiver circuit is shown in Fig. 
2. About 60 db of gain at 4000 me was provided by either two or three 
stages of traveling-wave tube amplifier depending upon the gain of the 
particular tubes used. It was necessary to provide very good shielding 
and also careful filtering of all power leads to eliminate the tendency for 
this amplifier to sing. he amplifier fed two crystal detectors through a 
hybrid tee junction. Each detector employed a silicon crystal of the 
IN23B type. 

Two indicator circuits are shown in Fig. 2. These circuits are very 
similar except that one employed a vertical amplifier coupled to a Dumont 
5XP2 CRO tube, whereas in the second the baseband output of the 
crystal was fed directly onto the deflection system of a traveling-wave 
type of CRO tube. The latter CRO tube, which has been described by 
J. R. Pierce in the November, 1949, issue of Electronics, has a very high 
deflection sensitivity and is used with a microscope to enlarge its trace; 
hence, no amplification was required between it and its driving crystal. 
The deflection system of this tube has a bandwidth of 500 to 1000 mc. 
The micro-oscilloscope was provided primarily for photographing pulses 
by means of a 35-mm camera attached to the microscope. (Exposure 
time was 5 to 15 seconds. The time recorded for each picture corresponds 
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Fig. 2—Receiving equipment. 
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to that at the end of exposure.) A second microscope made it possible 
to view and to photograph the screen of the tube simultaneously. The 
general procedure was to observe continuously during periods of dis- 
turbed transmission, taking pictures at regular intervals of 5 to 10 min- 
utes. When conditions were seen to be changing rapidly, pictures were . 
taken much more frequently. The large oscilloscope with its vertical 
amplifier had a bandwidth of about 150 mc and hence caused some 
deterioration of the pulse. It, however, was less tiring than the small 
scope, especially for long periods of observation and was watched to 
follow the general trend of events. It was capable of resolving the pulses 
resulting from two-path transmission when the path differences were 
large. 

The sweep circuits for the two indicator oscilloscopes were practically 
identical. The horizontal sweep voltage for each consisted of the linear 
portion of a sine wave which was generated by a c-w oscillator operating 
at one third of the pulse repetition frequency of 10 mc. Each oscillator 
was synchronized with the incoming pulses by means of a 10-mc voltage 
derived by amplifying the pulse energy through a narrow band amplifier. 
This circuit provided very satisfactory synchronization even during the 
times when signal amplitude was so low as to produce a very poor signal- 
to-noise ratio. Timing markers were provided on each roll of film by 
periodically photographing a series of pulses spaced by an interval of 9 
millimicroseconds. 


RESULTS OF THE EXPERIMENT 


The picture at the left of Fig. 3 shows the transmitted pulse. The 
right-hand picture shows the received pulse under what were considered 
to be normal transmission conditions. It is seen that, except for the addi- 
tion of noise and widening of the pulse due to passage through the ampli- 
fiers and other equipment, the pulse shape is unaffected. The time cali- 
bration on this and the following photographs are in millimicroseconds, 
each mark representing one millimicrosecond (0.001 us). 

During the summer of 1950, when this experiment was in progress, 





Fig.3—(Left) transmitted pulse (right) received pulse—normal transmission. 
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Fig. 4—Received pulses—disturbed transmission. 


there was comparatively little fading over the path in question at micro- 
wave frequencies. There were, however, a few nights of considerable 
activity and some interesting results were obtained. . 
The series of pictures on Fig. 4 show one good example of multi-path 
transmission where the path length difference was great enough to pro- 
duce complete resolution of pulses. At 11:05-20 there are two pulses, 
each 7 millimicroseconds wide at the base and with their peaks just 7 
millimicroseconds apart; in other words, the path difference was just 
sufficient to produce two pulses with no overlap. The pulse at the left 
is presumably coming by the main path and that at the right from some 
second path resulting from bending of the rays caused by atmospheric 
effects. At 11:07-20 the second path appears to have shortened, resulting 
in a path difference of only about 5 millimicroseconds. This may actually 
have been due to a change of length of the second path or it may have 
been due to distortion of the second pulse by energy coming by way of a 
third path. The pictures taken at 11:08-20 and 11:10-50 show evidence 
of transmission by a third path. In the first of these, for example, the 
width of the disturbance at the base line indicates the presence of the two 
original pulses spaced 7 millimicroseconds apart but the midpoint of the 
two no longer falls to the base line as was the case in the first picture. 
This could be accounted for by the presence of a third pulse coming over 
a path whose length was somewhere between that of the other two. Con- 
ditions obtaining at 11:10-50 could also be accounted for by the presence 
of pulses from three paths, that is, energy coming by way of a third path 
might cancel part of one pulse and at the same time add to the other. 
This could account for the fact that the spurious pulse is larger than the 
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Fig. 5—Received pulses—disturbed transmission. 


normal one. It is also possible that more than three paths were involved. 
On a number of other occasions the pulse coming by way of the second 
path appeared to be of greater amplitude than the one coming by the 
main path. This same effect has been observed by Mr. Crawford and his 
colleagues on the angle of arrival equipment. 

Information obtained from the above set of pictures shows that for 
a time-division multiplex system using the length of pulse used here 
(7 millimicroseconds at the base) and operating over this path, pulses 
would have to be spaced a minimum of about 14 millimicroseconds apart 
if it were desired to avoid distortion at all times. If very much shorter 
pulses were used the spacing might be reduced to 9 or 10 millimicro- 
seconds. However, the 7-foot path difference indicated by these pictures 
is about the maximum ever observed and occurs rather infrequently so 
that if somewhat closer spacings were employed troubles would result 
only a small percentage of the time. 

The next series of pictures, Fig. 5, taken July 8, show an example 
of a more common type of multi-path transmission. Here the path dif- 
ference is apparently less than for the last series. At 11:15 there are two 
distinct pulses with an apparent path difference of about six feet (6 milli- 
microseconds) if judged from the spacing between the peaks of the pulses. 
However, from the length of the disturbance at the base line, which we 
consider a better criterion, the path difference was more nearly four feet 
At 11:22 distortion of the trailing edge of the pulse was the only indica- 
tion of a second path. For the pictures taken at 11:29 and 11:44 the 
path difference is sufficiently small that there is almost complete can- 
cellation of pulses, only the leading portion of each pulse being present. 
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On the 11:44 picture there is just a trace of a second pulse. The next 
set of pictures (Fig. 6) were taken a little over an hour later on the same 
night and show about the same conditions, that is, pulse amplitude and 
shape change and other evidence of the presence of a second pulse delayed 
about 2 to 3 millimicroseconds. 

On the night of October 2, fading, which was apparently due to trans- 
mission by way of two paths with little path difference, was observed. 
Some of the results are shown on Fig. 7. At 7:49 two distinct pulses are 
evident, there being 6 millimicroseconds between their peaks. One might 
conclude from this that there was a second path about 6 feet longer than 
the main path but the total length of the disturbance along the base line 
and the shapes of the pulses indicate that the actual path difference was 
about 2 to 3 feet. 

Apparently we had here two pulses of r-f energy overlapping in time 
and involving a large number of frequencies. These pulses are capable 
of interfering with each other in a rather complicated manner, it being 
possible for some frequencies to add and others to cancel at the same time, 
depending upon their relative phases. Phase relationships of course de- 
pend upon frequency and path length differences. As a result pulses may 
be distorted and have their peaks shifted about by a considerable amount. 
We must, therefore, realize that the first picture of Fig. 7 does not really 
represent two distinct pulses as appears to be the case, but actually shows 
the resultant interference pattern of two overlapping pulses. Since 
a change of path difference of only about one and one-half inches is 
enough to produce a 180° change in relative phase at 4000 me, it is not 





Fig. 6—Received pulses—disturbed transmission. 
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Fig. 7—Received pulses—disturbed transmission. 


at all surprising that pulse shapes and amplitudes change very rapidly 
at times. . 

Looking again at the photograph, Fig. 7, we see that at 7:54-05 there 
was a complete fade as far as our system is concerned. To produce this 
degree of cancellation the. path difference must have been very small 
though still sufficient to give a relative phase angle of 180° at the radio 
frequencies involved. At 8:00-15 and 8:08-00 pulse distortion is the most 
noticeable effect of the “fading,” the pulses being considerably shorter 
than their normal value. Pictures, not shown here, taken between 7:54 
and 8:08 show definite evidence of two-path transmission with a path 
difference of 2 to 3 feet; therefore the pulses of 8:00-15 and 8:08-00 
are probably also the result of two-path effects. 

The first two pictures of Fig. 8 show another form of pulse distortion 
observed on a number of occasions. Here the pulse is flattened out on 
top probably due to energy coming in over a second path differing in 
length by only one or two feet from the main path. Each time this type 
of pulse was observed a check was made to be sure that the flattening 
was not due to overload in our equipment. The pictures presented up to 
now have all shown comparatively slow changes of conditions. Very 
rapid changes were, however, quite common. In many cases pulse shape 
or amplitude changed considerably during the 5 to 15 second exposure 
time ordinarily used. The picture taken at 2:20-45 A.M. on August 27 
is one example of such a rapid change, there being tio definite sets of 
conditions shown on the one photograph. The remaining picture on Fig. 
8 shows the pulses used for obtaining time calibration of the system. 
These pulses were spaced 9 millimicroseconds apart and by adjust- 
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ing sweep expansion so that succeeding pulses fell on proper parts of 
the scale and by keeping this expansion constant, it was possible to ob- 
tain a calibration. 


TWO-PATH SIMULATOR 


As an aid to interpreting the results obtained from the above experi- 
ment, particularly when the two pulses overlap and interfere, a circuit 
was set up in the laboratory to simulate two-path transmission. The 
equipment, as shown on Fig. 9, consisted of a wave guide hybrid junction 
with the r-f pulse energy being fed into the E plane arm. To each side 
arm was connected a variable attenuator in series with a few feet of wave 
guide fitted with a short circuiting plunger. Waves reflected from these 
two plungers recombine in the H plane arm where the detector is located. 
There are two separate paths through the hybrid as follows: (1) Input, 
side arm A, reflecting plunger A, side arm A to output. (2) Input, side 
arm B, reflecting plunger B, side arm B to output. By adjusting the 
attenuator in either branch the amplitude of the signal transmitted by 
way of that branch could be adjusted. In the same way by adjusting 
the position of the reflecting plunger in either branch the distance 
traveled by a signal in traversing that branch could be varied. 

If the path lengths were made the same and the amplitudes adjusted 
to be equal there would be perfect cancellation due to a phase turn-over 
in the hybrid junction and hence no output from the detector. If one 
plunger were now left fixed in the above position and the other moved by 
a quarter wavelength (to produce a total shift of half wavelength or 180°) 





Fig. 8—Received pulses and calibrating pulses. 
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the output signal would be at maximum amplitude due to addition of 
energy coming from the two paths. This amplitude is, of course, twice 
that of the signal from one branch alone. In the experiment the plunger 
in one branch was left fixed and the attenuator in that branch left set 
at zero. The path through this branch then represented the normal 
transmission path for an actual system. The path through the other 
_ branch could be made to correspond to spurious paths having different 
amounts of delay and attenuation simply by adjusting the position of the 
reflecting plunger and the setting of the attenuator. A series of photo- 
graphs were taken of pulses resulting from these different amounts of 
delay and attenuation. 

The first three pictures of Fig. 10 were taken with the path lengths 
exactly equal. When the amplitudes were also equal there was complete 
cancellation. As the signal in one branch was attenuated the amplitude 
of the resultant pulse increased until it became equal to that of the 
original pulse as shown in the third picture. Increasing one path by one- 
half wavelength brought the signals from the two branches into phase 
and they added up to double amplitude as seen in the fourth picture. It 
should be pointed out that although in our experiment we changed delay 
by 0.36 millimicroseconds in going from the first minimum to the first 
maximum, in free space a change of delay of only 0.125 millimicroseconds 
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Fig. 9—Apparatus to simulate two-path transmission. 
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Fig. 10—Simulated two-path transmission. 


would be required. The discrepancy lies in the large ratio between the 
phase velocity and group velocity in the wave guide used whereas in 
free space this ratio is, of course, equal to unity. In free space the amount 
of delay required to go from a maximum to a minimum signal corresponds 
to a change of path difference of only about one and one-half inches. 
With only this slight shift required to change conditions from those shown 
by the first picture of Fig. 10 to those shown by the last, it is not at all 
surprising that the received signal is very unstable during time of multi- 
path transmission. , 

Fig. 11 shows the effect of changing relative phases in 90° steps while 
keeping the amplitudes equal to each other. It is seen that even with the 
carriers in direct opposition cancellation is far from complete due to the 





Fig. 11—Simulated two-path transmission. 
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Fig. 12—Simulated two-path transmission. 


relative delay between the two component pulses. Flat topped pulses 
seem to be characteristic of conditions where the two carriers are in phase 
quadrature and about equal in amplitude. 

Fig. 12 shows a set of conditions with a constant delay difference of 2 
millimicroseconds (corresponding to a path difference in free space of 
about 2 feet). For the pulses shown on Fig. 13 there was a constant delay 
difference of 7.34 millimicroseconds, enough to provide complete separa- 
tion of the pulses. The carriers were in phase opposition but with this 
amount of separation there is no overlap of pulses and the results would 
have been the same if the phase had been aiding. Any increase of delay 





Fig. 13—Simulated two-path transmission. 
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beyond this point results only in moving the pulses farther apart and has 
no effect upon pulse shape or amplitude. 

The experimental set up just described proved to be somewhat un- 
satisfactory since it was not possible with it to produce phase opposition 
between the two carriers without having zero delay difference between the 
two paths or a difference of at least 0.66 millimicroseconds. For the length 
of pulse used this latter amount of delay difference is sufficient to prevent 
anything like complete cancellation of pulses. In fact the amplitudes of 
the two resultant peaks are only about 12 db below the peak amplitude 
of the original pulse. From this we know that for the natural path any 
fade which appeared to be complete must have resulted from path dif- 
ferences of less than 0.66 feet, in fact from differences of less than about 
one-half foot. 


SUMMARY 


The pulse experiment results indicate that over one particular path 
at least there is, at times, transmission of microwaves by at least two, and 
probably more than two, paths. Path differences involved are from a 
fraction of a foot up to about seven feet, differences of less than about 
three feet being the most common. These results agree with those ob- 
tained by other methods. These multi-path effects result in bad distortion 
of very short pulses and even in the presence of entirely separate spurious 
pulses. These effects put a definite lower limit on pulselength and spacing 
between pulses in a pulse transmission system. The limit depends upon 
the amount of distortion which can be tolerated and also upon the per- 
centage of time such distortion can be accepted. No statistical data were 
recorded. 

With the laboratory equipment for simulating transmission over two 
paths, many of the waveforms obtained over the natural path could be 
duplicated. There were times, however, when the waveforms received by 
way of the natural path were too complicated to be explained by trans- 
mission by as few at two paths. 
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Properties of Ionic Bombarded Silicon 


BY RUSSELL S. OHL 
(Manuscript received August 23, 1951) 


This paper deals with a new and very interesting technique by which the 
properties ot silicon surfaces are altered very materially by bombardment 
with ions of such gases as hydrogen, helawm, nitrogen and argon. The change 
in rectifying properties has been of special interest but there have been con- 
sidered also changes in the structural features of the material itself. The 
effects of bombardment on the rectifying properties are allustrated by a series 
of characteristic curves systematically arranged to bring out the effects of 
the several variables of experiment such, for example, as ion velocity, in- 
tensily of bombarding current, length of time of bombardment, kind of gas, 
and the temperature of the specimen during bombardment. The effect of 
bombardment on materials contaminated with impurities is also illustrated. 
It is of particular practical importance that silicon contaminated with boron 
to the point where it shows relatively little rectification can be modified by 
bombardment to make it even better than most wnbombarded materials. 


Some years ago, the writer discovered that the electrical properties of 
silicon surfaces could be greatly modified by bombardment with positive 
ions. The ions in question were generated in a low pressure discharge in 
some gas, like hydrogen, helium or nitrogen, and after passing through a 
perforated cathode were accelerated to a suitable velocity before imping- 
ing on the surface to be treated. This scheme may be contrasted with 
other methods subsequently reported for treating germanium! in which 
high-velocity ions were derived from radioactive sources. Preliminary 
results of the present research were described in a paper entitled Silicon 
Transistors, by W. J. Pietenpol and the writer, presented at an Elec- 
tronics Conference held at the University of Michigan, June 22, 1950. 
Since that time exploration has continued with a view both to learning 
about basic principles and about possible practical applications. 

Editorial Note—Since the resurgence of interest in point-contact rectifiers, 
considerable research has been carried on into the characteristics of silicon and 
germanium. The author of this paper.was a pioneer in this new field of study, as 
evidenced, for example, by Patent No. 2,378,944, applied for on July 26, 1939, and 
Patent No. 2,402,839, applied for on March 27, 1941. More recent work has been 
described in a large number of text books and technical papers such as Electrons 
and Holes in Semi-Conductors by William Shockley, D. Van Nostrand, 1950, and 
numerous papers by Lark-Horowitz published mostly in Physical Review. The 


work described in the accompanying paper is a continuation of this long research. 
1 Brattain and Pearson, Phys. Rev., 80, Dec. 1950. 
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The present paper gives the results of some more recent experiments 
made with improved equipment. Also described briefly are some related 
experiments in which silicon is bombarded with alpha particles derived 
from radioactive polonium. The overall results of this work indicate 
rather clearly that with suitable variations of bombarding voltage, target 
temperature and time of exposure as well as impurity content in the base 
material, it is possible to prepare to specification silicon surfaces having 
a wide range of properties. From the materials so treated it has been 
possible to construct improved forms of signal rectifiers, harmonic gener- 
ators, transistors, modulators, gating devices and also photo-electric 
cells. It is particularly significant that the voltage range over which these 
newer devices can be operated has been greatly extended, thus making 
them useful in places not previously regarded as possible. Since these 
new surfaces appear to be readily reproducible in large numbers and 
since they are electrically tough, chemically stable and show no unsatis- 
factory temperature or aging effects, it would appear that bombarding 
techniques should have considerable practical value. 

This paper is concerned mainly with the practical aspect of ion bom- 
bardment of silicon, namely its effect on the voltage current character- 
istics at low frequencies. Equally important, perhaps, are its theoretical 
aspects, particularly with regard to the interpretation of the rather 
pronounced changes in the properties in light of presently-accepted views 
of solid-state physics. These aspects are not covered in this paper. 


METHOD 


The bombardment process referred to above consists of exposing the 
silicon surface to ions that have previously been accelerated to energies 
in the range from about 100 electron volts to about 30 kilo-electron-volts. 
A recent setup is illustrated in Fig. 1. The electrons from a tungsten 
cathode are accelerated toward a grid which is at a positive potential with 
respect to the cathode. Many of the electrons pass a short distance 
beyond the grid and return for ultimate capture. Ionization due mainly 
to the impacts of electrons with gas molecules takes place in this turn- 
around region, producing amongst other things positive particles. Elec- 
trodes are so proportioned that this ionization is fairly uniform over the 
grid area. 

The silicon specimen to be bombarded is made negative with respect 
to the filament. This accelerates the positively charged particles toward 
the target. The latter rests on a graphite plate heated by a coil below, 
carrying high-frequency currents. A thermocouple with suitable con- 
nections to the exterior makes possible an adequate measurement of 
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temperature. The apparatus will accommodate circular surfaces as large 
as 15 inches diameter. The gases from which ions are derived are ad- 
mitted through the gas inlet. Thus far experiments have been made 
with hydrogen, helium, nitrogen and argon. The bombarding voltages 
have as already noted, been varied from 100 to 30,000 volts and the 
surface temperatures have been varied from about 20°C to 400°C. The 
effects of these several variables will be discussed more fully below. 


SAMPLE PREPARATION 


. The material to be bombarded is usually prepared in batches of about 
300 grams in fused silica crucibles roughly cylindrical in shape.? After 
solidifying, the cast is ground to a cylinder approximately 14 inches 
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Fig. 1—Bombardment apparatus. 
2 Seaff, Theuerer, and Schumacher, A.I.M.M.E., 185, pp. 382-392, 149. 
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diameter. This has the effect of removing some of the contaminating 
impurities derived from the crucible as well as providing samples of 
convenient size. This 14-inch cylinder is then sliced transversely into 
thin wafers which subsequently are polished on one side. Except as other- 
wise noted the material covered by this paper had an impurity content 
of about 0.1 per cent. The exception will be found in the data of column 
(a) of Fig. 8. 


BOMBARDMENT PROCEDURE 


The wafer, as prepared above, is placed in the bombarding apparatus 
with the rough face contacting the graphite support. The vacuum cham- 
ber is sealed by placing the ion generator in position and the whole 
assembly is evacuated. The sample is then heated to the proper tempera- 
ture and the desired kind of gas is admitted, the pressure being estimated 
from the ion current. When stable conditions prevail, the accelerating 
voltage is applied to the target and the bombardment is carried out for 
the proper length of time. A convenient current density is 5 micro- 
amperes per square centimeter of target area. The target area of our 
present apparatus, including the silicon and a portion of the graphite 
support, being 20 square centimeters, the ion current is generally around 
100 microamperes. The dosage is sometimes specified in microcoulombs. 

After bombardment, the sample is removed from the apparatus and 
the rough surface is covered with a thin layer of evaporated rhodium. 
For most of the tests outlined below the 15-inch diameter wafers were 
cut into $-inch squares, a size convenient for testing. 


GRAPHICAL REPRESENTATION 


In considering the merits of non-linear materials such as silicon, per- 
haps the simplest and most useful characteristic is the voltage-current 
relation. If this is plotted to a linear scale, it results in a smooth curve 
of the general form shown in Fig. 2a. Specific curves obtained in practice 
may depart widely from that shown but in general, all may be regarded 
as made up of two semiparabolas, one in the first quadrant and one in 
the third, joined by a nearly horizontal straight line. For present pur- 
poses, we shall further simplify this idealized characteristic by consider- 
ing it as made up of three straight lines. The first, AB, is associated with 
the reverse voltage current characteristic. The third, CD, is associated 
with the forward voltage current characteristic. These two character- 
istics are joined by the nearly horizontal line, BC. The slopes of these 
three lines correspond to resistances. The section BC for example, corre- 
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sponds to a region in which resistance is very high. The points B and C 
are particularly important for they represent points of inflection where 
the resistance undergoes rapid change and the material is departing 
most markedly from Ohm’s Law. Ideally they should be sharp but in 
practice there is usually considerable curvature. Though either inflection 
point could presumably be used in detection processes, the point to the 
right of the origin is for practical reasons, usually preferred. Point B 
defines a voltage Hz at which substantial backward currents flow. It is 
referred to simply as the reverse voltage. In a similar way, point C defines 
a forward voltage Er. The distance between B and C (Hz + Fy) will be 
referred to as the inflection interval. The difference in these quantities 
(i, — Ey) is also of interest. One-half of this voltage difference is re- 
ferred to as the self-biasing voltage. It is a significant quantity readily 
measured in practice by noting the d-c voltage across a large condenser 
placed in series with the crystal and a supply of 60 cycles AC. For de- 
tectors, point C should preferably be close to the origin and Ey should 





Loe I 
Fig. 2—Idealized characteristic curves. 
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be small. For certain kinds of voltage limiters, Hy should be large. In 
either case the inflection interval should be large. 

In an alternate graphical representation, see Fig. 2b, voltage-current 
data are plotted to a log-log scale. This form of representation is of 
particular value when large ranges of data are to be shown. It is also of 
value in determining the resistance (2,) at small voltages. Corresponding 
points on the two curves shown in Fig. 2 are identifiable by the letters A 
B, C, and D. Curves of both kinds are used interchangeably to show the 
effects of the several variables of the experiment. 


EFFECT OF CONTACT PRESSURE 


In point contact rectifiers,? pressure is of considerable importance. 
Usually the best pressure is a compromise between good electrical charac- 
teristics, usually obtainable only with light pressures, and good stability 
usually obtainable with higher pressures. Experiments have been per- 
formed with a range of contact pressures both on bombarded and unbom- 
barded materials. In general, the results are highly variable, particularly 
in the case of unbombarded material. From this wide range of data, 
however, two characteristics have been selected that may be regarded 
as typical for 10-gram and 60-gram pressure. They are shown in Fig. 3 
for silicon taken from nearby portions of the same sample. Significant 
points on these several curves may be compared with their idealized 
counterparts shown in Fig. 2. Although the samples chosen show some- 
what more than the usual intrinsic resistance typical of p-type silicon, 
the effects of contact pressure are nevertheless regarded as representative. 
As indicated in Figs. 3a and 3b, the effect of increased contact pressure,' 
particularly in the case of unbombarded material, is of reducing the low 
voltage resistance, R,, see Fig. 2b. The more desirable higher resistance 
is obtainable only with light contact, a condition unfavorable for high 
mechanical stability. In the case of bombarded material, the effect of 
contact pressure is less important. Thus it is possible in this case to 
incorporate in the design higher contact pressures and obtain thereby 
higher stabilities. For purposes of this paper a contact force of 10 grams 
has been accepted as standard. 

In addition to showing the effect of contact pressure, Fig. 3 shows some 
overall effects of bombardment. It will be noted, for example, that the 
effect of bombardment, see Fig. 3b, has been that of shifting the plots 
of Fig. 3a to the left by several orders of magnitude. Thus the resistance 
(R,) is increased by a factor of more than 10,000. It is to be noted also 


3 Scaff and Ohl, Bell System Tech. J., 26, Jan. 1947. 
* Realy contact force. 
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Fig. 3—Characteristic curves. (a) and (c) unbombarded silicon, (b) and (d) 
Silicon bombarded with 30-kv helium ions. 
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from a comparison of Fig. 3a with Fig. 3b that at the one volt level, the 
ratio of forward to reverse currents for the unbombarded case is about 
twenty, whereas that for the bombarded case, is more than 10,000. At 
other levels the difference is even greater. Referring particularly to Figs. 
3c and 38d, it will be seen that one effect of bombardment is that of 
separating the two significant points of inflection B and C. That is, the 
inflection interval has been notably increased. This increase is the result 
of a small increase in the forward voltage and a very substantial increase 
in the reverse voltage. 


EFFECT OF TYPE OF GAS 


Four high purity gases were tested as ion sources, namely, hydrogen, 
helium, nitrogen and argon, having atomic weights respectively of 1, 4, 
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Fig. 4—Characteristics sbowine effect of various gases all with a bombarding 
potential of 30 kv. 
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14 and 40. While all four gases worked very well, helium was the easiest 
to handle. In the course of the tests, identically prepared samples each 
z-Inch square, taken from a high-purity silicon melt, were bombarded 
with ions formed in the particular gas under test. Particular conditions 
known to be good for producing good rectifier units were adopted as 
standard for these tests. They corresponded to a total bombarding charge 
of 600 microcoulombs per sq. cm., a surface temperature of 395°C a 
contact force of about 10 grams and a bombarding potential of 30 kv. 

That the effect of bombardment varies with different gases is seen at 
a glance from the characteristic current-voltage curves shown in Fig. 4. 
Figs. 4b to 4e in particular indicate that as compared with an untreated 
sample, Fig. ta, the effect of bombardment is in general that already 
noted of separating the two significant points of inflection, B and C. A 
rather substantial increase in the forward voltage appears in the case 
of argon as compared with hydrogen, helium and nitrogen. In contrast 
with a small increase in the forward voltage resulting from the bom- 
bardment of helium, there is a very substantial increase in the backward 
voltage. Though substantial for all four gases, the effect of bombardment 
is largest for helium with progressively smaller effects noted respectively 
for argon, hydrogen and nitrogen. A particular characteristic of helium 
bombardment, as compared with argon, not readily appreciated from a 
linear scale, is shown in Fig. 4f. It will be noted that at the one volt level, 
the ratio of forward to reverse current for helium is about 130 whereas 
for argon it is about 25. At other levels the difference is even greater. 
At the moment helium is regarded as a preferred source of ions. 

The log-log current-voltage curves show as before that the lowest 
voltage at which substantial forward currents flow occurs in helium, while 
the highest forward voltages occur for argon. In a similar way the 
voltages for substantial reverse currents are highest for helium and 
lowest for argon. The sharp break in the reverse current characteristic, 
evident in these cases, has been observed so generally that it is now 
accepted as typical of bombarded surfaces. 


EFFECTS OF TEMPERATURE 


Investigations have been made of the properties of silicon surfaces 
as affected by the temperature at which bombardment was carried out. 
This has been done not only for surfaces used as rectifiers but surfaces 
used as transistors and as photo-electric cells as well. In the case of recti- 
fiers, a procedure was adopted similar to that used in the previous tests. 
Measurements were made at five different temperatures ranging from 
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Fig. 5—Characteristics showing effects of voltage and temperature variation. 


25°C to 395°C each at accelerating voltages of 1 kv, 3 kv, 10 kv and 30 
kv using helium gas. The data so obtained were useful not only for study- 
ing the effect of temperature but useful in the studies of the effect of 
ion velocity as well. The latter will be discussed in the following section. 

The results of the above measurements are plotted in Fig. 5. They 
are further summarized in Fig. 6a. The latter figure, in particular, indi- 
cates that as rectifiers, there is little choice of surface temperature be- 
tween about 250°C and 400°C. It has been found, however, that for 
temperatures below about 250°C the point contact seems to be more 
vulnerable to electrical shock. 
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EFFECT OF ION VELOCITY 


The effect of ion velocity (bombarding voltage) has been investigated 
for several types of silicon. The effects vary with the different types. 
Typical results are those given in Fig. 5 already referred to. It will be 
noted from a comparison of the data for a particular temperature, say 
300°C, that the principal effect of increased ion velocity is that of increas- 
ing the reverse voltage. Values of these reverse voltages /, and also the 
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corresponding forward voltages Hy have been scaled from the above 
drawings and have been tabulated below. Since they vary only slightly 
with surface temperature, only their mean values are regarded as signifi- 
cant. Mean values of sz are plotted in Fig. 6b. 

In order to isolate further the effects of bombardment we have sub- 
tracted from the mean values of Ez value of KH, for untreated silicon. 
Thus the curve marked E; represents the improvement in backward 
voltage that has accrued from bombardment alone. This is also tabu- 
lated as the mean Ey in Table I. 


EFFECT OF TOTAL CHARGE 


Tests havi been made to determine the effect of time of bombard- 
ment on the rectifying properties of silicon surfaces. In these tests, 
specimens taken from neighboring regions of the same melt were ex- 
posed for progressively longer periods all at the same bombarding po- 
tential of 30 kv and the same rate and density of application, 5 micro- 
amperes per square centimeter. Representative current-voltage charac- 
teristics are plotted in Figs. 7a and 7b for two neighboring regions. 
The results are summarized in Fig. 7c. The latter show a rather rapid 
improvement of back voltage Lz with total charge up to perhaps 50 
microcoulombs per square centimeter.’ Thereafter the improvement is 
small. For purposes of comparison, there is plotted as a vertical line a 
value of bombarding charge that would account theoretically for one 
positive ion in each unit crystallographic cell on the surface layer. This 
suggests that when all surface cells have been penetrated by a single 
ion, no marked increase in back voltage can be effected. 


5 Specifying results in terms of microcoulombs implies that bombarding ef- 
fects are independent of the rate of application. This is known to be true only 
between factor limits of 4 and 2 of the bombarding current. 
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Fig. 7—Effect of time and bombarding charge on back voltage. 
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EFFECT OF MATERIAL COMPOSITION 


Thus far discussions have centered around a single type of high purity 
material that was regarded as representative. It is of interest to examine 
the effect on other materials particularly those in which impurities have 
been added. For this purpose tests were made on comparable samples 
from four sources all bombarded for two minutes with 5 microamperes 
of current and each with five representative bombarding voltages. The 
results are illustrated by the curves shown in Fig. 8. The four columns 
correspond to progressively higher percentages of impurities beginning 
with (a) on the left as a material having an impurity content believed 
to be less than 0.01 per cent. The impurity content of (b) is not known 
accurately except that it lies between (a) and (c). The material repre- 
sented by column (c) was produced by adding 0.02 per cent boron® to 
a material illustrated in column (a). The last column (d) was produced 
by adding 0.1 per cent boron to the material illustrated in column (b). 

It is to be noted from Fig. 8 that marked changes in the voltage- 
current characteristic may be effected by bombardment for all degrees 
of the impurity content shown. It is especially interesting that in columns 
(c) and (d) corresponding to materials contaminated with boron to the 
point where nonlinearity is almost absent, rectification can not only 
be restored but indeed the product may be made better than the best 
unbombarded material. 


EFFECTS OF ALPHA-PARTICLE BOMBARDMENT 


The close relationship between helium ions such as generated above 
and alpha particles such as emanate from radioactive materials suggests 
that the latter may be used for the bombardment of silicon surfaces. 
A. few experiments of this kind have been made with results that are 
not only interesting but possibly useful. For these tests, four sources of 
alpha particles were obtained. They consisted of 3%; inch square pieces 
of nickel on which had been plated a thin coating of polonium followed 
by a covering of gold. The initial strength was 4 millicuries per square 
centimeter. The half-life of polonium is 140 days. 

The process of bombardment consisted simply of placing the polished 
surface of a standard silicon square against the layer of gold and examin- 
ing the same periodically. Tests of four samples were carried out simul- 
taneously. The results are given in Table II. The data for Sample No. 1 
departs so markedly from the mean that it may be disregarded. Since 


6 Boron is a particularly active agent in effecting changes in the properties 
of silicon. 
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Fig. 8—Effect of impurity content. (a) hyper-purity silicon, (b) high- -purity 
silicon, (¢) hyper-purity silicon plus 0.02 per cent boron, (d) high- -purity silicon 
plus 0. 1 per cent boron. 
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the data for self-bias are the result of a direct measurement while those 
for forward and reverse voltage are transcribed from a cathode-ray 
plot, they are perhaps the most significant. Samples 1 to 3 represent 
specimens of increasing degrees of purity. 

These alpha-particle bombardment experiments indicate rather defi- 
nitely that results may be obtained similar to those obtained from 
bombardment with gaseous ions and, like the ion bombardment, they 
tend to produce high resistance surfaces. 

















TABLE II 
Sample Number 
Days of Exposure 

1 2 3a 3b 

4 Self Bias-Volts <1 6 180 60 
4 Reverse Voltage <1 20 380 130 
4 Forward Voltage <0.5 | <0.5; <1 <1 
13 Self Bias-Volts <1 160 200 190 
13 Reverse Voltage 1 360 440 420 
13 Forward Voltage <0.5| <0.5|] <1 <1 
39 Self Bias-Volts <1 60 130 100 
39 Reverse Voltage 0.5 | 480 500 500 
39 Forward Voltage <0.2 | 350 180 180 
56 Self Bias-Volts <1 70 150 110 
56 Reverse Voltage 0.5 | 440 560 560 
56 Forward Voltage 0.3 | 320 320 240 
66 Self Bias-Volts 100 85 
66 Reverse Voltage 560 520 
66 Forward Voltage 200 320 


MECHANICAL EFFECTS OF BOMBARDMENT 


The marked changes in the electrical properties of silicon imposed by 
bombardment strongly suggest that bombardment may also impose a 
corresponding change in the lattice structure and that this might be 
detected by suitable optical methods. Attempts were made at an early 
date to detect such changes. To this end a mask of nichrome ribbon 5 
mils wide and 1 mil thick was laid over a sample of silicon during bom- 
bardment. An optical examination of the surface showed that after 
bombardment in the case of helium the surface on either side of the 
mask was elevated whereas in the case of argon it was depressed. This 
result has since been confirmed by one of the authors’s colleagues, 
Dr. F. W. Reynolds, who has found that in cases of prolonged bom- 
bardment by helium the adjacent surfaces may be elevated by as much 
as 225 Angstroms’ while in the cases of prolonged bombardment by 


7One Angstrom is 1078 cm. 
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argon the adjacent surfaces may be depressed by as much as 130 Ang- 
stroms. Further investigation of this phenomenon is under way. 


STABILITY OF BOMBARDED SURFACES 


No extended test has yet been made of the stability of bombarded 
surfaces but results extending over more than two years are encouraging. 
Similarly, rectifiers for the,millimeter wavelength range, mounted with- 
out the usual protective impregnation, show little or no change at the 
end of a year. _ 

In a few instances bombarded surfaces have been subjected to rather 
severe tests with results that suggest that under normal conditions they 
may be even more stable than surfaces activated by more conventional 
methods. For example, surfaces contaminated while cutting or while 
cementing them to their mountings have subsequently been cleaned 
with solvents such as alcohol and are substantially the same before and 
after treatment. In other cases, they have been heated in a flame to 
soldering temperatures with no appreciable effects. Even in the very 
severe case where the bombarded piece was heated to a cherry red and 
the superficially oxidized layer was removed with hydrofluoric acid the 
effects of bombardment were still evident. There was, however, con- 
siderable reduction in the tolerable reverse voltage. There is nothing 
in our experience to date to suggest that bombarded surfaces treated 
in accordance with the simple straightforward methods outlined above, 
are in any wise temporary in character. 


CONCLUSIONS 


The experiments reported above have shown that rather pronounced 
changes in the electrical properties of silicon may be produced merely 
by bombarding the polished surface with positive ions. The ratio of 
forward to reverse currents, for example, which for the usual untreated 
silicon is seldom more than a few hundred, can be made more than 10,000. 
Ixperiments show that the effect depends to some extent on the type 
of ion gas used, helium being a preferred medium. The effect depends 
also on the volocity of the bombarding particles, the total bombarding 
charge and to a lesser extent on the temperature of the specimen during 
bombardment. Good results are obtained from bombarding potentials 
of 80 kv with current densities of 5 microamperes per square centimeter 
for periods of one or two minutes. The temperature should preferable 
be about 300°C. 

Ordinarily the properties of silicon are materially affected by impurity 
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content. In the case of bombarded silicon the effect is much less. More 
particularly it-is possible to contaminate silicon with impurities such 
as boron to the point where its rectifying properties are almost com- 
pletely lost and by bombardment it is possible to convert the crystal 
into a very useful rectifier. It is possible to produce results similar to 
the above by exposing the crystal to radioactive polonium. Bombarded 
materials appear to be relatively stable. 


The writer wishes to express his appreciation of the encouragement 
and help of Dr. G. C. Southworth in the preparation of this paper, to 
A. J. Mohr, Jr., for his able assistance in the experimental work, and 
to numerous associates in Bell Telephone Laboratories for their assist- 
ance in preparing materials and in making special tests for which the 
author was not adequately equipped. 


Mechanical Properties of Polymers 
at Ultrasonic Frequencies 


BY WARREN. P. MASON AND H. J. McSKIMIN 
(Manuscript received October 25, 1951) 


Since the mechanical properties of solid polymer materials are largely de- 
pendent on the motions that segments of the polymer chains can undergo, 
to understand these properties one must use measuring techniques which can 
determine these motions. One.of the most promising methods is to measure 
the reaction of polymer materials to longitudinal and shear waves over a 
frequency spectrum wide enough to determine the relaxation frequencies due 
to thermal motions of the principle elements of the chain. The presence of 
relaxations is indicated by a dispersion in the velocity and attenuation 
constants of the material, or a dispersion in the characteristic impedance 
of the material of the attenuation ts too high to allow velocity measurements. 
A number of different types of measuring methods are described in this 
paper which make possible propagation and impedance measurements not 
only in solid polymers, but also in liquid polymers and in solutions of 
polymer molecules in typical solvents. 

When these techniques are applied to long chain polymers in dilute solu- 
tions, the three relaxations observed correspond to motions occurring in isolated 
molecules since as the dilution increases, the molecules seldom touch. The 
lowest relaxation corresponds to a configurational relaxation of the molecule 
as a whole, the highest relaxation corresponds to the twisting of the shortest 
segment—containing about 40 repeating units—while the intermediate re- 
laxation corresponds to a transient entanglement of chain segments. All 
three types of relaxations are present in pure polymer liquids but are spread 
out over a frequency range due to the perturbing effect of near neighbors of 
adjacent chains. The high frequency shortest chain relaxation can be traced 
in solid polymers of the linear chain type such as polyethylene and nylon 
and produces rubber-like response to mechanical shocks of very short duration. 


I. INTRODUCTION 


The mechanical properties of solid polymer materials are largely deter- 
mined by what motions, parts or segments of the polymer chains can 
undergo. Toughness, mechanical impact strength and ultimate elongation 
depend on the facility with which the polymer molecule can be displaced. 
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If only a small motion of the polymer chain can occur within the time of 
the measurement, the material has high elastic stiffness coefficients and 
acts similar to a rigid solid. On the other hand, if significant segments 
of the polymer chain can move at the frequency of measurement, the 
elastic stiffness 1s much lower and rubber-like behavior results. An inter- 
mediate case, which occurs when the significant motion of the polymer 
molecule is near the relaxation time at the frequency of measurement, is 
that of a damping material such as butyl rubber. Even “long time” 
qualities of plastics such as creep, stress relaxation and recovery depend 
on the integrated displacements of rapidly oscillating segments of the 
chain. 

One of the most promising methods for investigating these motions is 
to determine the reaction of mechanical waves on the polymer materials 
over a wide spectrum of wavelengths, eventually going to frequencies 
comparable with those of thermal vibrations of significant groups or 
segments in the macromolecules. 

If one wishes to understand the origins of these motions it is necessary 
to measure the molecules in the form of liquids or solutions since then 
the segments of the molecule are less restrained by their neighbors and 
can perform all the possible vibrations. Polymer liquids are also interest- 
ing in themselves as sources of damping material. To apply these results 
to rubbers and solid materials, one then has to measure the modifications 
of the polymer chain motion caused by the close approach of near neigh- 
bors, by measuring the mechanical properties of these materials. 

By using different types of techniques, these processes can be applied 
to molecules in solution, to liquid polymers and to solid polymers. The 
principal types of methods used for liquids are the torsional crystal, the 
torsional wave propagation system and the shear wave reflectance method, 
all of which are described in Section II. For solids an optical method and 
an ultrasonic method are described in Section V. All of these methods 
involve displacements of 10~-® cm or less so that non-linear effects are 
negligible. 

All of these methods depend on setting up shear or longitudinal waves 
in the medium and observing either the velocity and attenuation of the 
wave, or the reaction of the medium back on the properties of the 
transducer. If the attenuation of a wave in the medium under consid- 
eration is low enough to permit the wave parameters, i.e., the velocity 
and attenuation per wavelength to be determined, the relaxation of 
some significant part of the polymer molecule is determined by the dis- 
persion of the wave properties which occur, as shown by Fig. 1A, in the 
form of an increase in velocity and a maximum in the attenuation per 
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SQUARE OF VELOCITY => 
ATTENUATION IN NEPERS PER WAVELENGTH => 





LOG jg OF SOUND FREQUENCY, f 
Fig. 1A—Velocity and attenuation for a medium with one relaxation frequency. 


wavelength curve. If the variation in this relaxation mechanism is 
studied as a function of temperature and chain length, the type of seg- 
ment may be determined. If, however, the attenuation of the medium 
is so high that its wave properties cannot be determined, some informa- 
tion can still be obtained by determining the loading, or mechanical 
impedance, that such a wave exerts on the driving crystal or trans- 
ducer. If all the relaxations occur in the stress-strain relation, it can be 
shown that there is a reciprocal relation between the propagation con- 
stant [ = A + jB, and the characteristic impedance per square centi- 
meter Zo given by the equation 


ZV = (R + jX)(A + JB) = jop (1) 


where A is the attenuation and B the phase shift per centimeter, & the 
mechanical resistance and X the mechanical reactance per square centi- 
meter, w is 27 times the frequency and p the density of the medium. A 
typical two relaxation mechanism® is shown by the curves of Fig. 1B. 
By assuming values for the stiffness and dissipation factors and fitting 
a theoretical curve to the measured values, the relaxation frequency or 
frequencies can be determined. 

1 All the relaxation mechanisms discussed in this paper are represented in 
terms of equivalent parallel electric circuits in which the resistance terms repre- 
sent viscosities and the inverse of capacities represent shear elastic stiffnesses. 


In mechanical terms these correspond to a series of Maxwell models as discussed 
in a paper by Baker and Heiss to be published in the next issue. 
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The most information about the motions of isolated polymer chains 
can be obtained by investigating the properties of polymer solutions. 
This follows from the fact that in pure polymer liquids, and in solids, 
the mechanical properties are mainly determined by interactions be- 
tween chains on account of the close packing of the chains. If, however, 
one dissolves the polymer molecules in a solvent, the inter-chain and 
intra-chain reactions can be separated as the dilution increases. When 
the polymer is in the order of one percent of the solvent, the chains on 
the average touch very seldom and the mechanical properties of the 
solution are determined by the properties of single molecules. As dis- 
cussed in Section JII, three types of chain segment motion have been 
isolated, (1) a configurational relaxation of the chain as a whole, (2) a 
position change of the shortest segment and (8) twisting of the shortest 
chain segment. Above the frequency of relaxation of this chain segment 
the joints of the polymer molecule become frozen and the chain be- 
comes very stiff. These shortest chain relaxations occur also in pure 
polymer liquids, in rubbers and in non rigid solids with linear chain 
segments such as polyethylene. In pure liquids a lower frequency quasi- 
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Fig. 1B—Mechanical impedance loading for a medium with two relaxation 
frequencies. 
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configurational relaxation also occurs for chain lengths greater than 60 
elements but for chain lengths less than 40 elements this type of relaxa- 
tion disappears. From the difference between the high frequency shear 
elasticities measured for polyethylene and nylon 6-6 and those measured 
for static pulls, it appears that there may be lower frequency relaxations 
in these materials as well. 


II. METHODS OF MEASUREMENT FOR SOLUTIONS AND PURE POLYMER 
LIQUIDS 


To measure the mechanical properties of such dilute solutions, shear 
waves have been used since for longitudinal waves the added stiffness 
caused by the dissolved polymers is very small compared to the stiffness 
of the solvent alone. The velocity and attenuation of a longitudinal 
wave are given by the equations 


y= te, A= 2S ty 4 On (2) 


p pv? 





where \ and yp are the Lamé elastic constants, f the frequency, p the 
density, v the sound velocity, x the compressional viscosity and 7 the 
shear viscosity. Since for a one percent solution of polyisobutylene in 
cyclohexane the shear elasticity does not exceed 90,000 dynes/cm’, 
whereas the value of ) is in the order of 2 X 10" dynes/cm’, it is obvious 
that the longitudinal velocity would have to be measured to an accuracy 
of 1 part in 100,000 before the presence of polymer molecules could be 
ascertained. Attenuation measurements give some information on the 
added viscosity due to the chain molecules but since longitudinal attenu- 
ations are not easily measured below 1 megacycle, the most interesting 
frequency range is missed. 

A pure shear wave in a viscous liquid is propagated according to the 
equation? 


wip ; 
v = ne~V @ bt (3) 


where »v is the transverse particle velocity, p the density, f the frequency, 
n the shear viscosity, 7 = +/—1 and z the distance. For typical sol- 
vents, the attenuation is so high that wave motion cannot be measured. 
However the viscous wave produces an impedance loading on a crystal 
generating such a wave which can be measured by the change in the 
resonant frequency and the change in the resistance at resonance. The 
mechanical impedance per square centimeter caused by such a viscous 
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wave is equal to” 
Zo = Vafno lt 9) = Ru +jXm (4) 


This causes a change in resistance, and a change in frequency in a 
crystal generating a shear wave in the liquid equal to 


AR = Kika; Af = —K2X M (5) 


where K, and Ky are constants of the crystal which can be obtained 
approximately from the dimensions and piezoelectric constants of the 
_ crystal but which are more accurately obtained by calibration in known 
liquids. The constants K, and Ke vary slightly with temperature and 
should be calibrated over a temperature range. 

The first instrument to use a vibrational method for measuring vis- 
cosity was the vibrating wire method of Phillipoff.* In this method a 
wire was vibrated in a liquid and the damping rate was used as a measure 
of the viscosity. Another method also applicable in the low frequency 
range is the transducer method of Ferry.’ In this method wires are 
vibrated by electromagnetic transducers and the resistance and reac- 
tance drag on the wires are measured by the change in the electrical 
resistance and reactance of the transducer. From the constants of the 
transducer, the equivalent viscosity and stiffness of the liquid can be 
measured. 

In the medium frequency range a torsional crystal’? method was 
devised by one of the writers which has been applied in the frequency 
range from 10 to 150 ke. The torsional erystal is shown by Fig. 2. For. 
these types of measurements the crystal usually is made of quartz with 
four electrodes of gold evaporated on the surface. Four wires are soldered 
on the surface and serve as supports as well as electrodes. The motion 
is all tangential to the surface and tests at Bell Laboratories and at the 
Franklin Institute,’ where a precision study of the torsional crystal has 
been made, have shown no observable longitudinal waves from the 
crystal surface. The process of measurement consists in measuring the 

2W. P. Mason, Piezoelectric Crystals and Their Application to Ultrasonics, 
D. Van Nostrand, 1950, p. 340. 

3 W. Phillipoff, Physik. Zeits, 35, 1934, pp. 884-900. 

4T.L. Smith, J. D. Ferry and F. W. Schemp, ‘‘Measurement of the Mechanical 
Properties of Polymer Solutions by Electromagnetic Transducers,’’ J. App. 
Phys., 20, No. 2, Feb. 1949, pp. 144-153. 

5 W. P. Mason, ‘‘Measurement of the Viscosity and Shear Elasticity of Liquids 
by Means of a Torsionally Vibrating Crystal,’’ A.S.M.E., 69, May 1947, pp. 359- 
ae P. E. Rouse, Jr., E. D. Bailey, and J. A. Minkin, ‘Factors Affecting the 
Precision of Viscosity Measurements with the Torsional Crystal,’’ Laboratories 


of the Franklin Inst., Report 2048, presented to Am. Petroleum Inst., May 4, 
1950. 


128 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1952 


resonant frequency and resonant resistance of the crystal in a vacuum, 
then introducing the solution to be measured, the change in the resonant 
resistance Al?, and the change in resonant frequency Af are determined 
by an electrical bridge. Several short cuts are possible if the mechanical 
impedance is not too high. By measuring the capacity at a frequency 
considerably higher than the crystal frequency, the resistance at reso- 
nance and Af can be obtained by changing the frequency and resistance 
until a balance is obtained leaving the capacity unchanged. This method 
has been used to measure viscosity, and a recent precision study at the 
Franklin Institute’ has shown that it agrees with other methods to an 
accuracy of well under a per cent. 

The torsional quartz crystal has been successfully used to measure 
liquids having a viscosity up to 10 poise, but above this viscosity the 
electrical resistance gets so high that it is hard to measure it since it is 
shunted by the much smaller reactance of the static capacitance of the 
crystal. A crystal of higher electromechanical coupling such as am- 





Fig. 2—Cell and 80-ke crystal for shear viscosity and elasticity measurements 
of liquids. 
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Fig. 3—Photograph of torsional crystal and rod. 


monium dihydrogen phosphate (ADP) will cause the electrical resistance 
component to be smaller in comparison to the reactance of the static 
capacitance and hence can be used to measure higher viscosities. How- 
ever, since wires cannot be soldered to the surface but must be glued, 
the crystal is much more fragile than quartz and its use has been aban- 
doned in favor of another method which makes use of the phase and 
attenuation change in a torsional wave in a rod caused by the surround- 
ing liquid whose properties are to be measured. 

This method, devised by one of the writers,’ consists in sending a short 
train of torsional waves, periodically repeated, down a glass or metal 
rod. As shown by Fig. 3, the torsional wave is generated by a torsional 
quartz crystal soldered or glued to the end of the rod. These waves 
travel to the free end of the rod and are reflected back to the crystal 
where they are detected, amplified, and displayed on a cathode ray 


7 This method is described by H. J. McSkimin, in a paper before the Ac- 
oustical Soc. of Am. in October, 1951. 
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oscilloscope. Echoes due to end to end reflections also appear, being 
attenuated by normal acoustic losses until they are undetectable by the 
time the next pulse is applied. 

With only air surrounding the rod, a phase reference and amplitude 
reference are obtained for the first received wave (or subsequent echoes 
if greater sensitivity is desired). The rod is then immersed a definite 
length in the liquid to be measured, as shown in Fig. 4, with a resulting 
phase retardation and amplitude reduction. These are measured by em- 
ploying the experimental circuit shown by Fig. 5. In order to allow the 





Fig. 4—Photograph of complete torsional wave measuring equipment. 
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use of one crystal for both recelving and transmitting, the crystal is put 
in the bridge circuit of Fig. 5 where a resistance and capacity are used 
to balance out the transmitted pulse so that it will not overload the 
amplifier. The relatively weak voltages generated by the incoming 
acoustic waves pass through directly. The gate circuit provides pulses 
of radio frequency voltage at repetition rates in the range of 20 to 100 
per second with a synchronizing voltage supplied to the oscilloscope for 
the horizontal sweep. The frequency range of the device is from 20 to 
200 ke. Both glass and nickel-iron rods were used, the latter having a 
very low frequency-temperature coefficient. With a 100-kce quartz 
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Fig. 5—LExperimental pulsing circuit for measuring torsional impedance of 
liquids. 
torsional crystal, a rod length of 21 inches and diameter of 0.2 inch were 
used. The entire crystal-rod assembly is placed inside a glass tempera- 
ture control unit, as shown by Fig. 4, through which water can be cir- 
culated to provide temperatures in the range 0°C to 80°C. The test 
liquid is placed either directly into the inner bore of this water jacket, 
or in another tube which can be inserted from the bottom to surround 
the rod up to a fixed mark. 

In use, both phase and attenuator settings were adjusted to balance 
the first received pulse against the continuous wave component passing 
through the attenuator. Cancellation for the duration of the pulse was 
visually indicated on the oscilloscope. A plot of balance phase and level 
is made as a function of the temperature. When the liquid is introduced 
an attenuation change AA and a phase change AB are required to 


132 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1952 


re-balance the circuit. These are measured by the amount of attenuation 
in nepers (1 neper = 8.68 db) and the number of radians phase shift 
required to re-establish a balance. An alternate method of measuring 
phase shift is to measure the change in frequency required to re-estab- 
lish balance. If this method is used the phase shift change of the overall 
circuit with frequency has to be calibrated for the uncovered rod by 


E 


Fig. 6—High frequency shear reflection method for measuring shear imped- 
ances of liquids. 
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noting the frequencies for which 360° phase shifts (as measured by 
balance) occur in the circuit. 

It is shown in the appendix that the torsional impedance of the 
liquid per square centimeter is given by 


Z = (@)§ [AA + JAB] (6) 


where p is the density, and vp) the sound velocity in the rod, a is the 
radius and ¢ the covered length of the rod and AA and AB are respec- 
tively the change in attenuation in nepers and the change in phase shift 
in radians to re-establish balance. If a very viscous liquid is used it may 
be necessary to correct for the fact that the torsional impedance may 
differ from the plane wave impedance as discussed in the appendix. 

This device can measure liquids having dynamic viscosities from 10 
poise to 1,000 poise with an accuracy of the order of 10 per cent. The 
frequency range covered may be from 20 to 200 ke depending on the 
size of the crystal used to drive the rod. Hence it supplements the 
torsional crystal method for very viscous liquids. 

At frequencies above 500 ke, the torsional crystal becomes too small 
to be used practically and recourse is had to a high frequency pulsing 
method.’ As shown by Fig. 6, shear waves are set up in a fused quartz 

8 W. P. Mason, W. O. Baker, H. J. McSkimin and J. H. Heiss, ‘‘“Measurements 
of the Shear Elasticity and Viscosity of Liquids by Means of Ultrasonic Shear 
Waves,”’ Phys. Rev., 75, No. 6, March 15, 1949, pp. 936-946. See also H. T. O’Neil, 


‘Refraction and Reflection of Plane Shear Waves in Viscoelastic Media,’’ Phys. 
Rev., 75, No. 6, March 15, 1949, pp. 928-936. 
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rod by means of a Y-cut or AT cut crystals soldered to a silver paste 
layer baked on the fused quartz surface. The particle motion of the 
shear wave is parallel to the large reflecting surface and hence only 
shear waves are reflected from this surface. These impinge on a second 
shear crystal which is connected to an amplifier and oscillograph. Since 
the attenuation in fused quartz is so low, a long series of reflected pulses 
appear on the oscillograph. When a liquid, whose shear properties are 
to be measured, is placed on the fused quartz surface, this causes a 
change in the amplitude and phase of the reflected wave. By using the 
balance method shown by Fig. 7, in which two identical fused quartz 
rods are used, one of which has a liquid layer and the other does not, 
and by using a phase shifting network and an attenuator to balance out 
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Fig. 7—Method for obtaining resistance and reactance terms for high fre- 
quency shear reflection method. 


pulses, the shear impedance of the liquid can be determined. If R is the 
loss per reflection expressed as a current ratio, @ the change in phase 
angle required to rebalance the circuit and ¢ the angle between the wave 
normal and the reflecting surface, it can be shown’ that the shear im- 


pedance of the liquid is 
1 — R’ + QR sin | (7) 
1+ fF? + 2R cos 6 


where Z, is the impedance pv for shear waves in the quartz. This is 
equal to 


Zq = 2.20 X 3.76 X 10° = 8.27 < 10° mechanical ohms 


Zu = Ru + jXu = Zq COS | 


Since this impedance is much larger than that of the liquids that are to 
be measured, the sensitivity is increased by making ¢ large. In practice 
g was taken as 80°. This method is applicable from 3 mc up to 100 mc 
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and complements the other methods. Fig. 8 shows a photograph of the 
equipment. 


III. MEASUREMENTS OF POLYMERS IN SOLUTION 


When such methods are applied to a polymer solution, it is found 
that the resistance and reactance components are no longer equal but 
the resistance is invariably larger than the reactance. This indicates the 
presence of a shear elasticity in the solution. If the molecules have a 
single relaxation frequency, it has been found that the shear properties 
of the liquid can be represented by a stress-strain equation of the type 

a ee 
i oe (3) 

os bps 

NB OL 
where 7’ is the shearing stress, S the shearing strain, 7, the solvent 
viscosity, nz a molecular viscosity of some particular motion of the 
chain which disappears when the reactance of the chain stiffness uz of 
this motion is low enough so that the motion can follow the applied 
shearing stress at the frequency of the measurement. When this type of 
mechanism is present in the liquid, it has been shown’ that the impedance 
the liquid presents to the crystals is 


4 
punna + j | eon + PEP na + m) | 
w (9) 
nz + pp/o 
Fig. 9 shows a plot of the resistance and reactance components of an 
assumed solution having a single relaxation frequency, and a viscosity 
30 times the solvent viscosity. At very low frequencies, the resistance 
and reactance follow that of a solution, but for frequencies comparable 
with the relaxation frequency, the resistance becomes larger than the 
reactance while for very high frequencies the two come together on a 
line determined by the solvent viscosity. If there is more than one 
relaxation frequency, the resistance and reactance may coalesce for 
several intermediate stages. A continuous distribution would give a 
definite relation between the frequency dependence of resistance and 
reactance. . 
The torsional crystal and the shear wave reflection method have 
been applied to long chains of polyisobutylene dissolved in various sol- 


Zo = Ru + jXu = 





9W. P. Mason, Piezoelectric Crystals and Their Application to Ultrasonics, 
D. van Nostrand, 1950, p. 353. 
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Fig. 8—Dual reflecting block ‘assembly for measuring high frequency shear 
impedances of liquids. 


vents with concentrations ranging from zero per cent to 10 per cent. 
Polyisobutylene is a polymer molecule having the chemical formula 


i 
| 
CH; H 


Non-planar zigzag segments can be expected in the liquid state. Fig. 10 
shows measured curves for 20 ke of the resistance and reactance for 
solutions of viscosity average molecular weight of 3,930,000 in cyclo- 
hexane. Four values of concentration were used and two temperatures 
were measured. For pure cyclohexane, the resistance and reactance 
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Fig. 9—Resistance and reactance components of a solution having a single 
relaxation frequency and a solution viscosity 30 times the solvent viscosity. 
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Fig. 10—Shear resistance and reactance components of a solution of polyiso- 
butylene (molecular weight of 3,930,000) in cyclohexane plotted as a function of 
grams per cc of solution. 
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components are equal but as the percentage of polyisobutylene is in- 
creased, the resistance increases more rapidly than the reactance. 

By solving equation (9) for 74 , nz and wz in terms of R and X measured 
at one frequency and 74 + 72 the solution viscosity, we find 
| _2RX — (R* — X”)*/wp 


als wp wplna + 18) — ORX ’ ms = (na + te) — 14 


(R? — x’) ONB 
wp(na + nb) — 2RX 


Applying these formulae to the measured results, the curves of Fig. 11 
result. The shear elasticity is directly proportional to the concentration, 
the viscosity 74 is only slightly larger than the solvent viscosity while 
the main part of the measured viscosity resides in 7, the viscosity asso- 
ciated with chain motion. Fig. 12 shows these three quantities for a one 
per cent solution measured as a function of temperature. The apparent 
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‘Fig. 13—Shear elasticity for a one per cent solution of polyisobutylene as a 
function of frequency for 25°C. 
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stiffness decreases with increase in temperature for a single frequency 
measurement. However, when measurements were made at 20, 40, 80 
and 150 ke it was found that the elasticity was a function of the fre- 
quency, which indicates the presence of more than one relaxation and 
complicates the determination of the temperature relationships. Fig. 13 
shows measurements of the shear elasticity over a frequency range from 
2.5 ke! up to 14 megacycles for 25°C. There is a gradual rise up to 
about 300 kilocycles after which there is a sharp break to a stiffness of 
about 90,000 dynes/cm’ for a 1 per cent solution of the 3,930,000 molecu- 
lar weight polymer in cyclohexane. If one analyzes the frequency varia- 
tion of the elasticity he finds that it can be fitted by three relaxation 
frequencies, one having a frequency of 230 cycles, one around 66,000 
cycles and one around 4 megacycles. A possibility exists for a fourth 
relaxation. The lowest relaxation is thought to be a configurational re- 
laxation of all the elements of the chain. The highest one appears to be 
a relaxation of the twisting motion of the smallest segment of the chain 
such as a Kuhn segment. The intermediate relaxation appears to be due 
to the motion of the ends of the smallest chain segment from one posi- 
tion of entanglement to an adjacent position. This interpretation is 
based partly on the fact that the associated viscosity of the motion is 
very similar to that for the relaxation of the twisting motion of the 
smallest chain segment and partly from data presented in the next 
section on pure polymer liquids which shows a lower frequency relaxa- 
tion agreeing in frequency assymptotically with this one, which involves 
chain motions of approximately 30 to 40 chain elements. Temperature 
-variations of these elastic components show that the lowest relaxation 
mechanism has a stiffness that increases slightly with temperature in 
agreement with the kinetic theory of elasticity. The corresponding 
viscosity (2) which comprises most of the viscosity for a solution, when 
plotted against the reciprocal of the temperature, as shown by Fig. 14, 
indicates an activation energy. of 3.9 kilocalories per mole which is 
slightly higher than that of the solvent cyclohexane alone, which is 
about 3.2 kilocalories per mole. This difference of 0.7 kilocalories pre- 
sumably represents the added energy required to bend the chain in its 
configurational motion. Measurements with another chain length of 
1.18 X 10° molecular weight showed that the stiffness of the lowest (con- 
figurational) relaxation decreased from 310 to 160 indicating that the 
stiffness of this motion is approximately proportional to the square root 

10 The lowest frequency, 2.5 ke, was measured by means of a quartz crystal 


tuning fork which will be described in another paper. This instrument makes 
possible the direct measurement of configurational elasticities. 
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of the molecular weight. The viscosity decreased by a factor of 6.25 
and consequently the relaxation frequency increases from 230 to 660 
cycles. 

The second “entanglement” relaxation has a stiffness of about 3100 
dynes/cm’ for a 1 gram per cc solution of the 3.93 X 10° molecular 
weight solution and about 2650 dynes/cm’ for the 1.18 X 10° molecular 
weight solution. The variation with temperature, if any, is small. The 
corresponding viscosities 73; for the two solutions are nearly equal as 
shown by Fig. 14, and have an activation energy of 4.25 kilocalories 
per mole. The final high frequency “short segment” relaxation has a 
high stiffness of 83,000 dynes/cm’ for a 1 per cent solution of 3.93 x 10° 
molecular weight. The corresponding viscosities for the two solutions 
shown by Fig. 14 have nearly identical values and an activation energy 
of 4.25 kilocalories per mole, i.e. very closely equal to the “entangle- 
ment” relaxation viscosity. 

These upper two relaxations persist in pure liquid polymers as dis- 
cussed in the next section, although they are spread out over a small 
range of relaxation time values. The highest one can be traced in meas- 
urements of mechanical properties of solid plastics such as polyethylene 
and nylon which indicates that these materials should have rubber like 
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Fig. 14—Viscosities for the three components of the motion and for the solvent 
plotted against the inverse of the absolute temperature. 
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. Density Melt Viscosity Poise 
Average = : 
Liquid | Molecular oe 
Weight a atc 
15°C | 25°C | 35°C | 50°C | 15°C | 25°C | 35°C | 50°C 
A 904 | 0.83 | 0.826 | 0.821 | 0.814 | 1.10 | 0.67 | 0.33 | 0.16 16.2 
B 1,220 | 0.846 | 0.842 | 0.837 | 0.83 | 2.64 | 1.40 | 0.65 | 0.297 22.2 
C 1,590 | 0.856 | 0.849 | 0.846 | 0.840 | 8.4 | 4.0 | 2.1 | 0.87 28.4 
D 2,450 0.872 | 0.866 | 0.863 | 0.857 100 38 | 15.9 4.9 44.0 
E 3,520 0.881 | 0.877 | 0.872 | 0.866 580 216 95 30 63.0 
F 4,550 | 0.89 | 0.886 | 0.882 | 0.875 | 1,960 | 737] 320] 98 81.4 
G 5,590 0.897 | 0.892 | 0.887 | 0.882 | 4,600 | 1,840 740 | 220 100 
H 10,380 0.912 | 0.908 | 0.904 | 0.896 | 11,500 | 4,600 | 2,050 | 625 186 
uw dynes j . 
Measurement at 24 megacycles mie ee 
15°C 25°C 35°C 50°C 
Liquid 
B 7 H ” B ” B ” 
A 0.482 X 109 | 1.07 | 0.4 109 | 0.68 | 0.35 xX 109 | 0.39 | 0.1 x 109 | 0.13 
B 1.08 X 109 | 2.64 | 0.8 X 109 | 1.39 | 0.58 X109 | 0.7 | 0.17 X 109 | 0.293 
Cc 1.65 109 | 5.12 | 1.2 109 | 2.57 | 0.885 X 109 | 1.99 | 0.27 X 109 | 0.881 
D 2.71 X 109 | 15.7 | 2.0 X 109] 8.1 | 1.67 109 | 4.15] 0.98 x 109 | 1.77 
E 3.81 X 109 | 40 3.0 X 109 | 20 2.17 X 109 | 10.6 | 1.2 xX 109 | 4.62 
F 5.48 xX 109 | 73 4.22 X 109 | 38.3 | 2.8 XX 109 | 20.2 | 1.6 109 | 8.8 
G 5.9 xX 109 | 107 4.78 X 109 | 51 3.74 X 109 | 29.5 | 2.5 X 109 | 13.2 
H 6.9 X 109 | 119 5.55 X 109 | 65.8 | 4.22 xX 109 | 31.5 | 2.9 X 109 | 14.8 
Measurements at 14 megacycles 
A 0.38 X 10° | 0.975 | 0.8 X 109 | 0.6 | 0.22 X 109 | 0.35] 0.12 X 109 | 0.18 
B 0.475 X 109 | 2.53 | 0.35 X 109 | 1.4 | 0.25 x 109 | 0.74] 0.18 X 10° | 0.73 
Cc 0.68 X109 | 3.86 | 0.61109 | 3.4 | 0.8 x 109 | 1.62] 0.25 xX 109 | 0.56 
D 2.32 X109 | 16.3 | 1.7 xX 109 | 10 0.94 X 109 | 4.75] 0.39 X 109 | 1.86 
E 3.8 X 109 56.2 2.8 X 109 | 24.25 | 2.0 X 109 | 12.3 0.855 X 10° 5.7 
F 4.75 X109 | 90.4 | 3.6 X 109 | 48.3 | 2.65 xX 109 | 26.6 | 1.47 X 109 | 10.8 
G 6.03 X 109 | 136.5 4.6 X 109 | 81.4 3.24 X 109 | 39.6 2.0 % 109 | 15 
H 6.64 XX 109 | 160 5.3 X 109 | 93.5 3.98 X 109 | 54.2 2.3 109 | 18.5 
Measurement at 4.5 megacycles 
A 0.34 X 109] 1.18 | 0.19 X 109 | 0.64] 0.18% 109 | 0.34] 0.09 x 109 | 0.17 
B 0.55 X109| 2 0.21 X 109 | 1.33] 0.2 x 109 | 0.68 | 0.1 XX 109} 0.20 
Cc 
D 1.57 X 109 | 24 0.96 X 109 | 11.7 | 0.72 X 109 | 5.52 | 0.34 X 109 | 2.1 
E 2.79 X 109 | 76 1.9 X 109 | 37.4 | 1.43% 109 | 16.5 | 0.9 X109 | 5.15 
F 3.57 X 109 | 124 2.5 X 109 | 61.2 | 1.86 X 109 | 31 1.04 X 109 | 12.8 
G 4.4 X 109 | 176 3.4 X 109 | 98.5 | 2.6 x 109 | 54 1.6 X 109 | 22.6 
H 5.65 X 109 | 186 4.3 X 109 |124 2.8 X 109 | 76 1.8 XX 109 | 28.8 
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response to mechanical shocks of very short duration. The lowest fre- 
quency configurational relaxation is spread over a wide spectrum of 
relaxation times in pure liquids. 

Measurements” of these and other chains in various solvents have 
also been made and the results are discussed, from a chemical point of 
view, in a companion paper by W. O. Baker and J. H. Heiss. It is shown 
that the stiffnesses vary with the polymer chain and the solvent used. 


IV. MEASUREMENTS OF PURE LIQUID POLYMERS 


A. Shear Wave Measurements in Liquid Polymers 


Similar shear wave measurements have been made for pure polyiso- 
butylene liquids of molecular weights from 904 to 10,380 (i.e. from 16 
chain elements to 186 chain elements), by the techniques described in 
Section II. Some of these results have been discussed in reference (8) 
but the much more comprehensive measurements made since require 
some revisions of the original conclusions. 

The easiest data to interpret are the high-frequency data obtained by 
the shear wave reflectance method. The data of Table I give measure- 
ments of 8 liquids varying in average molecular weight from 900 to 
10,380, at three frequencies and four temperatures. If we plot for ex- 
ample the Maxwell shear stiffness and viscosity for the three frequencies 
and for 25°C as a function of the number of chain elements (here a chain 
element is taken as two adjacent carbon atoms one of which has two 
methyl groups attached and the other two hydrogens) the 4.5-mc meas- 
urements are shown by the triangles of Fig. 15. The 14 megacycle 
measurements are shown by the circles and the 24-mc measurements by 
the squares. 

An attempt was made to fit these measurements with a two relaxa- 
tion mechanism shown by the figure with two stiffnesses which are taken 
to be independent of the molecular weight and equal respectively to 
1.2 X 108 dynes/cm’ and 6 X 10° dynes/cm’. The best fit is obtained by 
taking the two viscosities 7; and 72 equal and these are adjusted for the 
different molecular weights in such a manner as to best fit the experi- 
mental curve. A fair agreement is obtained except for the range from 60 
to 90 chain elements where the two relaxation model gives too rapid an 
increase of stiffness with increase in the number of chain elements and 
at the high molecular weight viscosity range where the viscosity shows 
a dispersion in values but the model does not. The sum of the two vis- 

1 These results on the mechanical impedance of long chain molecules in sol- 


vents’ have been presented at the XIIth International Congress of Pure and 
Applied Chemistry by W. O. Baker, W. P. Mason and J. H. Heiss, Sept. 13, 1951. 
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cosities m and m2 assumed as a function of molecular weight is shown: by 
Fig. 16. The log of the viscosity starts proportional to the molecular 
weight but above a molecular weight of 2,400 the increase is very slow 
and becomes asymptotic to a value of 240 poises. An equation which 
fits the increase in viscosity with molecular weight is 


11-8 tanh Z/2370 
nD = Ke anh Z/ 


where Z is the molecular weight. The solid line shows a plot of this curve 
and the circles are the assumed values to obtain a best fit to the meas- 
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Fig. 15—Measured values of high frequency shear viscosity and elasticity for 
25°C and three frequencies plotted against molecular weight. Solid lines are best 
fit obtained by a two relaxation mechanism having the element values shown by 
the figure. 
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ured values. This equation indicates that when Z = 2,370 or 48 chain 
elements the viscosity increases only a small amount more by a chain 
articulation effect and hence in this high frequency range we are dealing 
with a chain length of about 40 elements or 80 carbon atoms. This is 
checked also by a comparison of the static and dynamic viscosity. The 
total dynamic viscosity due to the two relaxation mechanisms compared 
to the static viscosity does not differ markedly until the number of 
chain elements is more than- 40. Above this value other motions than 
that of the shortest chain segment can take place and can add to the 
dynamic viscosity. The static viscosity fits an equation of the same sort, 
but the indicated chain length for the viscous motion is about $ times 
that of the shortest segment. 

When a similar process is carried out over the temperature range the 
equations of Fig. 16 are obtained. The static viscosity has an activation 
energy of 16 kilocalories per mole, while the dynamic viscosity has an 
activation energy of about 11.2 kilocalories per mole. 

The relaxation frequencies for the two components are plotted as a 
function of the number of chain segments by the solid lines of Fig. 17. 
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Fig. 16—Triangles are measured static viscosities and circles are dynamic 
viscosities plotted as a function of molecular weight. Solid lines are a plot of the 
equations given for static and dynamic viscosities. 
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For long chain segments the values become asymptotic to 8 X 10° cycles 
and 160,000 cycles which are not far from the two highest relaxation 
frequencies obtained from the solution measurements of Section ITI. 
Hence it appears likely that these relaxations are due to the “entangle- 
ment”? motion and the twisting motion of the shortest chain segment. 
The increased activation energy is due to the fact that more energy has 
to be applied to the chain segment to break it loose from its equilibrium 
position when it is surrounded by adjacent polyisobutylene molecules 
than when it is surrounded by cyclohexane molecules. The stiffness of 
the chain is due more to the slope of the potential well than to any in- 
trinsic chain stiffness as is shown by Fig. 18, which shows the two stiff- 
nesses as a function of temperature. These values are obtained by fitting 
the 15°, 35° and 50° data in a similar manner to that used for the 25° 
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Fig. 17—Solid lines are mean values of relaxation frequencies for the two 
mechanisms plotted against molecular weight. Dotted lines indicate limits of 
regions assumed to obtain a better fit to the measured values. 
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data, and the activation energy of 11.2 kilocalories per mole is obtained 
in a similar manner. 

Due to the closeness of the surrounding polyisobutylene molecules, 
one would expect that the relaxation frequencies would not have discrete 
values but would be spread about the center value in some sort of a 
Gaussian distribution. If we approximate this by representing each 
region by two relaxations, one-half, and the other twice the frequency 
of the mean value, as shown by the dotted lines of Fig. 17, the agree- 
ment with the measured values of Fig. 15 is considerably better as 
shown by Fig. 19. A wider distribution yet is indicated. 

For molecular weights greater than 2,000 the shortest chain segment 
viscosity begins to diverge from the static viscosity indicating that 
there are other relaxations for these longer chains. Some data for the 
three longest chain polymers, F, G and H have been obtained by the 
torsional rod method and the results are given in Table II. These data 
are plotted on Fig. 20 as a ratio of dynamic to static viscosity plotted 
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Fig. 18—Variations of high frequency shear stiffness and ‘‘entanglement”’ 
stiffness plotted as a function of the temperature. 
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against the frequency times the static viscosity. All the viscosity data 
can be represented within the experimental error by a single curve, but 
the stiffness curves appear to require different curves for different tem- 
peratures. On analyzing the data in terms of a distribution of relaxation 
frequencies, a single curve for all temperatures could be obtained if the 
stiffness of each mechanism were independent of the temperature and 
the relaxation frequency were inversely proportional to the static vis- 
cosity, i.e., had an activation energy variation equal to that for the 
static viscosity for each relaxation mechanism. This condition holds 
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Fig. 19—Curves showing better fit to experimental data obtained by assuming 
the relaxation regions shown by Fig. 17. 


TaBieE II 














Srl 












































Measurement for Polymer F—Molecular weight = 4550 
15°C 25°C 35°C 45°C 
Frequency 
Kilocycles ——%_— > ie =|. <- - - — 
B 2 B 4 B ” B ” 
25 2.9 X 108 | 600 | 1.45 & 108 | 270 | 0.78 & 108 | 145 | 0.44 & 108 | 77 
40 3.8 X 108 | 550} 1.7 X 108 | 250 | 0.86 X 108 | 125 | 0.47 X 108 | 66 
62 4.0 X 108 | 470 | 1.9 X 108 | 220 | 0.78 XK 108 | 110 | 0.54 X 108 | 59 
140 5.8 X 108 | 380} 3.0 X 108 | 180 | 1.5 X 108 91 | 0.71 X 108 | 48 
Polymer G—Molecular weight = 5590 
24.5°C 65°C 
Frequency 
Kilocycles 
B 7 B 7 
27.5 2.0 X 108 | 510] 1.9 X 107 | 38.4 
40.96 2.95 X 108 | 480 
54.23 3.60 X 108 | 456 | 2.93 X 107 | 36.1 
Polymer H—Molecular weight = 10,380 
55°C 67°C 
Frequency 
Kilocycles 
B ” B ” 
32.2 7 | 150 y 94 


3.9 X 10 2.9 X 10 
40.5 4.8 X 107 | 187 | 3.5 X 107 86 
47.0 5.3 X 10 3.9 X 10 





cG6 AUVONVS ‘IVNUROL IVOINHOGL WALSAS TIGA AHL 


MECHANICAL PROPERTIES OF POLYMERS 149 


quite well for frequencies much lower than the relaxation frequencies 
of the smallest chain segment, but as the frequency approaches these 
relaxation frequencies, the stiffness of these polymers increases as the 
temperature decreases. 

A fair approximation to these measured values is obtained by assum- 
ing one more “configurational” relaxation frequency in addition to the 
two smallest segment relaxations discussed previously. Fig. 21 shows 
calculations of the ratio of dynamic to static viscosity and the shear 
stiffness for 65°C and 25°C. The lowest relaxation frequency is assumed 
to have a stiffness of 6.3 X 10° dynes/em’ and a viscosity of 20 poises 
at 65°C. For 25°C the stiffness of 2 X 10’ dynes/cm’ is assumed and an 
activation energy of 17.3 kilocalories gives the component a viscosity 
of 607 poises at 25°C, and a relaxation frequency of 5,250 cycles. The 
average value of 16 kilocalories for the static viscosity is a result of the 
sum of the variation due to the two components. Although the agree- 
ment can be improved by assuming distributions of relaxation fre- 
quencies centered around these three primary frequencies, there does 
not seem to be much doubt of the existence of these primary relaxation 
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Fig. 20—Plot of ratio of dynamic to static viscosity and the corresponding 
intermediate frequency shear stiffnesses as a function of temperature and prod- 
uct of frequency times static viscosity. 
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Fig. 21—Calculated values of stiffness and viscosity obtained by adding a 
single ‘‘configurational’’ relaxation frequency to the two short chain relaxations 
obtained from high frequency measurements. 
mechanisms which show up as discrete relaxations in long chain mole- 
cules in solution. 

Some measurements have also been made to determine the effect 
of chemical substitutions in the polymer chains. In a previous paper,* 
the high frequency properties of poly-a-methyl styrene were discussed. 
This material has the polyisobutylene chain but with one methyl re- 
placed by a phenyl so that its chain becomes 


ia 
se 
H 


For low molecular weights this material is liquid. The shear stiffness of 
this liquid is somewhat higher than for polyisobutylene but has about 
the same change with temperature. The variation of the high-frequency 
viscosity, however, is much larger for poly-a-methyl styrene than for 
polyisobutylene, and corresponds to an activation energy of 23.6 kilo- 
calories. The relaxation region for the shortest chain motion is much 


* See footnote on page 132. 
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narrower than in polyisobutylene, which correlates with the smaller 
steric hindrance. 

Measurements have also been made in the 70-kce to 140-ke range for 
two non polar liquids, polybutadiene and polypropylene and one polar 
liquid, polypropylene sebacate. The first two have the formulae 


a, i ee 
| 
ee a 
H H H H 
polybutadiene polypropylene 


while the polar liquid, polypropylene sebacate, has the formulae 


O H H H H HH H HH CH; H 
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The measured results are given by Table III. These data are plotted 
on Fig. 22 as a function of the product of frequency times static viscosity. 
For comparison the data for polyisobutylene polymer F is also plotted. 
By extrapolating to low values of the product-frequency times static 
viscosity it is seen that the low frequency “‘quasi-configurational’’ stiff- 
nesses of these liquids run from 10° to 1.5 X 10’ dynes per square centi- 
meter. Polyisobutylene has the greatest stiffness for any value of fre- 
quency times viscosity while polybutadiene has the least. This reflects 
the greatest steric hindrance of polyisobutylene and the smallest for 
polybutadiene which has only hydrogens connected to the carbon chain 
atoms. Another consequence of the larger steric hindrance of the CH; 
groups of polyisobutylene is that the viscosity associated with the short- 
est chain segment motion is largest for polyisobutylene and smallest 
for polybutadiene as can be seen from the ratio of dynamic to static 
viscosities for high values of frequency times static viscosity. 

The activation energy for static viscosity are for polypropylene, poly- 
propylene sebacate and polybutadiene respectively 21.2, 12 and 8 kilo- 
calories compared to 16 kilocalories per mole for polyisobutylene. The 
high-frequency activation energies as determined by the dynamic meas- 
urements are respectively 11.8, 4.6 and 1.4 kilocalories for polypropylene, 
polypropylene sebacate and polybutadiene. The differences between the 
activation energy for static flow and that for dynamic flow are respec- 
tively 9.6, 7.4 and 6.6 kilocalories, which values are all higher than 
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76.5 Ke Measurements 142.6 Kc Measurements 
Temp. °C | Density p vic Maxwell Maxwell 
1p/ns ing | t/t fns 
Np poise p dynes/cm2 | Np poise p dynes/cm?2 
Polybutadiene 
10 0.877 740 25.5 7.838 XX 10° | 0.0355 5.66 X 10? 14.65 | 9.96 > 106 | 0.0196 | 10.55 xX 107 
20 0.873 450 21.1 6.388 X 108 | 0.0469 3.44 xX 107 13.45 | 7.91 xX 108 | 0.0299 6.41 xX 10’ 
30 0.869 288 17.45 | 5.8 X 108 | 0.0605 2.21 x 10° 12.4 6.67 X 108 | 0.043 4.11 X 107 
40 0.865 189 15.3 4.67 X 108 | 0.081 1.44 X 10’ 11.05 | 5.96 X 108 | 0.0583 2.69 X 10’ 
50 0.861 128 11.9 3.92 X 10° | 0.093 0.98 xX 107 9.16 | 5.21 X 108 | 0.0715 1.82 X 107 
60 0.857 88 9.85 | 3.55 X 108 | 0.112 0.67 X 10’ 7.51 | 4.54 X 108 | 0.0855 1.25 X 10? 
Polypropylene 
76.5 Ke Measurements 142.5 Kc Measurements 
45 0.849 30,000 | 380 3.56 X 108 | 0.01265 | 2.8 X 109 | 354.5 6.69 X 108 | 0.0118 4.26 xX 10° 
55 0.842 10,900 160.2 1.4 X 108 | 0.01538 0.882 X 10° | 167 3.18 X 108 | 0.0153 1.55 xX 109 
65 0.836 4,200 89.8 0.786 X 108 | 0.0214 0.3821 * 109 79.1 1.46 X 108 | 0.0188 0.6 %X 109 
75 0.830 1,600 50.8 0.444 X 108 | 0.0318 0.128 X 109 46.6 0.998 & 108 | 0.0292 0.23 xX 10° 
85 0.824 700 31.2 0.264 X 108 | 0.0446 0.053 * 10° 28.4 0.51 XX 108 | 0.0405 0.10 xX 10° 
Polypropylene Sebacate 
77.1 Ke Measurements 141.5 Kc Measurements 
5 1.076 9,200 95.2 4.88 xX 107 | 0.0108 7.1 X 108 74.6 7.45 X 107 | 0.0081 1.38 X 109 
15 1.068 4,200 67.6 2.94 X 107 | 0.0161 3.24 X 108 54.9 4.28 X 107 | 0.0131 0.595 X 109 
25 1.060 2,060 51.4 2.01 X 107 | 0.0249 1.6 xX 108 41.3 2.82 X 107 | 0.02 0.292 * 109 
35 1.051 1,050 43.1 1.54 X 107 | 0.0411 0.81 X 103 33.2 1.97 xX 107 | 0.0316 0.148 X 109 
45 1.042 580 38.8 1.24 X 107 | 0.0669 0.45 X 108 26.9 1.69 xX 107 | 0.0464 0.082 X 109 
55 1.034 320 32.8 1.01 X 107 | 0.1028 0.25 XX 108 24.2 1.84 xX 107 | 0.0756 0.045 X 10° 
65 1.026 183 26.8 0.88 x 107 | 0.147 0.14 X 108 22.9 1.25 X 107 | 0.125 0.026 X 10° 
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Fig. 22—Ratio of dynamic to static viscosity and the shear stiffness for four 

polymer liquids plotted against product of frequency and static viscosity. 







RATIO OF DYNAMIC TO STATIC VISCOSITY 


the 4.8 kilocalories for polyisobutylene. This presumably indicates that 
there is more of a difference between the viscosity flow segment and the 
shortest chain segment in these materials than in polyisobutylene. Since 
no measurements are available over a range of molecular weights, no 
direct evidence has been obtained for the various chain lengths. 


B. Longitudinal Wave Measurements in Liquid Polymers 


Since the increase in shear elasticity for the highest relaxation fre- 
quency is so large, it should also appear in longitudinal wave measure- 
ments. Fig. 23 shows a calculation for the 5590 molecular weight liquid 
of the longitudinal velocity assuming that the Lamé } elastic constant 
is independent of frequency and that all the variation occurs in the shear 
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Fig. 23—Relation between measured longitudinal velocity for polyisobutylene 
of molecular weight 5590 and that calculated from shear stiffness measurements 
assuming the Lamé ) elastic constant is independent of frequency. 
constant as determined by the shear measurements. The points are 
velocities measured for longitudinal waves and as can be seen, the meas- 
urements agree closely with the calculated values. A slightly better 
agreement would be obtained if \ increased by a small amount as the 
frequency increased. As discussed in the next section there is some 
experimental evidence for an increase in A in nylon 6-6 and in poly- 

ethylene. 
~The question also arises as to how much of the attenuation is due to 
shear mechanisms and how much due to pure compressional effects. 
From longitudinal velocity and attenuation measurements at 30°C for 
the polymers E, F and G of Table I, the values of \ + 2u and x + 27 
can be determined and are shown by Table IV. The values of u and 7 
can be obtained from Table I by interpolation and are given in columns 



































Tasue IV 

Polymer + 2u dynes/cm? ae p» dynes/cm? 7” poise. \ dynes/cem? x 

5 megacycles 
E 1.92 X* 101° 62 0.17 X 101° 26 1.60 * 101° 10 
F 2.26 107 0.23 & 10% 46 1.70 X 101 15 
G 2.38 170 0.31 X 10% 75 1.76 X 10% 20 

8 megacycles 
E 2.01 X 101° 50 0.20 X 101° 22 1.61 X 101° 6 
F 2.26 93 0.27 * 10% 41 1.72 XK 10% 11 
G 2.56 155 0.34 * 101 69 |. 1.88 * 10” 17 
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4 and 5. Columns 6 and 7 show the values of \ and x, the compressional 
components. A definite longitudinal compressional viscosity is indica- 
ted which however is somewhat smaller than the shear viscosity 7. 


V. MEASUREMENTS FOR SOLID POLYMERS 


Recently two new methods have been devised for accurately measur- 
ing the properties of solid plastics in the ultrasonic frequency range 
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Fig. 24—Normal loss in db per centimeter measured as a function of frequency 


for several rubbers and plastics. 
moze] 





Fig. 25—Debye-Sears cell for making sound waves visible. 
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Fig. 26—Examples of refraction and focusing effects for sound waves. 


and these have shown relaxations in such plastics as polyethylene and 
nylon 6-6. The simplest method for measuring one of the properties for 
longitudinal waves, i.e., the attenuation, is to measure the change in 
loss between two transducers in a liquid such as water, caused by insert- 
ing a sheet of the material. This process, used during the last war, 
results in the losses in db per centimeter for several rubbers and plastics, 
shown by Fig. 24. Indications of relaxation mechanisms are given by 
the rubbers and the plastic tenite II which is a cellulose acetate butyrate. 
The first fairly accurate method for measuring longitudinal sound 
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velocities” in plastics was the method of observing the focusing effect 
of a cylindrical lens made of the plastic. Sound waves can be made 
visible by the Debye-Sears technique of using a sound wave as a phase 
diffraction grating. Here light from a slit S, is made parallel by the 
lens L, and passes through the cell parallel to the wave fronts of the 
sound waves as shown by Fig. 25. The compressed parts of the medium 
retard the light waves more than the rarefied parts do and hence the 
medium acts as a phase diffraction grating. If a second slit S_ is used 
which is small enough to pass only the zero order, a light valve action 
is obtained which modulates the light according to the sound wave 
intensity. If now the lens LZ; is used which focuses on the median plane 
of the tank, a picture of the sound beam is obtained as shown on Fig. 26. 
The bottom figure shows the focusing effect of a plastic and from the 
focal distance d and the radius of curvature r of the lense, one can cal- 
culate the velocity in a plastic compared to the velocity in the water 


by the formula 
r 
Up = we / (1 —- *) (11) 


This method gives velocities good to from 2 to 5 per cent depending on 
the attenuation in the lens. 

G. W. Willard’® has devised recently a more accurate method for 
measuring sound velocities as shown schematically by Fig. 27. Here a 
plastic to be measured is placed half way across the sound beam in the 
liquid and light is sent along the wave front occurring in both the plastic 
and the liquid. If the waves are in phase the retardation in the two light 
gratings, corresponding to sound propagation in both media, add up 
and for a slit selecting the zero order the darkest pattern occurs on the 
photographic plate. If the two waves are just out of phase, the retarda- 
tion is reversed in the two media and the lightest part occurs. With 
this relation it can be shown’’ that the spacing d of light and dark lines 


122W. P. Mason, Piezoelectric Crystals and Their Application to Ultrasonics, 
D. Van Nostrand, 1950, p. 404. It was used in this country by G. W. Willard as 
early as 1940. It was also used in Germany by J. Schaefer ‘‘Eine Neue Method 
Zur Messung der Ultraschallwellen in Festkorpern.”’ Diss Strassburg, 1942. By 
making the front surface part of a cylinder, Schaefer also measured the shear 
velocity in a solid. 

1G. W. Willard, J. Acous. Soc. Am., 23, Jan. 1951, pp. 83-94. The origin of 
this multiple path interference method goes back to the work of R. Bir (Helvetia 
Physica Physica Acta Bd 13 page 61 (1940)) who attached a piezoelectric crystal 
to a bar with a 45° end section and set up transverse and longitudinal waves, in 
the bar. These waves produced longitudinal waves in a surrounding liquid and 
by observing the interference pattern between them, the longitudinal and shear 
constants could be determined for an isotropic medium. Willard’s method as 
described above is much more direct and is capable of higher accuracies. 
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is related to the wavelength in the liquid \, and the wavelength in the 
solid As, by 


ft ee (12) 


This corresponds to a velocity in the solid compared to a venOeIy in 
the liquid given by 


Ue 


1 — vt/fa 


where f is the frequency. Fig. 28 shows a photograph of a series of lines 
in a transparent plastic and a transparent plastic in the form of a wedge. 
It is seen that beyond the edge of the plastic there is a dark interference 
band for each one in the transparent plastic. This phenomenon is caused 
by the refraction of the sound wave that has traversed the plastic and 
the dark lines are lines of equal phase of the two waves in the liquid. 
The angle of the dark lines is half the refraction angle. Hence the velocity 
can also be determined by counting the number of dark bands in the 
liquid beyond the plastic. This makes it possible to measure the veloci- 
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Fig. 27—Optical method for measuring sound velocities of plastics by com- 
paring their velocity with that of a liquid such as water. 
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“Ui = 2fs-sin @ 





Fig. 283—Photographs of interference patterns from sound waves in liquids 
and plastics. 


ties in opaque plastics. The accuracy of the method is better than 1 per 
cent if the attenuation is low enough to give a number of interference 
lines. For plastics of high internal loss, the method becomes somewhat 
inaccurate. 

Typical measurements using this system are shown in Table V. Small 
changes in chemical composition and plasticizer content are shown up 
as can be seen from the table. Of particular interest is the difference 
between nylon 6-6 and polyethylene. Chemically as shown by Fig. 29, 
the two are identical except for the dipoles occurring for every 6 units 
of the ethylene chain. These dipoles have the effect of bonding adjacent 
layers together and result in a higher shearing modulus. 

By attaching shear vibrating crystals to a right angled prism, as shown 
by the lower part of Fig. 28, with the direction of motion of the crystals 
parallel to the transmitting face, and setting up shear standing waves 
between the two crystals, the shear properties of the plastic can be 
measured. Longitudinal waves are generated in the liquid which inter- 
fere with one another and cause dark bands perpendicular to the plastic 
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TABLE V 


- Measured longitudinal velocity and attenuation at 25°C and 2.5 me. 
Longitudinal attenuation in DB/cm at 25°C and 2.5 mc except as noted. 





; Long | Shear : 
Material : are velocity. A DB/cm | Density 
; cm/sec cm/sec ; 
Dural Ve SE hac yrce- vi steeds Vee Ree 6.5 3.12 _ EY (ite: 
Brass, half hard..................0. 000005 4.7 2.11 — — 
Polystyrene «sis 2 seen te eich ees 2.35 1.12 | 2 1.05 
PION BIAS ost oh ons ecne ede hae eae eat esas 2.68 5 1.18 
Tenite II, (cellulose acetate butyrate), 2% 
plasticizer... 0.0... eee ees 2.08 9 1.23 
Tenite II, 18% plasticizer................. 2.02 10 1.21 
Polyvinal formal. ...................000005 2.68 10 1.24 
Polyvinlidene chloride............... -.ee.| 2.4 18 1.71 
Poly N-butyl methacrylate............... 1.96 5 1.05 
Poly I-butyl methacrylate................ 2.08 6 1.05 
NEOPrene...... 0... ccc ec cece et eee ee enes 1.51 20 .99 
Polyethylene soc. 026 hock cevnn se eh wees Meets 2.0 4.7 Fru .90 
Nylon 6-6 (8-80 megacycles).............. 2.68 1.0 F1-5 1.11 
Nylon 6-10 (8-30 megacycles)............. 2.56 1.0 F!-5 1.11 


surface. By determining the spacing of these lines the velocity of the 
shear waves can be determined. 

Another method has also been developed which is more applicable for 
high loss materials. This is a pulsing method and is a modification of the 
method proposed by one of the authors for measuring the properties 
of small crystal specimens." Here longitudinal or shear crystals are sol- 
dered to the fused quartz rod as shown by Fig. 30 and a sample to be 
measured is placed between these by means of a liquid such as poly- 
isobutylene which has a high shear elasticity. If the specimen has a 
small attenuation, this can be measured by taking the difference in the 
amplitude of successive reflections. If the specimen has a high loss, this 
does not work and another method has been used which consists in 
sending a pulse from both crystals.’° One crystal is then used to receive 
and it receives the wave sent through the sample and the wave reflected 
from the fused quartz-sample interface. By adjusting the amplitude 
until these two are equal and the frequency or phase of one channel 
until the waves cancel, a ratio of amplitudes and a frequency of half 
wavelength are accurately determined. From these the velocity and 
attenuation can be calculated. 

This method has been applied to measuring the longitudinal and shear 
velocities of polyethylene and 6-6 nylon. The polyethylene was of ‘‘equi- 
librium” crystallinity and average molecular weight corresponding to 

144. J. McSkimin, ‘‘Ultrasonic Measurement Techniques Applicable to Small 


Solid Specimens,’’ J. Acoust. Soc. Am., 22, No. 4, July 1950, pp. 413-418. 
18H. J. McSkimin, J. Acoust. Soc. Am., 23, No. 4, pp. 429-435. 
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an intrinsic viscosity in xylene of [n] = 0.89 at 85°C. Fig. 31 shows the 
longitudinal velocity of polyethylene plotted as a function of frequency 
and temperature. The velocity rises with frequency and a dispersion 
is indicated. This is confirmed by the attenuation per wavelength curve 
for two different frequencies plotted as a function of temperature, Fig. 
32. A definite dispersion is seen to occur with an activation energy of 
about 12 + 2 kilocalories per mole. This could occur in either the A 
constant or the shearing constants yu, but the data of Figs. 33 and 34 
show definitely that it occurs in the shear component. Fig. 33 shows 
the shear velocity for four temperatures plotted as a function of fre- 
quency. This can be fitted for 30°C with a single relaxation mechanism 
having a relaxation frequency of 8 megacycles. To agree with the meas- 
ured attenuation and velocity, there has to be a spreading of the single 
relaxation over a range as also occurs in liquids. The indicated shear 
stiffness below this relaxation frequency is 2.6 X 10° dynes/em’. Some 
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Fig. 29—Spatial structure of nylon 6-6. 
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data on the zero frequency shear modulus is obtained from the Young’s . 
modulus for a static pull which is from 30,000 to 50,000 pounds/square 
inch. Since the Young’s modulus is three times the shearing modulus, 
the zero frequency shearing modulus should not exceed 1.1 X 10° dynes/ 
em’. Hence one may expect that other relaxations will occur at lower 
frequencies. | 

Fig. 34 shows the attenuation per wavelength for shear waves. The 
solid line for 30°C represents the calculated attenuation per wavelength 
for the model assumed. If all the dissipation were due to shear mech- 
anisms, the calculated attenuations would occur as shown by the 30°C 
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Fig. 30—Ultrasonic pulse method for measuring the velocities and attenua- 
tions of highly attenuating plastics. 
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Fig. 31—Velocity of longitudinal waves in polyethylene plotted as a function 
of temperature and frequency. 
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Fig. 32—Attenuation per wavelength for longitudinal waves in polyethylene 
plotted as a function of temperature and frequency. ; 
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Fig. 33—Shear velocity of polyethylene as a function of frequency and temper- 
ature. Equivalent circuit shows elements necessary to account for the velocity 
and attenuation changes at 30°C. 
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Fig. 34—Attenuation per wavelength for shear waves in polyethylene. 
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Fig. 35—Velocity of longitudinal waves in nylon 6-6 plotted as a function of 
temperature and frequency. 
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Fig. 36—Attenuation per wavelength for longitudinal waves in nylon 6-6 
plotted as a function of temperature and frequency. 
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points of Fig. 32. Most of the loss is accounted for by shear mechanisms, 
but it appears that some compressional mechanisms may also be present. 

The mechanism causing the relaxation in the megacycle range for 
polyethylene appears to be the same as for polyisobutylene, namely 
the relaxation of the shortest chain segment that is free to move. The 
chain segment acting appears to be longer than six chain units for similar 
measurements of nylon 6-6 show no relaxations in this frequency range. 
Fig. 35 shows the longitudinal velocity and Fig. 36 the attenuation per 
wavelength for longitudinal waves. Since the attenuation per wave- 
length is still increasing for nylon 6-6 at 25 megacycles a still shorter 
chain segment may be operating for this material. The shear velocity 
and attenuation per wavelength for nylon 6-6 are shown by Figs. 37. 
and 38. 

Fig. 39 shows the shear stiffness of polyethylene and nylon 6-6 plotted 
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Fig. 37—Velocity of shear waves in nylon 6-6 plotted as a function of temper- 


ature and frequency. 
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Fig. 38—Attenuation per wavelength for shear waves in nylon 6-6 plotted as a 
function of temperature and frequency. 
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Fig. 89—Shear elasticity of polyethylene and nylon 6-6 plotted as a function 
of temperature and frequency. 
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Fig. 40—Value of Lamé \ elastic constant for polyethylene and nylon 6-6 
plotted as a function of frequency and temperature. 
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against 1/7’ where T is the absolute temperature. Both are plotted for 
8 me and 30 me. The dispersion in both materials is evident. Below 
30°C the shear elasticity of polyethylene varies exponentially with the 
temperature with an activation energy of 2.72 kilocalories per mole. 
Above this temperature a deviation occurs due to the approach to the 
melting temperature. Nylon has a smaller variation with temperature. 

Comparing the longitudinal and shear wave measurements one can 
calculate the Lamé ) elastic constant and this is shown plotted on Fig. 40 
for both polyethylene and nylon 6-6 as a function of temperature for 
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Fig. 41—Equivalent shear and compressional viscosities for polyethylene and 
nylon 6-6 plotted as a function of frequency for a temperature of 25°C. 


two frequencies. The dispersion of \ for polyethylene is small but is more 
prominent in nylon 6-6. This correlates with the larger compressional 
viscosity component present for nylon 6-6 which as shown from Fig. 41 
is as large as the shear viscosity. According to the structural rearrange- 
ment theory of compressional viscosity due to Debye,’® compressional 
viscosity can enter when some part of the chain can rearrange from one 
stable state to another stable state as a function of pressure. This re- 
arrangement occurs across a potential barrier and hence requires a finite 
amount of time to occur. This lag in the rearrangement results in a— 
compressional viscosity and as the frequency is increased, a frequency | 
is found for which the motion can no longer occur in the time of a single 


16 P, Debye, Z. Elektrochem., 45, 1939, p. 174. 
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cycle and the \ constant increases. It appears from these measurements 
that the dipole binding present in nylon 6-6 allows a greater structural 
rearrangement under pressure than can occur for polyethylene which 
has only linear chains. 


VI. CONCLUSIONS 


Measurements in dilute solutions, in pure polymer liquids and in non 
rigid solid polymers have all shown the presence of a shortest segment 
whose relaxation leads to a crystalline type of elasticity. In dilute 
polymer solutions the presence of a configurational type of relaxation 
and an entanglement relaxation of the shortest chain segment have 
been shown. For pure polymer liquids a quasi-configurational type of 
relaxation has been found for chain lengths greater than 60 segments, 
but for chain lengths less than 40 segments this type of relaxation dis- 
appears. From the difference between the high frequency shear elastici- 
ties measured for polyethylene and nylon 6-6 and the static measure- 
ment of Young’s modulus, it appears that there may be other relaxations 
in these materials for lower frequency ranges. . 

For pure polyisobutylene and for nylon 6-6 there appear to be struc- 
tural changes induced by pressure which account for a compressional 
viscosity and a dispersion in the d elastic constant. This effect is smaller 
for polyethylene. 


APPENDIX——EFFECT OF LIQUIDS ON THE PROPAGATION OF SHEAR WAVES 
IN RODS 


For radially symmetric rods, the tangential particle displacement. 2, 
in the rod is given by 


Up = Jr(kr)e?* (1) 


where 
2 
=u Pa ¢ (1A) 
Le 


All other displacements are zero. In this equation waves are considered 
to be travelling in the +z direction with a propagation constant @ = 
A + 7B, where A is the attenuation in nepers per cm and B the phase 
shift in radians per centimeter. u is the shear stiffness which may be 
complex to take account of the dissipation within the rod. 


MECHANICAL PROPERTIES OF POLYMERS 169 


From the defining relations for the stress strain equation 


ou U, 
Te = pS» = a (Se = “) 
or Tr 
and the tangential particle velocity du,/dt, one can calculate the im- 
pedance Z per square cm. of cylindrical surface at r = a. This relation is 


_ af jpk Be 2 
sma Olan eral ar 7 a @) 


Since only the first mode is excited, parameters can be adjusted to 
keep k quite small, i.e. (ka < .2) and equation (2) can be simplified by 
using power series expansions for the Bessel functions. Neglecting higher 
order terms this results in 





, T= —juak” (3) 
4w 
To evaluate the impedance of the liquid surrounding the rod, the tor- 
sional wave is first propagated along the length of the rod without the 
liquid, i.e. with Z = 0. Then from equation (8) k = 0 and from equa- 
tion (1A) 





ee 2 
05 = 7 = (Ay + jBo)’ (4) 


where Ap and Bo are respectively the attenuation and phase shift in 
the rod alone. With the small loss in metal and glass rods Ao can be 
taken equal to zero and 


Bo = w/WS n/p =, w/Vo (5) 


where v is the velocity of propagation in the rod alone. 
When the liquid surrounds the rod, however, 


2 
k? = oe + 6 = —(Ao + jBo)? + (A + jBY (6) 


For the usual case where (B + Bo) > (A + Ao), equation (6) approxi- 
mates 


k? = (B + Bo) (— AB + jAA) = 2By (—AB + jAA) 


where AB is the increase in phase shift per centimeter and AA the in- 
crease in attenuation per cm, both directly measurable quantities. The 
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final working equation is then given by 


Z ="" x 2BIAA + JAB] = E° (AA + JAB) =" (AA +jMB) (7) 
4a 29 2 
Since AA and AB are the attenuation and phase shift changes per unit 
length, then if f is the length of the rod, covered the total attenuation 
and phase shift changes will be AA and AB multiplied by 2. Hence if 
AAy and AB are the measured attenuation and phase changes, the 
impedance Z becomes 
i= ae (AAo + jABo) (8) 
This derivation neglects the change of phase occurring at the inter- 
section between the rod having no liquid and the rod surrounded by 
the liquid, but it can. be shown that this is small and moreover, the 
change in the wave on leaving is equal and opposite to that occurring 
on entering and hence this correction cancels out. However, if the 
liquid is viscous enough there is a correction due to the fact that the 
measured impedance of equation (8) is for a cylindrical surface, whereas 
the desired impedance is the characteristic plane wave impedance. 
Obviously if the radius of curvature is sufficiently large no correction to 
equation (8) need be made. To obtain a suitable criterion one may con- 
sider waves propagated into the liquid from the surface of the rod and 
solve for the impedance per square cm. of the cylindrical surface. This 
neglects the variation with z, but since the wavelength along the rod is 
quite large, little error results from neglecting variations with z. 
An outgoing cylindrical wave in the medium may be represented by 


Ug = [Ji(kr)’ oy 7V (kr) Je"! 7 H® (kr)! eit (9) 


where the primes refer to the wave in the liquid and 





, 2 2 2 , 
a ee a oe A (10) 
k 


where Z;, = ~/y’p’ is the plane wave impedance of the liquid. 
The shearing stress 


Ou U 
Tro = wSe = pw |—-— ; 
or Tr 
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and the tangential velocity may be obtained as before and the complex 
impedance over the cylindrical surface determined to be 


~ (2) 1 (2) June? 
Gow 8 = Hy a Lue ie ©) | / joH}” (kr)’ (11) 


U3 or 
Noting that for plane waves Z, = ~/p'u’, one may eliminate yp’ from 
(11) and evaluate Z at r = a with the result 
| AS? (ka)’ y | 
Z= | ae —- |Z (12) 
’ LHP (ka)’— (ay’ |" 


The above equation can be used to obtain a solution for Z;, in terms of 
the measured value of the cylindrical impedance Z for any set of para- 
meters that may apply. Except for very heavy loading of the rod, how- 
ever, the results of equation (8) may be used directly with little error. 
For example calculations indicate that for | ka|’ = 10, the multiplier 
of Z;, of equation (12) is approximately 1.07 Z — 7.4° for phase angle of 
(ka)’ of —25°. For (ka)’ < 10, the correction multiplier rapidly becomes 
important. This same correction is applicable for the torsional crystal, 
but since this is only used for dynamic viscosities less than 10 poises, a 
correction is seldom necessary. 


Relay Armature Rebound Analysis 


BY ERIC EDEN SUMNER 
(Manuscript received October 25, 1951) 


Rebound of mechanical structures subsequent to impinging on stops gen- 
erally has deleterious effects on their performance and should, therefore, be 
minimized. A considerable reduction in rebound can often be obtained by 
introducing additional degrees of freedom to the structure. 

A mathematical treatise of the dynamics of rebound motion of systems 
representing idealized relay armatures is presented. Normalized differential 
equations of motion and their solutions for the “free” and “impact” inter- 
vals are derived for systems having one, two, and three degrees of freedom, 
allowing the rebound behavior of a specific system to be calculated. The equa- 
tions of series of rebounds, and possible combinations of such series are con- 
sidered next for systems having one and two degrees of freedom. The field of 
possible rebound maxima is mapped for a practical range of mass distribu- 
tion constants, coefficients of restitution, and force ratios. A sufficiently broad 
optimum design region ts indicated. 

The results of this analysis have been checked closely on a model and have 
led to appreciable reduction of armature rebound in relay designs. 


I. INTRODUCTION 


In numerous types of mechanisms it is desirable to arrest the motion 
of a member at a particular point and to maintain it in this position. 
One of the simplest means of accomplishing this is to allow the moving 
member to impinge on a fixed member (stop) and to provide forces to 
tension it against this stop. Because the member to be arrested possesses 
kinetic energy and because the stop cannot generally absorb all of this 
energy, the moving member will rebound from the stop. The rebound 
motion generally deteriorates the performance of the mechanism and 
should be minimized. 

Investigation of this phenomenon has been stimulated by the armature 
rebound problem in relay operation, where rebound from the front stop* 
tends to reclose contacts and must therefore be compensated for by 
additional (waste) travel, resulting in deleterious effects on speed and 
~ * Among relay designers the front stop has been generally referred to as “‘back- 


stop’’. In this paper the terms front stop and heel stop have been used through- 
out for easier identification. 


172 


RELAY ARMATURE REBOUND ANALYSIS 173 


magnetic characteristics. Analysis in this paper will be directed towards 
relay armature systems, but it is also applicable to rebound in similar 
mechanisms. 


II. STATEMENT OF PROBLEM 


Analysis will be restricted to planar motion of armature systems 
having one, two, and three degrees of freedom as depicted in Figs. 1, 
2, and 3. Generally one stop must be provided for each degree of free- 
dom, although in the three-degree-of-freedom system of Fig. 3, two of 
the stops have been combined. 

Applied forces Ff; , F,, F'3, have been chosen so as to be most easily 
correlated with actual relay designs. 

The initial condition in all cases will be a pure rotation about the 
“heel” just prior to a “zero” impact at the “front” of the armature. 
The “zero” impact will be followed by rebound motion and impacts at 
the various stops eventually bringing the armature to rest. The object 


FRONT STOP YW 


FRONT STOP. 


fa 
ARMATURE 
HEELSTOP__ ——<—_*. 
CORE a Fi 
Fig. 1—Solidly hinged armature— Fig. 2—Loosely hinged armature— 
one degree of freedom. two degrees of freedom. 
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Fig. 3—Armature system—three degrees of freedom. 
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will be to minimize rebound motion at the front, since this is usually 
near the point actuating the relay contacts. 

The basic problem is then to find the response of the armature subject 
to aperiodic but well defined impulses, which are functions of the 
positions and velocities of the system. 


lil. ASSUMPTIONS 


In order to facilitate the solution of this problem,.the following 
modifying assumptions are made: 

(1) As mentioned in the previous section, analysis is restricted to 
planar motion. 

(2) The armature is assumed to be a rigid body. 

(3) Stops are assumed to be very stiff, massless springs capable of 
energy absorption during impact with the armature. The associated 
coefficient of restitution is assumed constant. Core and stop vibration 
are neglected. 

(4) The tensioning forces F; , F2 , F3 are assumed to be constant forces. 
(This is fairly closely true for moderate rebound amplitudes of practical 
relay structures.) 

(5) All displacements are small relative to the dimensions of the 
system and in particular the angular displacement 6 is sufficiently small 
so that 


cos @ = 1 


sin 6 = @ 


IV. DERIVATION OF EQUATIONS OF MOTION 


The derivation of the equations of motion resolves itself into the 
solution of two different types of intervals: 

(1) Free Interval: This is the period during which the armature is 
not in contact with any of its stops and only the tensioning forces are 
acting. 

(2) Impact Interval: During such intervals the armature is in contact 
with at least one of the stops. The stiffness of the latter is assumed so 
high that the tensioning forces during this interval may be neglected. 

The three-degree-of-freedom case will be considered first and the 
others subsequently deduced from it by allowing some of the constants 
to approach zero. 
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A. Free Interval 


The motion of the armature will be described by the displacement at 
the stop points: 21, 23, 23.* Let m be the mass and RF the radius of 
gyration of the armature about the center of gravity. The latter is 
located by the dimensions 4,R, &R, and &R relative to the stop points, 
i.e., the points on the armature which contact the stops in the rest 
position (Fig. 3). 

The equations of motion are derived in Appendix I and are put into 
dimensionless form: 





sailors] rm 
w—Bfoe + cold eco] (Ym (2) roe] 0 


= i (F3) ; t : 
y= 5 | Cs + Cy (F) + C33 a (: y + Yz0 (*) + Yz0 


where: 





hee Mi, sy, SE oe He 
Yi am aa ge : 
T (2) 
= tam 
Fy 


tq is the front velocity #1, just prior to the “zero” impact, and 


Cu an (Ci ae 1) Ci a Cs = hh 
Cy =(G+1) Cy = Ca = (1 — 46) (3) 
C23 = (f5 a 1) C23 = C2 = — bk 


Yi0, Yoo, Yo, are the initial velocities and yyw, Yoo, Yso the initial dis- 
placements for the free interval in question. 

The equations of motion for a two-degree-of-freedom system are 
obtained, if #3; = 0. Then for the two coordinates of interest: 


3 | Cu =A ia Cp (#)| (: ) + Yro (: ) + Yio 
1 
Yo = ; [oa 4 iG Co (*)| () “+ Yoo (‘) + Yoo 


* A summary of all notations used in this paper is given in Appendix IV. 


gus 


Yi 
(4) 
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For a one-degree-of-freedom system gz = Ca + Cx (7) = 0, whence 
1 


2 2 
Y= 5 | Cu at Ca (‘) + YH (‘) + on | (5) 


B. Impact Interval 


69) 
7 


The change of velocity at point due to an impact at “2” is, by 
definition of the coefficient of restitution ‘k’’, 


Az, Savor (1 + ki)ay 


It is assumed here that the action of the stops are true impacts, i.e., 
the changes in velocity take place while there is negligible motion of the 
body. The velocity changes then occur as instantaneous rotation about 
the conjugate axis, leading to the general relation for an impact at 
point “‘2’’: 


Yion = Yien—y + KyAien—» (6) 


The first subscript indicates the coordinate, the second subscript 
indicates the beginning (0) or the end (e) of the free interval described 
by the third subscript. The impact transfer coefficient K;; relating a 


con) 6699 
Ves 


velocity change at point “7” to an impact at point 


oe (+ id (7) 

Equations (1) through (7) allow any one specific case to be mapped, 
if the mass distribution and force ratio are known. A sample of such 
mapping of rebound motion for a rectangular two-degree-of-freedom 
armature appears in Fig. 4. 


Ky = = 





V. ANALYSIS OF REBOUND PATTERN-——ONE-DEGREE-OF-FREEDOM SYSTEM 


The rebound pattern for the one-degree-of-freedom system—as de- 
rived in Appendix II—consists of an infinite series of parabolic ares of 
diminishing amplitudes. The structure comes to rest after a finite time 
interval. The maximum rebound occurs during the first bounce and 
equals 

k? 
oer (8) 


RELAY ARMATURE REBOUND ANALYSIS 


where 
Che 
CS Ca 
ll Coo 
The system returns to rest at 
is 2 
r C(l—k) 
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(9) 


(10) 


VI. ANALYSIS OF REBOUND PATTERN—TIWO-DEGREE-OF-FREEDOM SYSTEM 


The reason for choosing a two-degree-of-freedom system over a one- 
degree-of-freedom system would be, in keeping with the philosophy of 
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Fig. 4—Front and heel motion of plate type armature. 
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this treatment, to reduce Y,, the greatest excursion at the front. In 
order to simplify mapping, this maximum excursion will be expressed 
as 2CY,, the ratio of Y; to Y as given by Equation (8) for the case of 
k = 1. Thus 2CY, is the ratio of the greatest excursion of the two- 
degree-of-freedom system under consideration to the greatest excursion 
of the corresponding perfectly elastic one-degree-of freedom system. 

We first introduce two basic constants which are functions of the mass 
distribution relative to the stop locations: 


Ce 
Ci; 





M;; = (11) 
This constant represents a mechanical coupling coefficient. As 
M;; = M;;, the two-degree-of-freedom system under consideration here 
has only one such non-trivial constant My . 

The second constant represents a force transformation factor from the 


“7” coordinate to the ‘“2”’ coordinate: 
Cs; 


Py i: Cs (12) 
In the analysis of the two-degree-of-freedom system only P;2 is important. 

If there is to be any heel motion, the “‘zero”’ impact at the front must 
impart a positive velocity to the heel. By Equations (6), (7), and (12), 
this requires that Py. be negative, which in turn implies that Gh > 1. 
For the limiting case of 42 = 1, Pe = My = 0 and no coupling exists 
between the heel and the front. Physically this means that the two 
stops are the centers of percussion of each other and the system will 
act as a simple hinge. 

With the above foundation, it is possible to analyze the patterns cf 
motion and maximum rebound amplitudes. 


A. Motion Immediately Following ‘‘ Zero” Impact 


After the ‘‘zero” impact at the front, both front and heel will lift off 
in accordance with impact Equation (6) and continue to move in ac- 
cordance with the free interval Equations (4). Whether the next impact 
occurs at the front or the heel depends on their respective periods, 
t, and f2 : 

Pr 


ip. 1+ Pof 1+ hy 
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where: 


A large value of t/t. will result in a series of heel impacts and the 
heel will come to rest while the front is still displaced from the stop. | 
This will be called a complete heel series. A small value of t/t results in 
a similar complete front series. If 4/t, is near unity, a limited number of 
impacts on one end are followed by an impact on the other end, etc. 
An analysis of front and heel series follows: 


B.. Front Series 


If t/t < 1a series of front impacts occurs. The impact velocities at 
the front are 


Yin = 1, k, b?, +--+, ke” (14) 
The corresponding time intervals are 
Qh Qk} 2k” 
aay ow cae (15) 
where 
A= (Cu a Cf) 
During this time, the heel velocity and displacement are given by 
; wes 2B 
Yoon) = Yoon + E = Pll a i) [on 
(16) 


2B 
Yointt) = Yoon + E Yin + inn ian 
where 
B= Cy + Cof 


The velocity and displacement at the heel after a given number of front 
impacts are obtained by a summation of Equations (16). For a com- 
plete front series n > ©, and 





_ Qhy E _p | 
Yre AG + hy? + hy? a 12 es 
. 1 2Bk 
Yre0 = = a A oa Py + i) | 


180 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1952 


In addition it is useful to set down energy equations in order to 
simplify evaluation of greatest rebound for the various groups of re- 
bound patterns. The kinetic energy function T is evaluated in Appendix 
I. A potential energy term V—the work done against F; and F: from the 
equilibrium position—is introduced. If 7» is the total energy of the 
system prior to the “zero” impact, then 


oo gi eae a 
Ts a peer PB . (18) 
-—2C (yr + frye) : 


The energy loss due to n front impacts is 











-a (FEY) | a — et  - ati (19) 
0 
For a complete front series n — «, and 
(FEY) 0 Mudie 20) 
0 
If a complete front series follows the ‘‘zero” impact, yio = 1 and 
A é r r) = (1 — My) (21) 
0 . 


After completion of this ‘initial’ front series, the system maintains 
only one degree of freedom (rotation about the front) until a heel 
impact occurs. By setting 71 = y: = ye = 0 in (21) we obtain the heel 
impact approach velocity ye = Py. 

Apparently energy loss due to n front impacts is a function of Mi , 
ky , and the approach velocity of the first impact. 


C. Heel Series 


An analysis similar to the above can be made for partial and com- 
plete heel series following the ‘‘zero” impact. This is demonstrated in 
Appendix III, yielding, for ky = k,* 





_ AP + aE kd = %) _ | 
RE Ra yet Be aR = ss 
_ ._l1+k[24Pr ka —k) | 
Yiewo = i af B GaeKy aaa My + i) | 


* The more general form k, ¥ ke can be obtained as indicated in Appendix III. 
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The energy relationships for heel series are 


Lik t+ ”) Sa 2") Mi(1 a My) ae (23) 
To Pie 





For a complete series n — , and 


a, (7 as ”) = M,(1 — Mi:) Yoco (24) 
0 


12 
If a complete heel series follows the “zero” impact, 20 = Pu(1 + kh), 
and 





i (Ft r) ey) One eee a (25) 


0 


Finally, for the special case where a complete heel series follows an 
initial complete front series Yeo = Pi, and 


T+V\_ 7 | 
-a( 7 ) = M1,(My — 1) (26) 


0 





It is to be noted that the energy loss due to a partial heel series is a 
function of Mi: , Pie, ke , and the approach velocity of the first impact, 
but that the equation for a complete heel series does not contain k, . 
Finally, a complete initial heel series is a function of only My and k, . 


D. Complete Mapping of Problem 


Equations (1) through (26) make it possible to completely map the 
two-degree-of-freedom rebound problem. The relative maximum ampli- 
tude 2CY, and the rebound pattern will be determined. 

Examination of the necessary equations, show that 2CY; is in all 
cases a function of four parameters: k, , ko, Mi and Pf. Of these, ke 
enters only if a partial heel series occurs prior to the time of maximum 
rebound. If it is assumed that for this limited group of cases kz = k; = k, 
the number of parameters is reduced to three: k, My, Puff. 

In Figs. 5 to 10, 2CY; is plotted against Pyf for the most useful range 
of 1/8 os My < 1/2.5, 0.38 <k < 0.6 and0 < Prof < 10. 

_As Pf is increased from zero to infinity (corresponding to an increase 
in the heel tension 2), the rebound pattern goes through some or all of 
five regions. The criterion for location in any one region is based upon 
the parameter 





1+ Perf & &k (27) 
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Region 1—Complete initial front series for 1 < Q < 1/k. Within this 


2 
region, if Mj. > See the maximum rebound occurs during the 
first bounce and 
(1 — Mi)k? 
2cy, = ————_— 28 
1 TF Pof (28) 


If the maximum rebound occurs later, it must occur during a com- 
plete heel series which follows the initial complete front series. From 
Equations (21) and (26) 


By comparing Equations (28) and (29), the critical requirement 
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Fig. 5—Relative maximum rebound amplitude for My. = 1/2.5. 
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Cee sud 2) ie 
12 
while the first rebound maximum, shown in solid lines on Figs. 5 to 10, 
is always realized, the later rebound given by (29) is an upper limit— 
shown in dashed lines—and is not always realized. In the dashed re- 
gions, phasing is extremely critical and small variations in the param- 
eters may cause large variations in maximum rebound. From an 

engineering standpoint these regions are essentially undesirable. 
Region II—Partial initial front series for 
1 l+k 
k Os eh 
This region is one of critical phasing, and attention is limited to 
special cases leading to maximum rebound. These cases occur when a 


for the latter case is that Prof > —1. It should be noted that 
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Fig. 6—Relative maximum rebound amplitude for My. = 1/3. 
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heel impact immediately follows the last front impact of the series. 
These cases occur at 


1—-—k 
Q= ae (30) 
and lead. to rebound amplitudes 
2CY, = My + (1 — My)[k" — Mu(1 — k")’] (31) 


In Figs. 5 to 10, these special points are plotted and connected with 
straight dotted lines, which therefore indicate upper limits to rebound. 
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Region I1I—Partial initial heel series for | 


1+k 1+k 
h i — bf Mok FB) ° 


This is a region of critical phasing, and values were determined only 
for the maximum cases, where a front impact just precedes the last 
impact of the partial initial heel series. Here: 


i hie 
+k. — b*)Me 


<Q< 


= face) (32) 
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and: 
2CY, =1— (1 — My) — FB) (14+ [hk — Mo(1 + BL — bP} 
— Mx(1 = My) + ky? {1 sea i” (33) 
+{[k—k” — Mall + kA — ke} 
Region IV—Complete initial heel series. A complete initial heel series 
implies that when the heel has come to rest, the front is still away from 
the stop. When the front finally meets its stop, the situation is identical 
with that just prior to the zero impact except that the energy content of 


the system is lower. The pattern must then repeat with diminished 
amplitudes. For this region we recognize two groups of cases. The first 


0.30 











le 


= 
R 


et 
a tL Medel 


° 
iN 






2 
fs) 
> 
> —, 
eae 


Sos 


RELATIVE MAXIMUM REBOUND AMPLITUDE, 2CY, 
o 
= s 
o 


feos 
ea 


») 
x 
- 








Ee sie eee etd eal dest 


Ory Jo 9 
ele | fe 


6 7 8 9 10 
Prof 


Fig. 9—Relative maximum rebound amplitude for My = 1/6. 
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group is that in which the front velocity is positive at the completion of 
the heel series. In that case 


2(1 + k) Sos (l+k) | 
kL — k) + Mu(l + by kK — k) + Munk + b? 


and 


(1 = My)k? 
205 SS 3 
i alae a, (34) 


For the second group the front velocity is still negative when the heel 
comes to rest from which point on the system acts as a one-degree-of- 
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Fig. 10—Relative maximum rebound amplitude for My» = 1/8. 
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freedom system until the next front impact. The requirement for this 
group is that 


(1 +B) 
> Fa —® + Mall op 


and the maximum rebound is given by 
2CY1 = My —(1 — Mu)[(1 — #) + Mu(1 + 4)’ (35) 


It is to be noted that in the upper part of Group 1 the amplitude in- 
creases with successive heel impacts. This can be explored through the 
use of Equation (22). For simplicity of mapping, however, the limit 
given by Equation (35) has been extended back from the lower boundary 
of Group 2 until it intersects the line marking the first rebound ampli- 
tude of Group 1. 

In Figs. 5 to 10 the respective regions have been identified by means 
of the symbols indicated below: 


Region I from Ppf = 0 to O 
Region II from O to A 
Region III from A to X 
Region IV, Group 1 from X to O 
Region IV, Group 2 from O to Ppf > © 


E. Discussion of Rebound Charis 


Aside from quantitative data contained in Figs. 5 to 10, the following 
general trends are of interest: 

For values of Miz > 4, and the values of & under consideration, most 
of the useful range of Py:f involves critical phasing and the rebound 
maxima are relatively high. 

For values of + < My < 4, consistently controllable rebound ampli- 
tude may be obtained. 

For values of Mi. < + rebound increases again and the structure 
approaches the one-degree-of-freedom case. 


VII. ANALYSIS OF REBOUND PATTERNS-——-THREE-DEGREES-OF-FREEDOM 
SYSTEM 


Rebound pattern analysis as in Parts V and VI has so far not. been 
performed for the three-degree-of-freedom system, partly because of 
complexity, and partly because for the system of Fig. 3 friction at the - 
hinging stop will greatly influence the motion. 
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However, it is felt that the approach and notation of the analysis 
presented here is sufficiently general to allow extension of the rebound 
pattern analysis to the three-degree-of-freedom case. At any rate, with 
the assumption of the magnitudes of frictional forces, the basic equations 
of Part IV may be used to plot any particular case. 


VIII. ARMATURE REBOUND MODEL 


In order to verify the formal analysis presented in Parts III, IV and 
V, a large model of a two-degree-of-freedom system was constructed. 
It consisted essentially of a large bar constrained to move in a plane, 
biased against two stops, and to the ends of which writing pens were 
attached. As rebound conditions were simulated by releasing the bar 
against its stops, chart paper moved at right angles to the bar motion 
and thus produced a record of end displacement versus time. 

By changing spring members and attaching masses to the bar, it was 
possible to vary the mass distribution and the biasing forces. 

The results obtained closely agreed with those suggested by the 
analysis. The maximum rebound amplitudes were generally somewhat 
lower probably due to frictional effects. 


IX. RELAY DESIGN CRITERIA RESULTING FROM ARMATURE REBOUND 
ANALYSIS 


A. Limitation of Analysis 


The assumptions which this analysis is subject to have been described 
under Part II. As applied to relays and probably the majority of mechani- 
cal structures, one assumption is most frequently and severely violated. 
The stops, which have been assumed to be stiff springs associated with a 
definite coefficient of restitution are, In practice, massive bodies which 
dissipate energy through excitation of high frequency modes of vibration. 
Accordingly, the assumption that the stops are at rest is violated, 
particularly if the mechanism is subject to repetitive (pulsing) im- 
pacts and the stop vibration does not decay greatly in the repetition 
period. | 
_ However, mechanisms designed in accordance with this analysis have 

performed well even under moderate pulsing conditions if the sensitive 
phasing region was avoided. In addition, every effort should be made to 
reduce the amount and duration of stop and mounting structure vibra- 
tion by making them stiff, massive, and dissipative, if possible. 
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B. Design Criteria 


1. Type of Armature Structure. 

The selection of the number of degrees-of-freedom for an armature 
structure depends on the expected coefficient of rebound as well as 
practical considerations. 

It can be shown without great difficulty that for very low coefficients 
of rebound the one-degree-of-freedom system is preferable. This is quite 
obvious when one considers the limiting value of k = 0. In this case the 
one-degree-of-freedom system will have no rebound whatsoever, while 
the two-degree-of-freedom system has a heel bounce followed by re- 
bound at the front. The value of k below which the one degree system is 
preferable varies with the mass distribution relative to the stop points, 
being 0.18 for a rectangular plate armature with stops located at its 
edges. 

Experience indicates that k in most practical relays and similar 
mechanical structures varies from 0.3 to 0.6. Hence the two-degree-of- 
freedom system is superior in all practical cases to the solidly hinged 
armature. 

As far as three and higher degree-of-freedom systems are concerned, 
it may be said that generally the greater the number of modes resulting 
in impacts, the quicker the rebound energy can be diverted and dis- 
sipated and the lower theoretical rebound values can be obtained. This 
consideration would favor systems containing many degrees of freedom. 
However, while multi-degree-of-freedom systems can reach very low 
rebound values, their motion (phasing) must be very closely controlled 
or they may prove to be inferior to simpler systems particularly under 
vibratory (pulsing) operation. It is this difficulty which makes it appear 
that the two-degree-of-freedom system offers the best promise with the 
three-degree system also quite promising. By the same reasoning, 
additional spurious rocking modes should be minimized. 

2. Armature Mass. 

The armature mass should be as low as possible. This will minimize 
stop and structure vibration. In addition, in. relay applications light 
armatures tend to increase magnetic ‘‘drag” losses of the armature 
during the release motion. 

3. Stops and Mounting Structure. 

As discussed before, it is desirable to reduce the amount and duration 
of stop and mounting structure vibration. 

4. Coefficient of Restitution. 

The coefficient of restitution should be kept low. Stops having low 
stiffness should, therefore, be avoided. 
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5. Biasing Forces. 

F, should be kept as high as practicable. 

For proper energy loss during impacts, the motion between impacts 
must occur outside the region of the compression, i.e., the armature 
and stop must separate. Therefore, because all practical stops have a 
finite stiffness, the biasing forces (Ff , Ff. , etc.) should produce a static 
deflection less than say, arbitrarily, 5 per cent of the maximum expected 
rebound amplitude. 

6. Design Parameters for Two-Degree-of-Freedom Systems. 

As clearly indicated in Figs. 5 to 10 for the practical range of coef- 
ficients of restitution, most consistently good results are obtained with a 
coupling factor My, = 5 to 4. This factor is most easily adjusted by 
correct placement of the front stop. 

For the above range of My. the force ratio F2/F, should be such as to 
make the product 


F 
Pua > 4 My = } 
>3 My = 4 
>3 My = ¢% 


(Note: For a rectangular armature structure with the stops placed at 
its edges Mi. = 4, Pw = } and force ratios in the neighborhood of 8 are 
desirable.) 
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APPENDIX I 


DERIVATION OF BASIC EQUATIONS OF MOTION THREE-DEGREE-OF-FREEDOM 
SYSTEM 


(1) Free Interval 


The motion of the armature will be described by the displacement at 
the stop points, 2 , X2 , 3. Let m be the mass and RF the radius of gyra- 
tion of the armature about the center of gravity. The latter is located 
by the dimensions 4R, £R, and GR relative to the stop points (Fig. 3). 
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The rotation and displacement of the center of gravity is then 





ae ae fs 
Lp = (t2 — 2X) nar + 23 
peer ee Cee ee (a) 
Ait h 
§= V2 — UW 
~ RA + &) 
From this the kinetic energy may be computed 
T = 4m(a, + a) + ¥mR7O" 
_ a+ +) + i+ G+ 1) + 4G + 6) 
2(& + 1)? (b) 
4 trial — fs — 1) — 2ivb(G + 4) (hi — te) 
2(4 + £)? 


Applying LaGrange’s Equation to the above, the equations of motion 


are obtained: 
d (7) _o? _p 








dt aqr OGr 
Py _ él + Gs + 1) + ale — & — 1) — ae + 6) 
m (a + 6) 
Fy _ filth — fe — 1) + alti + G+ 1 + a0(G + 8) (c) 
m (fA + 6)? 
m G+ 6) | 


The Equations (8) may be solved for #1, #2 , #3 and the results inte~ 
grated, yielding 


1 = = [Cul + Cok: + CisF alt? + tut + sro 
| | 
n= a [CaF + Cools + Coss]? + ctoot + 22 (d) 
| 
} 


3 = es [CaF a Cooks + Cask slé? +- teot + 220 
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where 
Cu= (+1) Cy = Cn = Ab 
Ce = (+1) Ce = Ca = (1 — 48) (3) 
Cu=(@ +1) Cy = Cn = —bh 


X19 , L290 , X30 are the initial velocities, xo , X20 , X30 the initial displace- 
ments for the free interval in question. Interpretation of the analytic 
results is simplified by the introduction of normalization. Let %. be 2; 
just before the ‘‘zero” impact and define 








ee __ 4 _ a 
Y= a7) em OE () Ue 
d(- | 
B, (2) 
_ £m 
T F, 


Dividing Equations (d) by éar yields the normalized equations of 
motion: 


HL oho] +a em 
no flere col] (OS rim() rm 0 


fe (Fs) 
Y= 5 | Cn + Co Gry + Ou | (E a] + iw(t ) + vs 


(2) Impact Interval 


The change of velocity at point “2” due to an impact at “2” is, by 
definition of the coefficient of restitution ‘‘k’’: 


At; = —(1 + kia; (e) 
Since this velocity change occurs as rotation about the conjugate 


point as an instant center of rotation, the impact relationships may 
be written, for an impact at point ‘1’’, 


Ag = —(1 + hi)t 


R 
(; - 6”) 
Ate = = 0+ kay 
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(4-1) Ciz 





=(1+ h)a (ap TD = Ta, (1 + ki)ay 
f, 

At; = —(1 — Ia) —“ 
LR + Pa 
4 

. fils Cis : 

= —(1+) = — — (1 k 
( aq: rt a 1) cr, | + ki)aa 


Similarly it can be shown that impacts at points (2) and (8) follow 
the same pattern. The general impact relations for impact at point “‘7”’ 
are then 

Yion = Yin + Kj¥ien—y (6) 


The first subscript indicates the coordinate, and the second subscript 
indicates the beginning (0) or end (e) of the free interval denoted by the 
third subscript. 


The impact transfer coefficient K;; relating a velocity change at 


point “7”? to an impact at point ‘2’: 


Ky= — 21 + by) (7) 
APPENDIX II 


ANALYSIS OF REBOUND PATTERNS—ONE-DEGREE-OF-FREEDOM SYSTEM 


The equation of motion of this system is 


Yin = $C” + Yront? + Yron | (f) 
where 
Cie 
C= Ue 9 
aa (9) 
pat 
. 


and is measured from the start of the particular interval of free motion 
in question. The impact relationship is 


Yin = —kyYtem—) 


The motion consists of a series of parabolic arcs having periods of 
2Qy10/C in general, or 2/C, 2k/C, 2k’/C, -+-, 2k” "/C. The time elapsed 
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is a convergent series and approaches, for a complete series: 


ike 2 " 2 
Lim = [1 +k +1 aoe et 
a a re Me a ea ea 
The maximum rebound amplitude in any interval is —Yiton/2C. 
The maximum excursion occurs during the first bounce at t’ = 1/C and 
equals —k’/2C. 


(10) 


APPENDIX III 
ANALYSIS OF REBOUND PATTERNS—TWO-DEGREE-OF-FREEDOM SYSTEM 
The equations of motion of this system are 


w= LAL” + Yront! + Yt0n 





Yo = $BU” + Yoont! + Yoon | . 
where A = Cy + Crf 
B= Cyt Cof- | 
_F: | (h) 
P; 
, _ ¢ measured from the start of the par- 
T ticular free interval in question. ] 
A. Complete Front Series 
At the beginning of a front series 
yi = 0 | 
Yi = Yreo | 4) 
Yo = Yre0 
io = tro | 


In a manner analogous to that for the one-degree-of-freedom system 
each front impact reduces 7; to —kyy; . Therefore, after the n** impact, 


Yron = —krYreo 
and the time elapsed in the n** interval is 


2hy , 
fe = ee Yreo : (j) 
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At the heel, from (g), the heel velocity preceding the n** impact is 
Yren—t) = Yoorn—1) + Bl (k) 

The velocity change during the (n — 1) interval is then equal to 
BT,_1. From Equations (6), (7) and (12), the change in velocity during 
the nb impact is —Pio(1 + ki)kf "Yeo - 

The total change of #2 between impacts is then 

Yoon — Yor = BTaa — Pel + kdkr 100 
Similarly in preceding intervals: 


Yoon—1) — Yroon—2y = BTn-2 — Pp(1 + kaki? Yeo 


Yoo2 — Yoo = BT, — Pr + b)khye 
Yoo. — Yreo = — Py(1 + ki)greo 
By addition of the above 


n—1 n—-1 


Yoon — Yreo = Bp, Tm — Pel + ki) 2 kTYr60 
2B : n—1 os ; n—1 2. 
= A Yieo 2 ky — Pll + ki) Yr60 a k 
The summations may be evaluated, yielding 


. _[2Bk — kr ed 
Yon — Yro = E eke Py + ka) =a Yieo (1) 


To evaluate the displacements at the heel, Equation (g) yields 





. 2 
Yon — Yon-1) = Yoorn—1) L'n—1 + 4BTn-1 


Adding these expressions for intervals 0 to n; the total change in y2 is 


n—1 n—-1 


Yrn — Yor = Dy GromT'm + $B DG Tn 
2h — kT) | 


a Atle ee 
2B(ki y Qo) 4. ki") (m) 
ae 1 1 1 
| A?(1 — ky)? 





_ 2Pywky(1 = kt = ki = i a 
Ad = kay? Yeo 
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Expressions for an initial series may be obtained by setting j0 = 1 


Yoo = Yee = O, and, finally, for an initial complete series m — © and 
hence k” — 0, and Equations (1) and (m) become 


a a i I Bla — Pp 
meee AC Ee et A - 


1 2Bk 
Yreo = | — Py»(l + md 


B. Complete Heel Series 


(17) 


For heel series, Equations (1) and (m) may be used by interchanging 
the initial velocities, accelerations, and impact transfer coefficients for 
those relating to heel motion: 


. _ sf sy E ke _ ke = My + ke) 1- a, (n) 
Yion Y1e0 7a ee i Sore as is Y200 
Qk(1 — ke) , 


Yon — Yio = ea Yreoy2co 
(0) 








ies = 2ke*t + ha") _ 2Mavkea(L — ke — ka + ml, 
BU — @) BP» — 9) i 


An initial heel series occurs when the heel strikes first after the “zero” 
impact. The first heel impact then occurs 7; = 2Py»/B(1 + k) after 
the zero impact and the initial conditions are 


You = Pll + ky) 


2AP 
Ya = —ky +A, = ee +h) - kh 


AP. 


Ya = —hkT; + LAT} am (1 + ky) | B 








ae Gaia | 


Substitution of the above into (n) yields 








ihn = —y + 74 mee 2 (1 — jy)" Mull + &)(1 — i) | (p) 
ke 


The corresponding expression for Yion 1s quite Involved. For the special 
case of k = ky = ke 
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APy(1 + k)® Ee (: = Ey _ kl — &”) 





ae B Bo ieok 1_-- iy 
q 
— Mel — R’)(k -— zi 
(1 — k)? 
If the initial series is a complete series, 7 — © and 
= AP,(1 + ky? APy» kA — k) “3 | 
Me BG aE | Be ae 

(22) 


B (+h) 





_1lt+k E= ka — k) 


Yieo = C5 = My2(1 + i) | 


C. Partial Front Series 


The worst rebound occurs when heel and front impacts occur nearly 
simultaneously, with the front hitting first. From Equation (m) for an 
initial front series, this requires that 

B 1— k" 


AP. - k—k (80) 


After n front impacts conditions are given by Equations (14) and 
(19) with y, = ye = O, and 


Dea 
To 


Mays _ 2M 2k Yo 
P 7 E 12 





1-(—- My)(1 — Kk") = }’” 4+. 


This may be solved for yz = Py(1 — k"). The maximum front excur- 
sion now possible is that for a complete series of heel impacts. The 
above value of # in Equation (24) yields 


2CY1 = Me + (1 — My)[k" — My(1 — k"y’ (31) 


D. Partial Heel Series 


- The worst rebound occurs again when heel and front impacts occur 
nearly simultaneously, with the front hitting first. From Equation (9) 
for an initial heel series, this requires that 


B 1 = Art 


ae ee zn - 4k — My 





(32) 
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After n heel impacts 72 = Py(1 + k)k” and from Equations (19) 
and (23) 


T+V 
To 





=1—(1 — My)(l — B) — Mal — My) + bP — 2B") 


= yi — 2Mn(1 + A) "G + Mu(1 + &)?k 


This may be solved for 7 = k — My(1 + k)(1 — k”) and after the 
front impact immediately following: 


Y. = Pll + k)k” = Ppl + k)[k = My(1 + k)(1 = k")] 


The maximum front excursion now possible is that for a complete 
series of heel impacts. The above value for #2 in Equation (24) yields 


20%, =1-—0=My)d =e) (1+ -— me + OHO = 27} 
— My(1 — My)(1 + kb)? {1 — hk (33) 
+ [k — kh” — Mp(l + AG — k)P} 


APPENDIX IV 


SUMMARY OF NOTATION 


A= Cu + Ce + Cof 

B= Cy = Cex + Cuof 
Cie 

C= Cy Cn 


Cu=1+h 2 Cy Ca = [1 — G6] 
Co =14+ 6 C3 = Cu = 4G 


I 


C33 =] + fs Cos = C2 — folo 


F, = front tensioning force 
F, = heel tensioning force 
ky = eeeticient of restitution at vertical front stop 


ke = coefficient of restitution at vertical heel stop 
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Ci. 
Kj: a en (1 + ki) 
Ci 


{ik = vertical front stop location relative to c.g. 
ff = vertical heel stop location relative to c.g. 


fk = horizontal heel stop location relative to c.g. 








m = mass of armature 
2 
Mu = Oo, 
P. 

ct Map Sei 

1+ Pruf tb kh 
Rk = radius of gyration of armature about center of gravity 
ae Lam 

Py 
t = time 
fos 

- : 
t; = basic period of front after “zero” impact 
tz = basic period of heel after “zero” impact 
Tn = duration of n** free interval 
2%, = vertical front displacement 
x2 = vertical heel displacement 
x3 = horizontal heel displacement 
ta = front velocity just prior to ‘‘zero” impact 
pa 


Let 


Y; = greatest excursion (rebound) of front 
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is no arcing at the surfaces; but with continued operation and accompanying 
arcing, the activating material is burned away, and the surface returns to the 
inactive condition if no activating vapor is supplied. 

Are voltages, which are independent of current and of ambient gas, as far as 
tested, have been measured for a number of metals and for carbon; the arc 
voltage for carbon is quite erratic in the range between 20 and 380 volts, but for 
each of a number of metals the arc voltage is steady. 

Arcing at noble metal surfaces, similar to that induced by carbonaceous 
material from organic vapors, can be produced also by insulating particles or 
insulating films. The active condition gradually disappears with continued arc- 
ing, unless there is a steady supply of insulating material to the surface. 

The minimum arc current has been measured to be 0.6 amp for active silver’ 
and for carbon, and 0.03 amp for inactive silver. These are the currents at 
which an established arc is extinguished. 


Iron-Stlicon Alloys Heat Treated in a Magnetic Field.* M. Gorrrz!, 
Jl. Appl. Phys., 22, pp. 964-965, July, 1951. 
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429-434, July, 1951. 
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silica buffers; and waves generated by quartz crystals at the ends of the assem- 
bly are transmitted simultaneously through the specimen in both directions. 
Comparison of transmitted and reflected components by means of a special 
balancing circuit provides information from which the complex propagation 
constant can be calculated, and hence dynamic rigidities and viscosities. 

Illustrative data for polyethylene and Nylon are given. 
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J. Opt. Soc. Am., 41, pp. 454-456, July, 1951. 

The semi-metals germanium and silicon have high transparency and high 
refractive indices over a wide range of infrared wavelengths and are stable to 
normal atmospheres. Their relative indices and dispersions make achromatic 
combinations possible; and designs are given for axially corrected doublets of 
relative apertures f{:2 and f:1. The optical homogeneity of the materials is dis- 
cussed: compositional variations are not considered an optical hazard, but there 
is evidence of structural imperfections in some specimens whose effect on optical 
properties remains to be evaluated. 


On the Motion of Gaseous Ions in a Strong Electric Field. I.* G. H. 
WannieEr!. Phys. Rev., 83, pp. 281-289, July 15, 1951. 

This paper applies the Boltzmann method of gaseous kinetics to the prob- 
lem of positive ions moving through a gas under the influence of a static, uni- 
form electric field. The ion density is assumed to be vanishingly low, but the 
field is taken to be strong; that is, the energy which it imparts to the ions is 
not assumed negligible in comparison to thermal energy. Attention is focused 
upon the computation of velocity averages, and the drift velocity in particular, 
rather than a complete knowledge of the entire velocity distribution. It is shown 
in Sections C and E that the problem so formulated is completely soluble if 
the mean free time between collisions of ions and molecules is a constant; this 
is the case for the so-called polarization force between ions and molecules which 
predominates over other forces at low temperature. A method for obtaining 
averages to any desired accuracy in the general] case is developed in Section D. 
The method is applied to the hard sphere model for the high field range and 
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mass ratio 1. An application of the resulting formula (48) to experimental mate- 
rial has been published earlier. 


Evidence for the Noncubic High Temperature Phase of BaTiO3;. E. A. 
Woobp!. Letter to the Editor. J. Chem. Phys., 19, p. 976, July, 1951. 


Seven-League Oscillator.* F. B. ANpERSON!. Proc. I.R.E., 39, pp. 881- 
890, Aug., 1951. 

A bridge-type RC oscillator is described which is continuously adjustable 
over a frequency range of 20 cps to 3 mc in one sweep of a two-gang linear 
potentiometer control. Tracking requirements of the two-gang control are not 
severe. The output is available in four phases, and the frequency is an approxi- 
mately logarithmic function of the linear potentiometer setting. Practical limits 
of the frequency range are tentatively 0.01 cps and 10 me. Accuracy of setting 
of the order of one per cent is attainable with ordinary components. Frequency 
stability is of the order of 2 per cent per db of tube gain variation. 


Semi-Conductor Surface Phenomena.* W. H. Brarratn!. Semt-Con-. 
ducting Materials, H. K. Heniscu, ed., pp. 37-46. Proceedings of a con- 
ference held at the University of Reading (July 10-15, 1950) London, 
Butterworths, 1951. 

Developments in the understanding and interpretation of phenomena occur- 
ring at the surface of a semi-conductor are reviewed. The development starts 
with the Mott-Schottky theory of the space charge layer. Bardeen’s concept of 
a space charge layer due to ‘surface states’ explained the independence of recti- 
fication on contact potential and Meyerhoff’s small values for the contact poten- 
tial between n- and p-type silicon. Shockley and Brattain observed that this 
contact potential difference increased with impurity concentration. Illumination 
of a silicon surface produced hole and electron pairs in the space charge layer. 
The potential of the surface changed until the photocurrent was balanced by a 
conduction current. The relation between photo-current and potential change 
was of the same form as a forward characteristic for a rectifying contact. In an 
experiment similar to Becquerel’s, using water as an electrolyte, the surface 
may be biased in the reverse direction. When so biased the response to modu- 
lated light gives the differential resistance of the space charge layer. This re- 
sistance increases rapidly with reverse bias and the time constant of the layer 
increases, both agreeing qualitatively with theory. This experimental method 
of measuring changes in surface potential caused by illumination permits deter- 
mination of the properties of the space charge layer at the free surface of a semi- 
conductor, 


The Calculation of Traveling-Wave-Tube Gain.* C. C. Curtmr!. Proc. 
T.R.E., 39, pp. 914-917, Aug., 1951. 
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Essential information for claculation of traveling-wave-tube gain is sum- 
marized and condensed in this paper. The important relations are documented, 
presented in a concise form for simplified computation, and developed as a 
nomograph. The conclusions have been found to be in agreement with measure- 
ments on six different tube designs. 


Thermal Variation of Young’s Modulus in some Fe-Ni-Mo Alloys.* 
M. E. Fine! and W. C. Exuis!. References. J/. Metals, 3, pp. 761-764, 
Sept., 1951. 


A Broad-Band Transcontinental Radio Relay System.* T. J. GRIESER! 
and A. C. Prrrrson!. lec. King., 70, pp. 810-815, Sept., 1951. 

Spanning the continent from coast to coast, this microwave relay system 
provides six channels, each of which can carry one television circuit or hundreds 
of telephone circuits. Some features of this vast network are descirbed. 


An Improved Telephone Set.* A. H. Ineuts' and W. L. Turrne.w, 
Elec. Eng., T0, pp. 770-775, Sept., 1951. 

The familiar telephone set has undergone numerous changes which will pro- 
vide better service at lower cost than do present models. Increased transmitting 
and receiving gain, better sidetone control, broader frequency response, faster 
dialing, simple ringing control, and a trim appearance are some of the features 
of the new design. 


The Institutes for Basic Research—Their Contribution to National 
Strength. M. J. Keuiy', pp. 11-23. Applied Research is Not Enough, 
(booklet). Addresses at the Dedication of the Institutes for Basic Re- 
search, The University of Chicago, May 16, 1951. 


The Crystal Clock. W. A. Marrison!. Science Marches On, JAMES 
Sroxuey, ed., N. Y., Ives Washburn, Inc., 1951, pp. 303-309. 


Observations of Zener Current in Germanium p-n Junctions. K. B. 
McAree!, E. J. RypEr’, W. SHockxiey!, and M. Sparks'. Letter to the 
Editor. Phys. Rev., 83, pp. 650-651, Aug. 1, 1951. 


Experimental Radio-Telephone Service for Train Passengers.* N. 
Monk). Proc. I.R.E., 39, pp. 873-881, Aug., 1951. 

Experimental public radio-telephone service for train passengers was in- 
augurated by the Bell Telephone System several years ago. Initial installations 
operated in conjunction with a series of urban mobile base stations. More re- 
cently, highway mobile systems have been used for this service, and this paper 
describes a typical train installation operating through a highway channel. All 
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of these early systems utilized an attendant on the train. The cost of providing 
an attendant has, in some cases, been found excessive. Consequently, experi- 
ments have been initiated in which a coin box is used on the train. The arrange- 
ments for this purpose are also described. 


The Magneto-Resistance Effect in Oriented Single Crystals of Ger- 
manium.* G. L. Pearson! and H. Sunu!. Phys. Rev., 83, pp. 768-776, 
Aug. 15, 1951. 

This paper describes an extensive study of the magneto-resistance effect in 
germanium as a function of crystal orientation. Experimental measurements 
establish the constants involved in the dependence of the effect on orientation 
of magnetic field and electric relative to the crystal axes. The measurements 
are internally consistent with existing phenomenological theory based on cubic 
crystal symmetry, in which terms involving the magnetic field to higher than 
the second order are neglected. It is shown that such deviations as do occur 
arise from higher terms in the field, since an extension of the phenomenological 
theory to the fourth order predicts their symmetry. Relations are established 
between the experimentally observed phenomenological constants and those 
constants appearing in existing magneto-resistance electronic theories. It is con- 
cluded that no electronic theory yet worked out is entirely consistent with 
experiment. The present electronic theories are special cases of a very general 
theory recently proposed by Shockley, and it is possible that agreement can 
be obtained as soon as the computational difficulties of the latter theory are 
overcome. 


New Phenomena of Electronic Conduction in Semi-Conductors. W. 
SHockuey!. Semi-Conducting Materials, H. K. Heniscu, ed., pp. 26-36. 
Proceedings of a conference held at the University of Reading (July 
10-15, 1950) London, Butterworths, 1951. 

The semi-conductors silicon and germanium may be discussed as insulators 
the electronic structure of which is disturbed. Excess electrons, which act as 
negatively charged current carriers, may be present as may be ‘holes’ or places 
where electrons are missing from the valence-bond structure. Holes act as 
positively charged current :carriers. In ordinary electronic conduction the flow 
of current carriers is substantially incompressible so that the density of carriers 
remains constant. When a new transistor phenomenon known as ‘carrier injec- 
tion’ occurs, however, the total density of holes and electrons may be greatly 
increased and this modulation of the electronic structure may be used both for 
scientific investigation and for practical amplification. In particular, carriers 
may be injected at a predetermined time and place into a known uniform elec- 
tric field and their transit time to another place accurately timed by detecting 
their arrival with a “collector” point. Drift velocities and mobilities may be 
measured precisely in this way with a directness unattainable by pre-transistor 
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methods involving conductivity and Hall effect. These new experiments and 
their related theories may furnish a foundation for a new engineering science of 
transistor electronics. 


Domain Patterns on Nickel.* H. J. Wiuurams! and J. G. WALKER}. 
Phys. Rev., 83, pp. 634-636, Aug. 1, 1951. 

Domain patterns have been observed on two single crystals of nickel cut in 
the form of hollow parallelograms. The length of the sides were parallel to the 
(111) directions in one specimen and to the (110) directions in the other. The 
crystals show domain structures with the three types of domain boundaries 
which are to be expected from a material having the directions of easy magneti- 
zation along the (111) directions. Domain boundary movement under the influ- 
ence of an applied magnetic field was observed. 


Polymorphism in Potassium Niobate, Sodium Niobate, and Other ABO; 
Compounds.* EK. A. Woop’. Bibliography. Acta Cryst., 4, pp. 353-362, 
July, 1951. 

The first part of this paper presents the results of optical and X-ray studies 
of the perovskite-type crystals, potassium niobate and sodium niobate. Potas- 
sium niobate is orthorhombic at room temperature, changing to tetragonal at 
about 225°C. and cubic near 485°C. Sodium niobate is orthorhombic at room 
temperature, changing to tetragonal at about 370°C. and to cubic at about 
640°C. 

The second part of the paper discusses relations among the structures of the 
ABO; compounds. 


Subjective Sharpness of Additive Color Pictures.* M. W. Baupwin!. 
Proc. I.R.E., 39, pp. 1173-1176, Oct., 1951. 

This is a report on the first numerical results to come from a laboratory 
experiment on the subjective sharpness of additive three-color pictures. The 
sharpness factor is isolated by using out-of-focus projection (of slides) instead 
of actual television transmission. 

An observer’s acuity for defocus is pauses for the green component and least 
for the blue component, in an additive three-color picture. When the same 
picture is reproduced in monochrome (white, red, green, or blue) at the same 
brightness, the observer’s acuity for defocus is equal to that found for the 
green component. | 


Frequency-Modulation Terminal Equipment for the Transcontinental 
Relay System.* J. G. CHarrse! and J. B. Maaato!. Elec. E'ng., 70, pp. 
880-883, Oct., 1951. 

To meet the exacting requirements of the new transcontinental microwave 
relay system, specially designed frequency-modulation terminal equipment was 
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constructed. The terminal transmitter converts either message or television 
signals to a frequency-modulated signal centered on 70 me and the terminal 
frequency-modulation receiver recovers these signals, thus providing a link 
between the relay system and other telephone facilities. 


Observer Reaction to Low-Frequency Interference in Television Pictures.* 
A. D. Fowuter'. Proc. I.R.E., 39, pp. 1332-1336, Oct., 1951. 

This paper presents results of tests to determine how much low-frequency 
interference can be tolerated in black-and-white television pictures. Various 
levels of single low-frequency interference were superimposed on a locally trans- 
mitted television picture. Observers viewed the picture and rate the disturbing 
effect of each level of the interference. Ratings were made in terms of preworded 
comments ranging from “not perceptible” to “unusable.” Interfering frequencies 
from 48 to 90 cycles per second were employed. 

Just visible interference appears as a flicker. The rate of flicker is the differ- 
ence between interfering and 60-cycle field frequencies. The most disturbing 
interference produced a flicker rate of 5 or 6 cycles per second. To be tolerated, 
peak-to-peak amplitude of this interference had to be 54 db weaker than the 
peak-to-peak amplitude of the television signal (including synchronizing pulse). 
For flicker rates of 0.5 and 12 cycles per second, the amount of interference 
which could be tolerated was larger by 14 and 3 db, respectively. 


Arcing at Electrical Contacts on Closure. Part II. The Initiation of an 
Arc.* L. H. Germer'. Jl. Appl. Phys., 22, pp. 11383-1139, Sept., 1951. 

The capacity of the plates of an oscilloscope charged to 35 or 40 volts is dis- 
charged repeatedly by approaching electrodes of carbon, active silver, and inac- 
tive silver. Facts about the discharges, which are arcs of very short duration, 
are inferred from resulting open circuit potentials and calculated electrode 
separations. 

The separation at the first arc varies in different experiments but corresponds 
on the average to a nominal electric field of 0.6 X 10° volts/em for carbon or . 
active silver and to 2 X 105 volts/em for inactive silver. Each arc is initiated 
by a very small number of field emission electrons. The hypothesis that a single 
electron may perhaps be sufficient is consistent with observations at later stages 
of each closure when the electrodes are closer and the field much higher. 

The earlier observation, that the potential across a short arc is constant and 
independent of current, is not true if the are time is sufficiently short. For 
active silver a time comparable with 2 * 1078 sec is required to establish the 
steady arc voltage characteristic of later stages of arcs which last longer than 
this. The initial time during which the potential is decreasing toward its final 
steady value is 100 times the transit time of a silver ion across the gap. 
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Computation of Control Limits for p-Charts When the Samples Vary in 
Size. H. L. Jonus®. Ind. Quality Control, 8, pp. 26-27, Sept., 1951. 


The Design of Switching Circuits. W. Kurster!, A. E. Rircuin', and 
8S. H. Wasuspurn'. NV. Y., Van Nostrand, 1951. 556 pp. (Bell Telephone 
Laboratories Series). 

This is the first published textbook in its field. It presents, first, the funda- 
mental design principles of switching circuits composed of discrete-valued 
switching elements. Most of the discussion concerns two-valued elements, with 
greatest emphasis placed on electromagnetic relays. Chapters cover basic circuit 
paths, the logical interpretation of requirements, and the techniques of con- 
structing networks to fulfill these requirements. The symbolic methods of 
Boolean algebra and its application-to the design of combinational and sequen- 
tial circuits is covered. Later chapters cover various unifunctional circuits such 
as selecting, connecting, translating, counting, and lockout. Final chapters dis- 
cuss methods of synthesising unifunctional circuit building blocks into larger 
circuits and systems. 


Measurement of the Elastic Constants of Silicon Single Crystals and 
Their Thermal Coefficients. H. J. McSximin!, W. L. Bonp!, E. BuEHLER!, 
and G. K. Trau'. Letter to the Editor. Phys. Rev., 83, p. 1080, Sept. 1, 
1951. 

Interest in the properties of silicon single crystals arising from their use as 
semiconductors has led us to make measurements of the elastic constants of 
two single crystals. Measurements of velocities of propagation for both shear 
and longitudinal waves were made in the crystals as described in a recent paper 
by McSkimin. Frequencies in the range 8-12 mc/sec were used. 

The three independent elastic constants were evaluated, a density of 2. 331 
(measured by pycrometer) being used. Data and formulas used are summarized 
in Table I. Two crystals were measured—as indicated—with data obtained 
from the larger one being used to determine the elastic constants. Check meas- 
urements were made for the smaller crystal; and despite the less accurate ‘‘pulse 
overlap” technique used for two of the measurements, velocity agreement to 
within 0.15 per cent was obtained. 

Both crystals were of a high degree of crystalline perfection as shown by 
etching and X-ray tests. 


Domain Wall Relaxation in Nickel. W. P. Mason!. Letter to the 
Editor. Phys. Rev., 83, pp. 683-684, Aug. 1, 1951. 


Phase Transitions in Ferroelectrics. B. Matrutias!. National Research 
Council, Comm. on Solids. Phase Transformations in Solids. Ed. by 
R. Smoluchowski, J. E. Mayer, W. A. Weyl. N. Y., Wiley, 1951. 660 pp. 
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Under the name ferroelectrics are classified those materials which exhibit 
dielectric anomalies phenomenologically similar to the magnetic behavior of 
the ferromagnetics. Perhaps it would have been more logical to use the term 
Rochelle-electrics, thus emphasizing the similarity in the dielectric behavior to 
that of Rochelle salt. 

In this paper the known ferroelectrics are listed first, and then there follows 
a discussion of the various theories which have been created to explain them. 


Data on Random-Noise Requirements for Theater Television.* P. 
Mertz!. Jl. S.M.P.T.E., 57, pp. 89-107, Aug., 1951. 

Provisional evaluation of permissible random noise for theater television is 
considered from several sources of information. These cover broadcast television 
experience and the graininess in motion picture film; the requirements deduced 
from the various sources generally agree. For broadcast television, a frequency 
weighting and limit on weighted noise power have been used. The finer picture 
detail of theater television presumes a lower permissible random noise. Changes 
in weighting curve are discussed. A limit figure of noise is suggested, which is 
comparable to graininess effects in motion pictures, though slightly more severe 
than present published performance on camera, tubes. 


A Spatial Harmonic Traveling-Wave Amplifier for Six Millimeters 
Wavelength.* 8. Miutman’. Proc. I.R.E., 39, pp. 10385-1048, Sept., 1951. 

This paper describes a traveling-wave amplifier in which the electron beam 
interacts with a spatial harmonic of an electromagnetic wave propagating along 
an array of resonator slots. The result is a considerable reduction in operating 
beam voltage for a given physical separation of the circuit elements. This type 
of amplifier operating at about 1,200 volts has yielded net power gains of about 
18 db in the 6-mm wavelength region. A magnetic field of about 1,600 gauss is 
sufficient for proper beam focusing. Aside from small variations of gain with 
frequency that is caused by internal reflections, the bandwidth is of the order 
of 3 per cent. 


Form of Transient Currents in Townsend Discharges with Metastables.* 
J. P. Mounar’. Phys. Rev., 88, pp. 933-940, Sept. 1, 1951. 

The form of the current is calculated for a Townsend discharge stimulated 
by a pulsed light beam, with particular reference to the current component 
initiated by metastable effects. The calculation is directed particularly to the 
development of methods for quantitative interpretation of current patterns 
observed experimentally. 


Studies of y-Processes of Electron Emission. Employing Pulsed Town- 
send Discharges on a Millisecond Time Scale.* J. P. Mounar!. Phys. 
Rev., 83, pp. 940-952, Sept. 1, 1951. 
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The relative amounts of electron emission from the cathode in a Townsend 
discharge caused by ions, photons, and metastables have been studied experi- 
mentally for several cathodes in argon, using pulsed-light stimulation of the 
discharge. The current initiated by metastables exhibits a slow build-up and 
decay, thus permitting easy separation from the faster rising effects of gas 
ionization and electron emission by photons and ions. Time constant studies of 
the slow component yielded a diffusion constant for metastable argon atoms of 
45 cm? sec! at one millimeter pressure. The efficiencies of electron emission by 
metastables and ions was found to be closely the same, while the quantum 
yield for photon emission was found to be generally smaller. 


Electrical Properties of aFe,03 and aFe,03 Containing Titanium.* F. J. 
Morin}. Phys. Rev., 83, pp. 1005-1010, Sept. 1, 1951. 

Electrical conductivity, Hall effect, and Seebeck effect have been measured 
on two sets of polycrystalline samples of aFe,O; and ale.O; containing from 
0.05 to 1.0 atomic per cent titanium (n-type impurity). One set of samples con- 
tained 0.6 atomic per cent excess of iron (n-type impurity), the second set con- 
tained 0.6 atomic per cent deficiency of iron (p-type impurity). 

The conductivity of pure aFe.0; is independent of this amount of stoichio- 
metric deviation. The slope of the log conductivity vs reciprocal temperature 
plot is 1.17 ev and the intercept at 1/T = 0 is 2.1 X 10! ohm cm. Room 
temperature conductivity varies from —10-* ohm=! cm@ (extrapolated) for 
pure aFe,0; to 0.38 ohm- cm- for aFe.O; containing 1.0 atomic per cent 
titanium. 

The measured Hall voltages seem to result entirely from magnetization of 
the samples, which are weakly ferromagnetic, and disappear above the ferro- 
magnetic Curie temperature. 

The temperature variations of the Fermi level are determined from Seebeck 
data. The temperature variations of carrier concentration are determined from 
Fermi level and of mobility from carrier concentration and conductivity for 
some samples. Carrier concentration results indicate that each added titanium 
ion donates approximately one electron to the conduction process. Mobilities 
are found to be less than 2.0 em?/volt sec, suggesting that conduction involves 
electrons in the d level of iron. 


Acceptance Inspection of Purchased Material.* J. HE. Patmrr’ and E. 
G. D. Paterson!. Ind. Quality Control, 8, pp. 15-19, Sept., 1951. 

This paper describes some of the principles and procedures employed in the 
inspection of purchased material in the form of components or finished products. 
The authors’ experience has been largely with procedures used in the Bell Sys- 
tem, and the illustrations have therefore been drawn from this source. It is felt, 
however, that considerations leading to the choice of specific inspection tech- 
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niques will be generally applicable even though the number and volume of items 
purchased and the number of suppliers involved may in some cases differ widely. 
In this presentation, stress has been placed on a discussion of the broader gauge 
factors underlying the engineering planning of inspection procedures rather than 
on specific sampling and control techniques. 


Analysis of Audio-F'requency Atmospherics.* R. KX. Porrer'. Proc 
I.R.L., 39, pp. 1067-1069, Sept., 1951. 

Sound portrayal techniques used in studies of speech and noise reveal the 
structure of atmospheric disturbances well known to long-wave radio and ocean- 
cable engineers as ‘‘whistlers,” “‘swishes,’’ and tweeks.”’ It is suggested that 
renewed investigation of these effects, using modern analyzing tools, might yield 
information of considerable scientific interest. 


Reflection of Electromagnetic Waves from Slightly Rough Surfaces.* 
S. O. Ricr!. Communications on Pure and Applied Math., 4, pp. 351-378, 
Aug., 1951. 


Color Television and Colorimetry.* W. T. AN TEIN GE Proc. I.RE., 
39, pp. 11385-1172, Oct., 1951. 

The high lights of the hietow of color measurement and of color photography 
are reviewed. Following this introduction, the principles of modern three-color 
colorimetry are developed from a hypothetical experiment in color matching. 
The conventional theory of “perfect color reproduction” by color television is 
built up from colorimetric background. Some of the difficulties to be expected 
in applying colorimetry to color television are brought out. 

Finally, there is some discussion which tends to show that colorimetry may 
not be a sufficiently powerful tool to provide answers to all of the questions which 
will arise in the reproduction of scenes in color by television. The advantage of 
colorimetry as a background is indicated, however. 
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